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1.  Introduction 


High  temperature  behavior  of  metal  matrix  composites  (MMCs)  has  been  reviewed  recently 
by  Taya  and  his  co-workers  [1-3],  where  creep  behavior  and  thermal  cycling  resistance  of  MMCs 
are  discussed.  Creep  resistance  of  a  MMC  is  superior  to  that  of  the  unreinforced  metal,  while 
thermal  cycling  resistance  of  a  MMC  is  normally  degraded  more  than  that  of  the  unreinforced 
metal.  High  temperature  behavior  of  a  MMC  includes  that  of  the  MMC  at  room  temperature  which 
underwent  a  high  temperature  exposure.  The  key  mechanism  of  a  MMC  that  is  subjected  to 
excursion  to  and  from  high  temperature  environment  is  the  relaxation  of  the  otherwise  high  stresses 
induced  by  the  difference  in  the  coefficient  of  thermal  expansion  (CTE)  between  the  metal  matrix 
and  ceramic  or  refractory  filler. 

This  final  report  will  discuss  a  number  of  subjects  related  to  high  temperature  behavior  of 
MMCs.  The  experimental  work  includes  design  of  creep/thermal  cycler,  thermal  cycling  resistance 
of  SCS-6/Ti-15-3  matrix  composite,  dimensional  change  of  SiCp/Al  composite,  all  of  which  are 
discussed  in  chapter  2.  The  analytical  study  consists  of  four  subjects:  dislocation  punching  from 
ceramic/metal  interfaces,  chapter  3,  dimensional  change  of  MMCs,  chapter  4,  stress  relaxation  by 
plastic  flow,  interfacial  sliding  and  diffusion  in  an  inclusion  bearing  material,  chapter  5,  and 
dislocation  punching  in  functionally  graded  materials,  chapter  6.  Finally  conclusion  is  given  in 
chapter  7,  followed  by  list  of  publications  in  chapter  8. 


References 

1.  M.  Taya,  1991,  "Creep  Behavior  of  Metal  Matrix  Composites,"  in  Metal  Matrix  Composites: 
Mechanisms  and  Properties,  eds  by  R.  K.  Everett  and  R.  J.  Arsenault,  Academic  Press,  pp.  189- 
215. 

2.  M.  Taya,  1993,  "Mechanical  Behavior  of  Metal  Matrix  Composites  at  High  Temperatures,"  in 
Mechanical  Properties  of  Metallic  Composites,  ed.  S.  Ochiai,  Marcel  Dekker  Inc.,  1993,  pp.  553- 
611. 

3.  R.  J.  Arsenault,  S.  G.  Fishman  and  M.  Taya,  1994,  "Deformation  and  Fracture  Behavior  of 
Metal-Ceramic  Matrix  Composites,"  in  Progress  of  Mater.  Sci.,  Vol.  38,  eds  J.  W.  Christian  and 
T.  B.  Massalski,  pp.  1-157. 
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2.  Experimental  Study 


We  have  performed  several  experimental  works,  (1)  design  of  a  creep/thermal  cycler,  (2) 
thermal  cycling  resistance  of  SCS-6  fiber/Ti- 15-3  matrix  composite,  (3)  measurement  of 
dimensional  change  of  SiCp/Al  composite  subjected  to  combined  creep/thermal  cycling  using  the 
equipment  designed  in  (1),  and  (4)  fracture  toughness  of  PSZ-SS/PSZ  composite-PSZ  laminate. 
The  analytical  model  are  also  developed  to  explain  the  experimental  results  of  (2)~(4),  some  of 
which  are  detailed  in  chapters  3  and  4. 


2,1  Design  of  a  Creep/Thermal  Cycler  for  High  Temperature  Capability* 

We  reported  in  the  first  year  annual  report  a  summary  of  designing  a  creep/thermal  cycler 
with  low  to  intermediate  temperature  capability  (Tmax  up  to  400°C)  where  a  capability  of  a  constant 
stress  loading  was  emphasized.  For  the  second  year,  an  emphasis  is  placed  upon  designing  a  new 
creep/thermal  cycler  with  high  temperature  capability  (Tmax  up  to  1 10CTC)  where  a  constant  load  is 
used.  This  is  because  for  high  temperature  MMCs,  dimensional  change  (or  accumulated  strain)  is 
expected  to  be  small  (less  than  a  few  percent),  thus  justifying  use  of  constant  load  device  which  is 
much  easier  to  construct. 

This  creep/thermal  cycler  is  designed  to  apply  temperature  fluctuations  up  to  1 100°C  under  a 
dead  load.  A  constant  temperature  furnace  travels  up  and  down  while  the  specimen  is  subjected  to  a 
tensile  load.  The  furnace  is  shuttled  from  one  station  to  the  other  with  wire  cables  driven  by  a  1/4- 
horsepower  12-V  AC  electric  motor.  The  station  hold  time  can  be  set  independently  in  2-minute 
increments  up  to  a  maximum  of  32  minutes,  the  construction  of  the  creep/thermal  cycler  is  half 
done  and  will  be  finished  in  a  few  months. 

Test  Stand 


The  test  stand  is  necessary  for  applying  a  load  to  a  test  specimen  for  an  extended  period  of 
time  under  thermal  cycling.  This  structure  can  be  divided  into  two  main  sets;  load  train  and  load 
frame  which  consists  of  the  front  and  rear  frame.  The  front  frame  is  102  mm  high  by  52  mm  wide 
by  6.35  mm  wall  thickness  A500  grade  steel  rectangular  tube,  while  the  rear  frame  is  L50.8mm  x 
50.8mm  x  6.35mm  steel  angles.  Fig.  2.1  shows  the  creep-thermal  cycling  system  schematically. 

*The  results  of  this  work  were  published  in  J.  Experimental  Mechanics,  Vol.  35, 1995,  pp.  66-70, 
"Creep  and  Thermal  Cycling  Fixture  Design  for  Metal  Matrix  Composites,"  by  Taya,  M.,  Lee,  J. 
K.,  Lui,  M.,  and  Armstrong,  W.  D.. 
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Furnace 


The  machine  is  equipped  with  one  300  mm  high  by  450  mm  diameter  resistance  element 
furnace  as  shown  in  Fig.  2.2.  The  furnace  is  mounted  on  a  moving  frame  which  slides  on  guides 
attached  to  the  front  frame.  The  furnace  is  powered  by  eight  356  mm  long  by  25.4  mm  diameter 
silicon  carbide  heating  elements.  Each  heating  element  has  a  102  mm  hot  section.  The  maximum 
steady  operating  temperature  of  the  furnace  can  go  up  to  1100°C.  The  resistivity  of  silicon  carbide 
is  so  low  that  a  transformer  is  required  in  most  cases.  Here  two  sets  of  transformers  are  used  in 
parallel.  As  it  is  heated,  the  resistivity  of  silicon  carbide  decreases  first(negative  temperature 
coefficient  of  resistance)  and  then  increases  (positive  temperature  coefficient  of  resistance).  The 
maximum  recommended  operating  temperature  of  the  heating  elements  is  MOOT.  To  prevent  the 
heat  convection  and  radiation  from  the  top  and  bottom  sides  of  the  furnace,  two  hole  caps  are 
installed. 

The  temperature  of  furnace  is  controlled  by  a  Barber  Coleman  560  controller  and  a  phase  SCR( 
Silicon  Controlled  Rectifier)  power  switch.  The  SCR  modulates  small  increments  of  power  to  the 
load,  unlike  conventional  mechanical  control,  and  eliminates  the  overshoot  and  undershoot 
associated  with  contactor  control.  In  this  design,  the  type  R  thermocouple  is  used  to  monitor  and 
control  the  temperature  of  the  furnace.  Ordinarily  the  controllers  are  powered  by  the  same  current 
as  the  furnace. 

The  external  jacket  of  the  furnace  is  made  of  stainless  steel.  Inside  the  jacket  the  furnace  is 
insulated  with  450  mm  diameter  disks  of  high  temperature  Cerachem  insulating  blanket  and 
Fiberfrax  insulating  GH  board  which  can  stand  vibration,  mechanical  stresses  and  strong  erosive 
forces.  The  maximum  temperature  at  which  the  insulating  blanket  can  be  used  is  1400T. 

Furnace  Transport  System 

An  arrangement  of  the  furnace  transport  system  is  also  shown  in  Fig.  2.1.  The  furnace 
moves  vertically  to  place  its  hot  section  over  a  specimen  while  the  specimen  is  kept  stationary.  The 
furnace  hold  time  at  each  section  can  be  set  independently  by  a  hold  time  setter.  This  results  in 
thermal  cycling  of  a  specimen.  The  furnace  is  installed  on  a  moving  frame  which  slides  on  steel 
rods  attached  to  the  front  frame.  The  moving  frame  is  driven  by  a  1/4  horsepower  AC  electric 
motor  with  two  take-up  reels.  The  take-up  reels  wind  3.2  mm  diameter  swaged  steel  cables  which 
are  routed  through  pulleys  attached  to  the  top  of  the  front  frame.  Ten  seconds  are  required  to  move 
a  furnace  from  the  hold  station  to  the  specimen  station. 
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Fig.  2. 1  Creep/thermal  cycler  with  higher  temperature  capability 
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Fig.  2.2  Furnace 


Fig.  2.3  Geometry  of  high  temperature 
MMC  specimen 


Load  Train 


The  load  train  is  a  set  of  parts,  connecting  the  dead  weight  with  the  specimen.  The  load  train 
must  sustain  a  large  tensile  load  while  being  subjected  to  the  time- temperature  conditions  imposed 
by  the  moving  furnace.  The  load  train  is  composed  of  weight  pan,  pulley,  chain,  two  sprockets, 
one  coupler,  two  long  bars,  two  inconel  grips,  and  a  test  specimen.  The  chain  is  made  of  stainless 
steel  which  can  sustain  200  °C  temperature.  The  chain  connects  one  side  of  tensile  bar  with  the 
sprocket  attached  to  the  pulley  to  transmit  the  load  from  the  rear  to  the  front  frames.  The  load  is 
applied  to  the  specimen  through  a  pulley  (radius  rj),  a  sprocket  (radius  r2),  a  chain,  and  a  rod. 
The  load  magnification  factor  due  to  the  cam  and  sprocket  is  rj/  r2.  In  this  design,  the  radius  of  the 
pulley  is  5  times  larger  than  that  of  the  sprocket.  The  coupler  is  made  of  mild  steel  and  has  a  one 
inch  horizontal  hole  through  which  a  hardened  steel  pin  runs.  The  steel  rods  are  connected  with 
high  temperature  inconel  grips  which  are  tapped  and  attached  to  the  specimen  holder. 

Specimen 

Target  metal  matrix  composite  that  will  be  tested  by  the  new  creep/thermal  cycler  is  W-ThC>2 
fiber/FeCrAlY  matrix  composite,  which  is  durable  up  to  1 100°C.  The  above  MMC  specimens  have 
been  processed  through  powder  metallurgical  processing  route.  Due  to  the  limited  amount  of 
specimens,  the  specimen  geometry  of  Fig.  2.3  will  be  used. 
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2.2.  Thermal  Cycling  Resistance  of  SCS-6/Ti-15-3  Metal  Matrix  Composite 

2.2.1  Introduction 

In  order  to  achieve  high  performance  aerospace  structures,  materials  possessing  high 
strengths  with  low  densities  (high  specific  strengths),  combined  with  high  temperature  resistance, 
are  required.  This  is  one  reason  for  the  relatively  recent  development  of  metal  matrix  composites 
(MMC)  systems.  Ceramic  fibers  in  typical  MMC’s  possess  a  coefficient  of  thermal  expansion 
(CTE)  much  lower  than  the  metal  matrix.  This  mismatch  of  CTE’s  has  the  potential  to  cause 
cracking  whenever  this  material  is  thermally  cycled  due  to  the  strain  on  the  fiber-matrix  interface. 
This  has  been  observed  in  various  studies  which  have  shown  that  accumulated  damage  can  occur 
due  to  thermal  cycling  and,  in  certain  cases,  have  resulted  in  the  reduction  of  their  mechanical 
properties  [1-6].  In  these  particular  studies,  the  MMC’s  were  thermally  cycled  at  elevated 
temperatures  without  external  loading  which  caused  a  measurable  elongation  in  the  longitudinal 
fiber  direction  creating  porosity  and  cracking  at  the  fiber-matrix  reaction  zone.  Both  of  these 
damage  types  were  determined  to  be  the  direct  result  of  CTE  mismatch  between  the  fiber  and  the 
matrix. 

In  most  thermal  cycling  studies  the  specimens  were  exposed  to  the  maximum  temperatures 
for  short  time  periods  by  using  either  “sine  wave”  or  “saw  tooth”  temperature  profiles.  But  in 
actual  applications,  these  materials  would  likely  be  subjected  to  a  thermal  profile  resembling  a 
“square  wave”  with  a  dwell  time  at  the  maximum  cycling  temperature.  This  would  then  subject  the 
specimens  to  a  combination  of  sustained  high  temperature  effects  and  thermal  cycling  and 
presumably  give  rise  to  severer  damage. 

2.2.2  Experimental  procedure 

The  SCS-6  fiber/Ti- 15-3-3-3  (Ti- 15-3)  MMC  system  was  used  for  this  study.  This  alloy  is 
generally  limited  to  applications  below  550°C  -  600°C  because  of  alloy  embrittlement  due  to  oxygen 
diffusion  through  the  surface  [7].  The  SCS-6  fibers  were  manufactured  at  Textron  using  a 
chemical  vapor  deposition  (CVD)  process  starting  with  a  33  pm  diameter  carbon  filament,  upon 
which  SiC  was  deposited.  After  the  SiC  was  deposited,  a  complex  outer  coating  of  SiC  and 
Carbon  giving  the  fiber  and  overall  diameter  of  145  pm  [8].  This  outer  coating  is  applied  to  reduce 
the  surface  flaws,  allow  for  easier  handling,  and  to  limit  the  reaction  zone  between  the  fibers  and 
the  titanium  matrix. 

The  reaction  zone  appears  in  almost  all  SiC/Ti  MMC’s  as  an  interface  which  initially  occurs 
as  a  result  of  processing  or  consolidation  of  the  metal  and  fibers  into  MMC’s.  The  metal  matrix 


6 


consolidation  is  usually  done  at  high  temperatures  (1 100°C)  where  the  outer  coating  of  the  fibers 
reacts  with  the  titanium  in  the  matrix  creating  a  zone  which  can  start  as  low  as  727°C.  This  reaction 
zone  is  composed  of  small  particulates  which  consist  of  TiC  adjacent  to  the  fiber  and  Ti5Si3  next  to 
the  matrix.  In  addition,  some  small  “islands”  of  reaction  products  usually  exist  within  the  matrix 
adjacent  to  the  fibers.  This  reaction  zone  is  the  fiber-matrix  interface  that  transfers  the  loads  from 
the  matrix  into  the  fibers  and  can  greatly  influence  the  mechanical  behavior  of  the  MMC  [8]. 

The  tensile  specimens  were  fabricated  from  panels  manufactured  by  Textron  and  donated  by 
Boeing  for  this  study.  The  panels  consisted  of  eight  alternating  foils  of  Ti- 15-3  with  SCS-6  fibers 
in  unidirectional  orientation.  The  panels  were  determined  to  have  a  fiber  volume  fraction  of 
approximately  36%  -40%  based  on  previous  analysis  by  Boeing.  No  heat  treatment  was  given  to 
the  panels  prior  to  fabrication  or  testing  of  the  specimens.  Specimens  were  cut  from  the  panels 
using  electrical  discharge  machining  (EDM)  because  usual  machining  processes  could  cause 
excessive  fiber  pullout  in  the  curvatures  of  the  specimen  shape.  No  attempt  was  made  to  grind  off 
any  heat  affected  zone  that  may  have  developed  during  the  EDM  process.  The  tensile  specimens 
were  in  a  flatbar  type  configuration  with  the  fiber  oriented  in  the  longitudinal  direction.  They  were 
76.2  mm  long  by  12.7  mm  wide  with  a  gage  region  5  mm  wide  by  12.7  mm  long. 


Figure  2.4  Typical  thermal  cycling  temperature  profile 

The  specimens  were  thermally  cycled,  unloaded,  into  and  out  of  a  100  mm  diameter  vertical 
furnace  tube  heated  using  four  SiC  heating  rods.  Two  thermocouples  were  used  for  temperature 
control;  one  was  placed  adjacent  to  the  furnace  wall  for  power  supply  control,  and  the  other 
attached  to  the  cycling  specimens  to  determine  the  appropriate  settings  for  the  furnace  power 
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supply.  The  specimens  were  attached  to  a  rod  that  was  moved  into  and  out  of  the  furnace  by 
means  of  a  motor/pulley  system  using  a  mechanical  timer. 

Each  thermal  cycle  combined  a  heat  cycle,  approximately  5  minutes,  with  a  7  minute  non- 
forced  air  cooling  cycle.  The  overall  thermal  cycle  temperature  profile  resembled  a  “square  wave” 
having  an  approximately  3  minute  dwell  at  the  maximum  thermal  cycle  temperature  (Tmax ),  Fig. 
2.4.  The  specimens  were  exposed  to  Tmax  for  25,  50  and  75  hours  for;  500,  1000,  and  1500 
cycles,  respectively. 

Four  maximum  elevated  cycling  temperatures  were  used;  500°C,  600°C,  650°C,  and  700°C. 
Three  specimens  were  cycled  at  each  maximum  cycling  temperature  with  the  removal  of  one 
specimen  after  every  500  cycles.  In  addition,  dimensions  of  the  specimens  were  determined 
periodically  during  the  thermal  cycling  testing. 

2.2.3  Results 

The  results  of  the  tensile  testing  are  shown  in  Fig.  2.5  with  each  point  representing  the 
results  of  one  tensile  specimen  test.  The  failure  stresses  appear  to  decreased  consistently  with 
increasing  maximum  cycling  temperatures  by  about  250  MPa  for  each  different  setting  of  Tmax- 
Samples  of  specimens  were  examined  metallographically  revealing  a  greater  change  in  the  matrix 
microstructures  with  increasing  maximum  cycling  temperatures  when  compared  with  the  uncycled 
specimens.  The  uncycled  specimen  exhibited  a  typical  as-received  Ti- 15-3  microstructure  (Fig. 
2.6a)  which  contained  a  few  dark  blotchy  regions  distributed  throughout  the  matrix.  These  regions 
were  probably  the  solute  lean,  coherent  (T  phase  [9]. 

Starting  with  the  500°C  specimens,  the  a  phase  began  to  precipitate  in  a  very  fine  needle  form 
from  the  |T  phase.  (Fig.  2.6b).  These  precipitates  gradually  increased  in  dimensions  as  seen  in  the 
600°C/1500  cycle  microstructure  then  became  extremely  coarse  as  the  maximum  cycling  and 
temperatures  were  increased  to  700°C/1500  cycles.  The  largest  a  phase  precipitates  were  located 
adjacent  to  the  grain  boundaries  and  the  fibers.  However,  no  cracks  or  other  damage  was 
observed  in  the  matrix  of  any  of  the  specimens. 

In  addition,  an  oxide  layer  covering  the  specimens  was  observed  which  increased  in  thickness 
with  higher  maximum  cycling  temperatures  and  cycles.  But  even  at  the  most  extreme  case,  700°C 
at  1500  cycles,  the  oxidation  did  not  penetrate  beyond  the  first  fiber  layer. 

The  fiber-matrix  reaction  zones  do  not  appear  to  change  with  either  increasing  cycles  or  maximum 
temperatures.  In  addition,  no  debonding  or  damage  was  observed  at  the  reaction  zone.  The  fibers 
were  neither  affected  by  the  thermal  cycling,  other  than  cracks  observed  in  a  few  fibers  which 
corresponded  directly  to  the  mechanical  test  machine  grip  impressions  made  on  the  surface  of  the 
specimens. 
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Figure  2.5  Residual  tensile  strength  of  thermal  cycled  specimens 

Evaluation  of  the  fracture  surface  morphology  of  the  tensile  specimens  revealed  no  significant 
differences  regardless  of  the  cycling  conditions.  The  matrix  adjacent  to  the  fibers  formed  sharp 
ridges  consistent  with  plastic  flow.  Dimpling,  also  indicative  of  plastic  overload,  was  present  on 
all  of  the  matrix  surfaces  whether  it  was  on  the  ridges  adjacent  to  the  fibers  or  between  the  fibers. 
The  reaction  zones  on  the  other  hand,  appeared  to  have  fractured  in  a  britde  mode.  These  materials 
did  not  adhere  consistently  to  either  the  outer  coatings,  on  the  fibers,  or  the  matrix.  Furthermore,  a 
fracture  region  between  the  carbon  fiber  core  and  the  SiC  region  was  clearly  observed.  This 
interface  was  not  easily  characterized  in  the  metallographic  examination.  It  appears  to  be  similar  to 
the  fiber-matrix  reaction  zone  in  that  it  adhered  to  both  fracture  surfaces  inconsistently. 

The  metallographic  mounts  of  the  uncycled  and  the  700°C/1500  cycled  specimen  were  placed 
in  a  scanning  electron  microscope  for  further  evaluation  of  the  reaction  zone.  Like  the  optical 
examination,  no  significant  microstructural  changes  could  be  detected  at  a  7500X  magnification. 
The  only  notable  change  was  the  configuration  of  the  outer  layer  of  the  fiber.  The  uncycled  coating 
on  the  fiber  appeared  to  be  in  three  layers  probably  composed  of  equal  compositions  based  on  the 
signal  strength  detected  by  the  SEM  which  corresponding  to  their  atomic  numbers.  The 
700°C/1500  specimen  reaction  zone  appeared  to  have  only  two  such  layers  but  it  had  the  same  total 
thickness  (approximately  4.5  pm)  as  the  un-cycled  reaction  zone.  The  reaction  zones  also  appeared 
to  have  same  average  thickness  and  possessing  similar  jagged  appearance. 

Most  of  the  specimens  developed  overall  dimensional  changes  during  the  cycling.  All  of  the 
recorded  specimen  dimensions  generally  increased  with  increasing  maximum  cycling  temperatures 
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Figure  2.6  a)  Microstructure  of  the  uncycled  specimen  and  b)  microstructure  of  the  700°  C/1500 
cycled  specimen 
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and  cycles.  The  maximum  dimensional  change  was  0.215  inch  on  the  700°C/1500  specimen  at  the 
gage  thickness.  However,  the  length  or  longitudinal  direction  measurements  did  not  increase 
consistently  with  the  other  measurements.  A  comparison  of  the  micrometer  measurements  with  the 
oxide  measurements  from  the  micrographs,  indicated  that  the  specimen  growth  was  due  to  the 
oxidation  products. 

2.2.4  Discussion 

Since  no  general  damage  or  porosity  was  observed,  this  would  suggest  that  a  microstructural 
or  chemical  change  occurred  within  the  composite  material  which  would  cause  a  reduction  in  the 
tensile  properties.  There  are  three  primary  material  components  within  the  MMC;  the  fibers, 
reaction  zone,  and  the  Ti- 15-3  matrix,  which  can  affect  the  properties  of  the  composite  system. 

From  observations,  the  fibers  did  not  appear  to  have  been  affected  either  by  the  increase  in 
the  number  of  cycles  or  the  maximum  cycling  temperature.  No  cracks  or  microstructural  changes 
were  observed  in  the  cross  sections,  either  in  the  core  or  the  bulk  of  the  fibers,  which  would 
suggest  that  the  thermal  cycling  did  not  effect  the  fiber  material  to  any  noticeable  degree.  This  is 
consistent  with  other  studies  that  have  shown  that  these  fiber  types  by  themselves  exhibited  little 
change  in  their  mechanical  properties  up  to  approximately  1100°C  [10].  Since  these  specimens 
were  cycled  at  considerably  lower  temperatures  it  would  appear  that  the  fibers  themselves  were  not 
affected. 

The  fibers’  strength  within  the  matrix  can  be  greatly  affected  by  the  properties  of  the  fiber- 
matrix  reaction  zone.  Yet  observations  of  the  reaction  zones  did  not  reveal  any  microstructural 
changes  and  these  regions  were  probably  not  effected  by  the  thermal  cycling.  In  other  work,  [11] 
it  was  reported  that  the  reaction  zone  of  the  SCS-6/Ti-15-3  system  can  actually  increase  in  strength 
under  specific  thermal  aging  conditions..  In  a  SCS-6/Ti-15-3  composite  system  the  interfacial 
shear  strength  increased  from  approximately  118  MPa  to  131  MPa,  at  800°C  from  12  hours  to  50 
hours.  This  would  increase  the  overall  strength  of  the  fibers  and  therefore  the  composite  as  it 
would  allow  for  better  transfer  of  the  loads  to  the  fibers.  But  since  the  composite  material  in  this 
study  exhibited  a  decrease  in  strength,  any  possible  changes  to  the  interfacial  shear  strength  of  the 
reaction  zone  did  not  seem  to  have  any  affect. 

Since  the  fibers  and  the  reaction  zone  did  not  appear  to  change  the  mechanical  properties  of 
the  MMC,  the  only  component  of  the  composite  that  most  likely  contributed  to  the  reduced  strength 
in  the  composite  was  the  matrix.  Heat  treatment  studies  of  SCS-6/Ti-15-3  have  shown  that  the 
matrix  will  decrease  in  strength  over  time  due  to  phase  changes  [9].  The  study  compared  un¬ 
reinforced  Ti-15-3  alloy  properties  to  SCS-6/Ti-15-3  composite  properties  by  aging  the  materials  at 
temperatures  ranging  from  300°C  to  700°C  for  24  hours.  It  revealed  that,  based  on  hardness 
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measurements,  the  strength  of  the  MMC  would  increase  to  a  plateau  between  400°C  and  500°C  and 
then  decrease  with  increasing  aging  temperatures.  In  a  related  study  [12],  this  MMC  was  aged  for 
the  times  and  temperatures  shown  in  Table  2.1  with  the  resulting  ultimate  tensile  strengths: 


Table  2.1 


T  emperature/Time 

UTS  (MPa) 

As  Fabricated 

1572 ± 125 

500°C/50  hr 

1728  ± 130 

500°C/100  hr 

1607  ± 104 

800°C/12  hr 

1331  ±41 

800°C/50  hr 

1269  ±  212 

The  thermal  treatments  shown  above  most  likely  causes  the  solute  rich  alloy  to  separate  into 
solute-rich  p  and  a  solute-lean  p’  phases  prior  to  the  precipitation  of  the  a  phase  from  the  P‘ 
crystals.  The  a  phase  provides  strengthening  to  the  alloy  but  requires  long  aging  time  due  to  the 
slow  formation  kinetics.  Ti-15-3  is  also  susceptible  to  grain  boundary  formation  of  the  a  phase. 
This  is  due  to  the  nucleation  kinetics  of  the  a  in  the  grain  boundaries  which  is  not  as  sensitive  to  p 
stabilizers  as  the  kinetics  for  the  formation  of  intergranular  a  [13].  Additional  aging  coarsens  the 
a  phase  which  in  turn  reduces  its  strain  to  failure  and  eventually  reducing  the  overall  strength  of  the 
alloy. 

Since  the  residual  strengths  of  the  specimens  doesn’t  change  appreciably  for  a  given 
maximum  thermal  cycling  temperature  level  after  the  first  500  cycles,  it  appears  that  the  reduction 
in  strength  of  this  MMC  occurred  during  the  first  500  cycles. 

In  comparison  to  a  similar  thermal  cycling  study  [14],  the  tensile  specimens  exhibited  no 
appreciable  change  in  the  residual  strength  even  after  undergoing  more  cycles  than  specimens  had 
undergone  for  this  paper.  In  that  study,  two  sets  of  tensile  specimens  were  cycled  from  149°C  to 
427°C  and  from  149°C  to  649°C,  for  up  to  16,750  cycles  yet  the  tensile  strengths  varied  between 
1750  MPa  and  1900  MPa  with  no  correlation  to  the  number  of  cycles  or  the  maximum  cycling 
temperature  level.  The  composite  lay-up  of  eight  plies  was  similar  to  the  one  used  for  this  study 
however  their  tensile  specimens  were  rectangle  in  shape,  12.7  mm  by  177  mm,  with  no  specific 
gage  length  or  necked  region.  In  addition,  these  specimens  started  cycling  in  the  “as-received” 
condition  and  in  an  unloaded  condition.  However,  the  specimens  were  cycled  using  a  48  second 
cycle  “saw  tooth”  temperature  profile  with  no  dwell  time  at  either  the  high  or  low  temperatures.  It 
appears  that  their  specimens  were  probably  not  overaged  due  to  the  limited  time  at  the  elevated 
temperature  and  therefore  did  not  suffer  a  reduction  in  tensile  strength.  The  different  cycling 
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thermal  profiles  used  in  these  studies  may  be  responsible  for  the  discrepancy  of  the  results. 

In  another  study  [15]  using  a  similar  thermal  cycling  temperature  profile,  33.3  hours  at  Tmax, 
the  reduction  in  the  residual  strength  of  the  specimens  (SCS-6/Ti-6Al-4V)  was  attributed  to 
changes  of  the  matrix,  primarily  oxygen  embrittlement,  which  subsequently  caused  a  degradation 
of  the  fibers. 

A  calculation  was  performed  on  the  tensile  data  using  the  rule  of  mixtures,  Eq.  (2.1),  to 
determine  whether  the  reduction  in  the  matrix  strength  could  be  directly  correlated  to  the  overall 
composite  strength: 


gc  =Vfof  +(1-Vf)cm 


(2.1) 


The  matrix  and  fiber  sttength  values  used  for  the  calculations  were  obtained  from  previous  studies 
[9,16].  As  seen  in  Fig.  2.7,  the  results  do  not  correspond  generally  for  the  data  beyond  the  600°C 
strength  levels.  There  are  two  possible  reasons  for  this  discrepancy.  One  is  that  the  rule  of 
mixtures  is  a  simple  model  that  may  not  take  into  account  the  interactions  of  other  variables  that 
could  affect  this  specific  composite,  such  as  the  interfacial  shear  strength.  Another  reason  would 
most  likely  be  the  nonuniform  matrix  changes  such  as  those  observed  in  Fig.  2.6b.  The  precipitate 
lean  region  surrounding  the  fibers  would  exhibit  a  lower  mechanical  strength  due  to  fewer  a 
precipitates  which  would  reduce  the  overall  matrix  strength.  The  effective  matrix  flow  stress  that  is 
influential  in  governing  the  fracture  of  a  MMC  is  that  of  the  matrix  adjacent  to  the  fiber.  Thus  if 
<7m  in  equation  (1)  is  replaced  by  the  matrix  flow  stress  of  the  soft  region  surrounding  the  fiber, 
then  the  rule  of  mixtures  suggests  that  the  strength  of  the  composite  is  reduced  due  to  the  existence 
of  the  soft  region. 


Figure  2.7  Plot  of  the  rule  of  mixtures  predictions  compared  to  the  experimental  data 
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2.2.5  Conclusions 


1.  In  regards  to  the  fibers,  matrix  and  the  reaction  zone,  only  the  microstructure  of  the  Ti-15-3 
matrix  exhibited  a  change  which  can  be  related  primarily  to  the  exposure  of  the  specimen  to  a 
duration  of  accumulated  thermal  cycles  at  the  maximum  cycling  temperature. 

2.  The  microstructural  changes  of  the  Ti-15-3  matrix  is  the  most  likely  explanation  for  the 
reduction  in  residual  strength  of  the  specimens  and  due  to  prolonged  exposure  to  the  higher 
temperatures  used  in  this  study. 
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2.3  Dimensional  Change  of  SiCp/AI  Composite 

The  present  MMC  specimens  made  of  SiC  particulate/aluminum  (A359)  composite  were 
processed  by  using  the  DURALCAN  molten  metal  mixing  technology  by  Alcan  International  Ltd., 
followed  by  extrusion  with  the  ratio  of  9:1.  The  volume  fraction  of  the  particulates  (f)  is  0.2.  The 
range  of  the  longer  diameter  of  angular  shaped  SiC  particulates  is  from  10  jam  to  15  pm.  As- 
extruded  SiC  particulate/Al  (SiCp/Al)  composite  were  subsequently  all  solution-treated  at 
temperature  530  °C.  The  same  processing  route  was  also  used  to  fabricate  control  unreinforced 
matrix  material[l].  As-extruded  bar  was  further  machined  to  a  buttonhead  specimen  with 
shouldered  ends  whose  loading  axis  is  parallel  to  the  extrusion  direction.  The  dimensions  of  the 
gage  section  of  a  cylindrically  shaped  test  specimen  are  26.99  mm  in  length  and  3.81  mm  in 
diameter. 

Thermal  cycling  creep  tests  were  performed  on  SiCp/Al  composite  specimens  for  two 
temperature  profiles  of  Tmax=300°C  and  350°C  by  using  a  thermal  cycler  with  constant  stress 
capability,  consisting  of  a  movable  furnace  of  a  constant  temperature,  temperature  controller,  load 
train  and  a  cooling  system[2].  Since  the  creep  thermal  cycling  tests  of  MMCs  are  expected  to  result 
in  a  large  deformation,  a  constant  stress  cam  was  employed  to  compensates  the  instantaneous 
change  of  the  cross-sectional  area  of  specimen  so  that  the  specimen  can  be  loaded  always  under  a 
constant  stress  condition  during  tests.  A  specimen  is  thermally  cycled  by  the  moving  furnace  and 
the  forced  cooling  system  where  the  heating  and  cooling  durations  are  kept  as  6.3  and  3.2  min., 
respectively.  A  typical  temperature  profile  was  measured  by  the  K-type  thermocouple  embedded  in 
the  specimen,  shown  in  Fig.  2.8  where  (a)  and  (b)  denote  the  cases  of  Tmax=350  °C  and 
Tmin=115  °C,  and  Tmax=306  °C  and  Tmin=100  °C,  respectively.  The  detail  explanation  of  the 
creep  thermal  cycler  is  given  elsewhere[2]. 

The  dimensional  change  of  SiCp/Al  composite  has  been  known  to  be  proportional  to  the 
number  of  thermal  cycle  (N)  for  a  constant  applied  stress  up  to  several  hundred  percent  extension 
without  any  fracture  and  necking  so  that  tests  were  performed  up  to  the  strain  of  about  0.5[3,4]. 
The  dimensional  change  of  a  specimen  was  obtained  by  measuring  the  gauge  length  change  after 
cycles  with  a  digital  vernier  caliper  of  resolution  10  pm  and  experimental  results  are  shown  in  Fig. 

2.9.  With  the  initial  specimen  gauge  length  of  26.99  mm  the  smallest  resolvable  strain  is 
approximately  4  x  10"4.  Although  the  strain  measurement  is  intermittent  and  its  resolution  is  not  so 
high,  this  strain  measurement  system  is  considered  to  be  adequate  because  the  purpose  is  to 
measure  the  dimensional  change  at  different  thermal  cycles  (N),  for  example,  N  =  100,  500, 1000. 

The  microstructural  studies  were  also  conducted  using  scanning  electron  microscope  (JEOL- 
T330A)  to  compare  the  microstructural  change  in  terms  of  particulate  aspect  ratio  and  particulate 
orientation  before  and  after  tests.  The  typical  micrographs  are  shown  in  Fig.  2.10,  where  (a),  (b) 
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and  (c)  denote  as-received  state,  as-tested  states  under  the  applied  stress  of  o°=24  MPa  and 
Tmax=300  °C  and  Tmax=350  °C,  respectively.  The  two  tested  specimens  of  Figs.  2.10  (b)  and  (c) 
were  deformed  to  the  elongations  of  56%  and  77  %  without  fracture,  respectively.  Based  on 
micrographs  having  the  total  number  of  604  particulates,  a  two-dimensional  analysis  was 
performed  to  determine  the  average  aspect  ratio  and  orientation  of  particulates.  The  distributions  of 
particulate  aspect  ratio  are  shown  in  Fig.  2.11,  where  average  aspect  ratios  are  2.15,  2.28  and  2.34 
for  Figs  2.10  (a),  (b)  and  (c),  respectively.  The  results  of  the  microstructural  study  are  used  as 
input  data  for  the  analytical  study. 


9  10  20  30 

Applied  stress  (MPa) 


Figure  2.9  The  experimental  results  of  SiCp/Al  composite  subjected  to  creep/thermal  cycle  loading: 
strain  per  thermal  cycle  vs.  applied  stress 
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Figure  2.10  Typical  microstructure  of  SiCp/Al  composite  observed  by  SEM:  (a)  as-received  state, 
(b)  as-tested  under  the  applied  stress  of  a°=24  MPa  and  Tmax=300  °C  and  (c)  Tmax=350  °C, 
respectively,  where  the  horizontal  direction  (<->)  coincides  with  the  directions  of  the  extrusion  and 
applied  stress. 


Figure  2.11  Histograms  showing  the  distribution  of  fiber  aspect  ratios  from  the  micro  structural 
studies  on:  (a)  the  as-received,  (b)  as-tested  under  o°=24  MPa  and  Tmax=300  °C  and  (c) 


Tmax=350  °C,  respectively.  Average  aspect  ratio  is  determined  to  be  2.25. 


19 


References 


1.  D.  J.  Lloyd:  In  High  Performance  Composite:  Commonalty  of  Phenomena,  edited  by  K.  K. 
Chawla,  P.  K.  Liaw  and  S.  G.  Fishman,  The  Minerals,  Metals  and  Materials  Society,  1994,  pp. 
465-473. 

2.  M.  Taya  et  al,  "Creep  and  Thermal  Cycling  Fixture  Design  for  Metal  Matrix  Composites,"  L 
Exp.  Mech..  35  (1995),  66-70. 

3.  M.  Y.  Wu  and  O.  D.  Sherby,  "Superplasticity  in  a  Silicon  Carbide  Whisker  Reinforced 
Aluminum  Alloy,"  Scripta  Metall.,  18  (1984),  773-776. 

4.  S.  M.  Pickard  and  B.  Derby,  "The  Deformation  of  Particle  Reinforced  Metal  Matrix 
Composites  During  Temperature  Cycling,"  Acta  Metall.,  38  (1990),  2537-2552. 


20 


2.4  Fracture  Toughness  of  PSZ-SS/PSZ  Composite-PSZ  Laminate* 


Ceramics  has  not  only  high  stiffness,  high  compressive  strength  and  high  hardness,  but  also 
good  fatigue  resistance,  good  creep  and  corrosion  resistance.  These  material  properties  make 
ceramics  very  useful  in  engines,  cutting  tools  and  some  structure  designs  in  aerospace  engineering. 
Despite  the  above  advantages,  the  toughness  or  crack  growth  resistance  of  ceramics  is  low, 
typically  2~5MPaVm .  This  disadvantage  really  weakens  the  load  carrying  capacity  of  a  ceramic 
structure  and  impedes  the  wide  application  of  this  material.  One  of  the  methods  to  improve  the 
fracture  toughness  is  to  reinforce  the  ceramics  with  fillers  (fibers  and  particles)  and  form  a  ceramic 
matrix  composite  (CMC).  The  microstructure  in  a  CMC  is  normally  homogeneous,  i.e.  the  volume 
fraction  of  fillers  is  uniform  spacially.  It  would  be  interesting  to  examine  the  fracture  toughness  as 
a  function  of  crack  extension  (so  called  R-curve)  of  CMC  if  the  microstructure  of  the  CMC  is  not 
uniform.  An  example  of  such  a  CMC  with  non-homogeneous  micro  structure  is  shown  in  Fig.  2.11 
where  pure  ceramic  (partially  stabilized  zircornia,  PSZ  in  short)  sandwiches  stainless  steel/  PSZ 
composite  (ceramic  matrix  composite,  SS/PSZ  composite  in  short).  It  would  be  valuable  to  assess 
the  effect  of  the  SS/PSZ  composite  layer  on  the  R-curve  behavior  of  the  laminate  which  is  expected 
to  exhibit  a  transition  behavior  at  the  PSZ  ceramic-SS/  PSZ  composite  interface  as  schematically 
indicated  in  Fig.  2.11.  The  purpose  of  this  study  is  to  assess  the  R-curve  behavior  of  the  PSZ- 

SS/PSZ  composite-PSZ/laminate  composite  experimentally. 

Pre-crack 


Figure  2.1 1  PSZ-SS/PSZ  composite-PSZ  laminate  and  expected  R-curve* 


*  The  results  of  this  work  were  published  in  J.  Amer.  Ceram.  Soc.,  Vol.  78  (6),  1995,  pp. 1633- 
1639,  "Experimental  Study  of  the  Fracture  Toughness  of  a  Ceramic/Ceramic-Matrix  Composite 
Sandwich  Structure,"  by  Z.  Liu,  M.  Taya,  M.  L.  Dunn,  and  R.  Watanabe. 
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2.4.1  Experiment 


Loading  Device  and  WL-DCB  Specimen 

The  loading  device  shown  in  Fig.  2. 12  can  load  or  unload  a  specimen  respectively  by  clock¬ 
wise  or  counter  clock-wise  turn  of  the  arm  which  is  sensitively  controlled  by  turning  of  the  drive 
screw  handle.  The  load  is  determined  by  the  strain  in  the  load  cell  which  is  connected  with  the 
strain  gage  indicator — model  3270,  whose  output  numbers  in  the  display  are  linearly  related  to  the 
virtual  load  (P)  on  the  specimen. 

Unlike  some  other  loading  device  in  which  the  load  is  controlled  by  weight,  this  load  is 
controlled  by  the  displacement  at  the  load  point  between  loading  pin  and  surface  of  the  specimen. 
When  the  load  increases,  the  displacement  is  increased,  in  the  meantime,  the' increased 
displacement  will  reduce  the  load  a  little.  This  device  is  called  displacement  controlled  loading 
device,  or  COD-feedback  loading  device,  since  this  displacement  is  related  to  crack  opening 
displacement  (COD). 

In  order  to  control  the  speed  of  crack  propagation,  wedge-loaded  double  cantilever  beam 
(WL-DCB)  specimen  with  Chevron  notch  is  selected  (Fig.  2.13(a)).  The  cutting  process  is 
accomplished  by  using  the  diamond  blade  in  the  machine  LEEMATTC  2000  surface  grinders/slices. 

Some  improvements  are  made  in  the  geometry  and  dimension  of  the  specimen.  The  Chevron 
notch  is  replaced  by  the  triangular  notch  (Fig.  2.13(b)),  the  V  notch  angle  is  reduced  from  60°  to 
45°  and  the  width  of  the  specimen  is  cut  significantly  into  a  narrow  one.  By  doing  so,  one  can 
obtain  a  pre-crack  more  easily. 

Moire  Interferential  Optical  System 

In  order  to  measure  the  fracture  toughness,  one  can  use  data  of  COD  of  a  specimen  and  the 
results  of  finite  element  analysis  (hybrid-experimental-numerical  method).  Moire  interferential  is  a 
reliable  displacement  measurement  method  which  can  measure  COD  of  the  specimen  used  in  this 
experiment.  The  details  of  the  principle  of  the  Moire  interferometry  can  be  found  elsewhere  (Post, 
1987). 

Moire  interference  can  be  produced  by  several  different  optical  systems.  The  main  part  of  the 
work  is  to  obtain  two  coherent  parallel  lights  from  the  right  angle.  In  this  experiment,  a  simple 
optical  set-up  as  shown  in  Fig.  2.14  is  used,  where  the  second  beam  of  coherent  light  is  obtained 
with  assistance  of  mirror  2  instead  of  two  laser  radiators.  The  function  of  the  special  filter  and  lens 
is  to  get  a  disk  of  parallel  light,  which  should  cover  Mirror  2  and  the  specimen.  The  red  laser  beam 
comes  from  the  Model  159  lOmw  Helium-Neon  laser  generator.  Due  to  low  power,  there  is  no  fan 
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in  the  radiator.  This  is  very  important,  because  if  the  vibration  is  brought  from  the  fan,  it  will  affect 
the  fringes.  In  this  experiment,  the  grating  frequency  is  1200  lines  per  millimeter,  and  the 
wavelength  of  the  He-Ne  red  laser  light  is  6328nm  which  is  equal  to  6328xlO-10m,  hence  the  angle 
of  the  incident  light  is  given  by 


.  ,fl  .  ,  1200  x  6328  x  10'7  „„  ^ 
a  =  sin- 1  =  sin' 1 - 55 - =  22.3 


Since  the  fringe  pattern  directly  represents  the  displacement  field,  it  is  used  to  measure  CODs 
which  are  in  turn  used  to  compute  elastic  energy.  Following  the  principle  of  the  moire  fringe  in 
deformed  field,  the  x-direction  displacement  ( u )  can  be  found  in  the  fringe  pattern,  and  given  by 

u=NP  (2.3) 


Mirror  1 


(b) 


Figure  2.14  Moire  interferometry  system  for  displacement  measurement. 
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Figure  2.15  An  example  of  Moire  pattern. 

Where  N  is  fringe  order  and  P  is  pitch  of  the  grating.  There  exists  a  difference  between  geometry 
Moire  and  Moire  interferometry.  In  the  Moire  interferometry,  the  pitch  of  the  grating  must  be  used 
as  half  of  the  real  pitch.  For  example,  if  the  pitch  of  the  grating  (P)  is  1/1200  millimeter,  when  the 
displacement  is  calculated  by  Eq.  (2.2),  the  effective  pitch  (Peff)  should  be  1/2400  millimeter.  Fig. 
2.15  is  one  of  the  moire  fringe  patterns  observed  in  this  experiment. 

Basically,  the  zero-order  fringe,  which  has  no  x-direction  displacement  on  the  fringe,  is  a 
bright  straight  line  at  the  specimen's  symmetric  center  connected  with  the  crack  tip.  Next  to  the 
zero-order  fringe  is  the  dark  half-order  fringe,  which  represents  displacement  of  the  half  pitch  of 
the  effective  grating. 

u  =  NPeff  =  N ^ f  - \xjm  =  °-2083  x  10'3  (mm)  (2,4) 

Similarly,  the  fringe  order  at  load  point  can  be  counted,  for  example  for  N  =  56.5,  the 
displacement  is 

u  =  NPeff  =  56.5  x  ^qq  =  23.542  x  10*3  (mm)  (2.5) 

Since  the  load  and  geometry  of  the  specimen  are  symmetrical,  only  a  half  of  the  specimen  is 
used  to  compute  the  R-curve.  However  in  reality,  due  to  the  error  of  the  grating  coating,  the  fringe 
pattern  is  often  not  symmetrical.  Thus  the  average  number  is  used  to  represent  x-displacement. 
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Experimental  Procedure 


First,  one  must  adjust  the  angle  of  the  mirror  2  in  Fig.  2.14  to  obtain  horizontal  parallel 
fringe  lines  when  the  load  is  zero.  The  fewer  these  lines,  the  more  accurate  the  optical  system  is 
adjusted.  With  the  increase  of  the  load,  both  the  numbers  of  the  fringes  and  the  shape  of  the 
fringes  are  changed.  In  order  to  keep  the  fringe  pattern  symmetrical  in  this  procedure,  the  angle  of 
mirror  2  must  be  adjusted  every  time  when  an  increase  of  load  is  applied.  But  only  one  direction  of 
the  angle  needs  to  be  adjusted,  either  up  and  down  or  left  and  right;  otherwise,  CODs  of  the 
specimen  are  not  accurately  shown  by  the  fringe  pattern.  In  fact,  the  absolute  symmetric  fringe 
pattern  is  very  hard  to  maintain  during  the  experiment,  so  the  CODs  measured  at  both  edges  are 
averaged. 

Without  a  successive  stable  crack  growth,  the  R-curve  can  not  be  obtained.  Thus  the  load 
must  be  increased  very  slowly. 

Computation  of  Energy  Release  Rate 

Energy  release  rate  is  defined  by 


r  d(F-U) 

a  —&r 


(2.6) 


where  F  is  work  done  by  external  force,  U  is  elastic  energy  contained  in  the  specimen  and  a  is  the 
crack  length.  In  the  actual  computation,  the  derivative  in  the  Eq.  (2.7)  is  approximated  by  the 
incremental  form: 


Gn  = 


AF-AU 

Aa 


(Fn+1  —  Fn)  —  (Un+i  —  Un) 
(an+l  ~  an) 


(2.7) 


This  requires  the  experimental  procedure  of  continually  recording  the  load,  crack  length  and 
displacement  from  which  the  R-curve  is  computed. 

Since  PSZ  and  PSZ/SS  composite  are  both  brittle  material,  only  elastic  energy  is  considered 
in  this  case,  and  the  analysis  and  computation  can  be  based  on  linear  elastic  fracture  mechanics 
(LEFM).  In  this  fracture  test ,  load  in  y  direction  is  read  by  the  strain  gage  indicator  while  the 
displacement  in  x  direction  at  loading  point  by  the  Moire  fringe  pattern. 

External  work  and  elastic  energy  can  be  computed  by  a  finite  element  method  (FEM).  For  a 
crack  extension  Aa  in  the  specimen,  a  new  mesh  is  generated  and  new  external  work  and  elastic 
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energy  are  computed.  Then  by  the  use  of  Eq.  (2.4),  the  energy  release  rate  is  computed.  The 
contribution  to  the  energy  change  by  Y  component  of  the  force  (Fy)  is  so  small,  thus  it  can  be 
neglected  compared  with  the  contribution  by  the  displacement  in  x-direction.  Then  the  Eq.  (2.4) 
can  be  simplified  as 


Gn 


Un~Un+l 
^n+1  — 


-  Rn 


(2.8) 


In  this  equation,  an+i  is  always  larger  than  an  and  Un+i  is  always  less  than  Un,  hence  the  Rn  is 
always  positive. 

Due  to  the  symmetry  of  the  specimen  geometry  and  the  load,  only  a  half  of  the  specimen 
geometry  needs  to  be  considered.  Along  the  crack  growth  line,  the  denser  mesh  is  needed  and  the 
node  position  around  the  crack  tip  should  be  changed  to  the  quarter  node.  By  doing  so,  a  more 
accurate  computation  result  can  be  achieved. 


2.4.2  Result  and  Discussion 

R-curve  of  PSZ-SS/PSZ  Composite-PSZ  Laminate 


The  experiment  is  performed  on  three  specimens:  one  specimen  made  of  PSZ  only  and  two 
PSZ-SS/PSZ  composite-PSZ  laminate  specimens  with  two  different  volume  fraction  of  SS  in  the 
SS/PSZ  composite  area.  The  dimensions  of  the  specimens  are  identical  and  they  are  shown  in  Fig. 
2.16.  The  mechanical  properties  of  SS,  PSZ  and  SS/PSZ  composite  are  shown  in  Table  2.2, 
where  those  of  the  composite  are  estimated  by  using  the  Eshelby's  model  [6]. 

In  order  to  examine  the  effects  of  SS  phase  on  the  R-curve  behavior  in  PSZ-SS/PSZ 
composite-PSZ  laminate,  two  specimens  are  used:  one  contains  10%  SS,  another  one  20%  SS  in 
SS/PSZ  composite  layer.  As  a  reference,  a  pure  PSZ  specimen  is  also  used.  The  results  of  R-curve 
(energy  release  rate)  of  these  three  specimen  as  a  function  of  Aa  are  shown  in  Fig.  2.17,  where 
SS/PSZ  composite  layer  (CMC)  region  is  indicated  by  horizontal  bar. 

It  is  often  useful  to  convert  R-curve  to  KR-Aa  relation,  where  Kr  is  the  critical  stress 
intensity  factor  or  fracture  toughness  at  a  given  crack  length  Aa.  In  plane  stress  condition 

Kr  =  VER  (2.9) 
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where  E  is  modules  of  elasticity  and  R  is  crack  growth  resistance.  Note  in  this  formula,  different 
values  of  E  should  be  used  in  the  different  region.  Fig.  2.18  is  the  KR-Aa  relation  obtained  from 
Fig.  2.17  and  Eq.  (2.9). 


44.5° 


Figure  2.16  Dimensions  of  PSZ-SS/PSZ-PSZ  laminate  specimen. 


Table  2.2  Mechanical  Properties 


E  (GPa) 

Poisson's  ratio  v 

CTE  (10'6/C°) 

ss 

196 

0.26 

14.0 

PSZ 

168.8 

0.3 

11.1 

10%SS+90%PSZ 

171.4 

0.296 

11.3  j 

20%SS+80%PSZ 

174.0' 

0.292 

11.6 

It  is  clear  from  Fig.  2.17  and  2.18  that  the  energy  release  rate  (R)  and  fracture  toughness 
(Kr)  which  are  small  during  early  stage  of  crack  propagation,  jump  to  a  relatively  large  value  once 
the  crack  advances  the  CMC  region.  The  jump  in  terms  of  Kr  is  almost  twice  as  large  as  that  of 
unreinforced  PSZ  specimen  .  Fig.  2.17  and  2.18  also  show  the  gradual  increase  in  R-Aa  or  KR-Aa 
curve  in  early  stage  of  crack  propagation  characteristics  with  brittle  materials. 
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Crack  Growth  Resistance  GR  (Pam) 
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Figure  2.17  R(Gr)  -  Aa  relation. 
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Figure  2.18  Kr  -  Aa  relation. 
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3.  Analytical  Study  I:  Dislocation  Punching  from  Ceramic/Metal  Interfaces* 

3.1  Introduction 

When  two  different  materials  are  bonded,  misfit  strain  exists  at  the  interface  at  and  around 
which  high  stress  field  is  induced.  This  misfit  strain  at  the  interface  is  due  to  the  mismatch  in  the 
stiffness  of  the  materials  and  in  coefficient  of  thermal  expansion  (CTE)  under  temperature  change. 
The  order  of  the  misfit  strain  due  to  CTE  mismatch  can  be  quite  large  for  larger  temperature 
changes.  Microscopically  this  misfit  strain  at  the  interface  is  accommodated  by  dislocations  near  the 
interface.  When  the  magnitude  of  the  stress  field  becomes  too  large,  relaxation  of  the  stress  of  large 
magnitude  is  more  likely  to  take  place.  Punching  of  dislocations  from  the  interface'  is  such  a 
relaxation  mechanism  and  has  been  observed  in  dispersion-hardened  alloys  [1,2]  metal  matrix, 
composites  [3-5],  and  metal  coating/ceramic  substrate  system  [6]. 

The  dislocation  punching  model  was  first  used  by  Ashby  and  Johnson  [2]  to  explain  the 
strengthening  mechanism  of  dispersion-hardened  alloys  subjected  to  shear  (plastic)  strain  Y  where 
the  punching  of  prismatic  dislocation  loops  along  the  secondary  slip  planes  is  assumed.  The  work¬ 
hardening  predicted  by  the  Ashby’s  model  is 

x  -  Xy  =  0.24|iVbfW  (3.1) 

where  x  and  xy  are  the  flow  stress  and  initial  yield  stress  in  shear,  q,  b  and  Y  are  the  shear 
modulus,  the  Burgers’  vector  and  plastic  shear  strain  of  the  matrix,  respectively,  and  f  and  d  are 
the  volume  fraction  and  diameter  of  particles.  This  model  predicts  well  the  parabolic  dependence  of 
x  on  Y  observed  in  the  experiment.  The  misfit  strain  in  the  Ashby’s  model  before  punching  is 
represented  by  arrays  of  edge  dislocations  piled  up  at  the  particle/matrix  interface.  Tanaka  and  Mori 
[7]  solved  the  problem  of  misfit  strain  at  the  interface  by  using  the  Eshelby’s  model  [8]  where  the 
misfit  (plastic)  strain  at  the  particle/matrix  interface  is  smeared  out  to  become  “transformation 
strain”  [8]  or  “eigenstrain”  [9].  The  stress-strain  relation  of  a  composite  predicted  by  the  Tanaka- 
Mori  model  is  given  by 


Go  —  CTy  hep 


(3.2) 


*The  results  of  this  work  were  published  in  J.  Eng.  Mater.  Tech.,  Vol.  116,  July  1994,  pp.  408- 
413,  "Dislocation  Punching  from  Ceramic/Metal  Interface,"  by  M.  Taya  and  T.  Mori. 
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where  G0  and  Gy  are  the  flow  stress  and  initial  yield  stress,  h  and  ep  are  the  work-hardening  rate 
and  plastic  strain  along  the  loading  direction.  The  linear  relationship  between  G0  -  <3y  (work¬ 
hardening)  and  plastic  strain  ep  in  the  above  equation  agrees  with  the  experiment  for  small  range  of 
plastic  straining.  For  larger  plastic  straining,  however,  stress-strain  curve  deviates  from  the  linear 
work-hardening  predicted  by  the  Tanaka  Mori  model.  Tanaka  et  al.  [10]  proposed  a  model  to 
account  for  non-linear  work-hardening  at  larger  strains  by  using  a  dislocation  punching  model. 

Modeling  of  dislocation  punching  in  a  metal  matrix  composite  subjected  to  temperature 
change  AT  was  first  studied  by  Arsenault  and  Shi  [11]  who  used  an  equiaxed  particulate  (its  length 
t)  as  a  reinforcement  and  obtained  the  average  dislocation  density  in  the  matrix  metal,  p  given  by 


-  12feT  1 
p"b(l  -f)  t 


(3.3) 


where  f  is  the  volume  fraction  of  particulates,  ej  is  CTE  mismatch  strain  (ATAa),  b  is  Burgers’ 
Vector.  Once  p  is  calculated,  the  increase  in  the  matrix  flow  stress  due  to  CTE  mismatch  strain, 
AC>cte  may  be  estimated  by 


Aocte  = 


(3.4) 


where  P  is  a  constant  of  order  1,  p.  is  the  matrix  shear  modulus. 

Since  1986,  more  rigorous  dislocation  punching  models  have  been  developed  by  Taya  and 
his  co-workers,  and  they  are  aimed  at  various  types  of  reinforcement  geometry:  spherical  particle, 
short  fiber,  disc-shaped  filler  and  continuous  fiber.  In  this  paper,  we  shall  review  recent  studies  on 
dislocation  punching  models  for  these  types  of  reinforcement,  and  also  the  dislocation  punching 
model  for  thin  film/substrate  system. 


3.2  Dislocation  Punching  From  a  Short  Fiber 

Taya  and  Mori  [12]  studied  analytically  the  dislocation  punching  from  a  short  fiber  in  a  short 
fiber  metal  matrix  composite  (MMC).  A  short  fiber  is  simulated  by  a  prolate  ellipsoidal 
inhomogeneity  to  make  use  of  the  Eshelby’s  method  and  the  CTE  mismatch  strain  ex  is  initially 
adhered  to  the  fiber-matrix  interface  as  surface  dislocations  (prismatic  dislocation  loops).  Fig. 
3.1(a).  Due  to  the  fact  that  the  stress  at  fiber-ends  is  highest  and  the  punching  of  the  dislocation 
loops  are  favorable  along  the  fiber  axis  (x3-axis),  some  of  the  dislocation  loops  are  punched  along 
the  X3-axis,  Fig.  3.1(b). 
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Figure  3.1  Analytical  model  used  by  Taya  and  Mori  [12]  to  study  the  relaxation  of  CTE  mismatch 
strain  at  short  fiber/metal  matrix  interface  by  dislocation  punching:  (a)  before  punching  (unrelaxed 
stage),,  (b)  after  punching  (relaxed  stage) 


The  stress  field  in  and  around  a  short  fiber  and  the  total  strain  energy  of  a  short  fiber  MMC 
can  be  calculated  if  the  following  eigenstrains  are  given  in  the  domains  of  fiber  (Qi)  and  of  the 
punched  region  encompassing  the  fiber  (O2): 
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(3.6) 


where  c'  and  c  are  the  punching  distance  and  half  of  short  fiber  length  along  the  x3-axis, 
respectively.  Fig.  3.1.  The  criterion  for  dislocation  punching  is  given  by 
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(3.7) 


_  au  aw 

”  3c'  "  3c' 

where  U  is  the  total  strain  energy  of  the  composite  of  Fig.  3.1(b)  and  its  explicit  expression  is 
given  elsewhere  [12],  and  W  is  the  total  energy  dissipation  due  to  punching,  i.e.,  plastic  work 
required  for  the  motion  of  the  dislocation  loops  during  punching  and  given  by 


where  c/a  is  the  fiber  aspect  ratio,  k  is  the  friction  stress  of  the  matrix  metal  against  dislocation 
motion. 

Eq.  (3.7)  provides  the  relations  between  punching  distance  (c')  and  several  constituent 
parameters.  Fig.  3.2  shows  such  an  example,  i.e.  the  punching  distance  (c')  normalized  by  half  of 
fiber  length  (c)  as  a  function  of  fiber  aspect  ratio  (c/a).  Fig.  3.2  implies  that  the  punching  along  the 
fiber  axis  would  become  more  difficult  as  fiber  aspect  ration  increases.  In  this  case,  the  punching  is 
more  likely  to  take  place  along  the  direction  transverse  to  the  fiber  axis,  as  discussed  in  Section 
3.3.  The  relaxation  by  dislocation  punching  results  in  lowering  the  stress  field  in  around  a  short 
fiber.  Table  3.1  shows  the  normal  stress  components  (an,  G22  and  033)  normalized  by  2q  where 
(l  is  the  matrix  shear  modulus  for  relaxed  (after  punching)  stage,  Fig.  3.3(a)  and  non-relaxed 
(before  punching)  stage.  Fig.  3.3(b).  The  stresses  are  calculated  for  several  locations:  within  a 
fiber  (A),  just  outside  the  fiber  at  equator  (B)  and  just  outside  the  fiber-ends  (C).  Table  3.1 
indicates  that  the  magnitude  of  the  stress  field  in  and  around  a  short  fiber  in  the  relaxed  stage 
indeed  is  smaller  than  that  of  the  non-relaxed  stage.  The  data  used  to  compute  the  stress  field 
(Table  3.1)  and  punching  distance  (Fig.  3.2)  referred  to  SiC  short  fiber/Al  system  and  are  given 
elsewhere  [12]. 


Table  3.1.  Stress  at  points  A,  B  and  C  (see  Fig.  3) 


^xlO-2 

2|lXlU 

For  relaxed  state 

For  non-relaxed  state 

11 

22 

33 

11 

22 

33 

A 

-0.225 

-0.225 

-0.84 

-0.245 

-0.245 

B 

-0.225 

0.198 

0.069 

-0.245 

0.213 

C 

0.0035 

0.0035 

-0.84 

-0.149 

-0.149 
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Figure  3.2  Punching  distance  (c')  normalized  by  half  of  short  fiber  length  vs.  fiber  aspect  ratio 
(c/a)  for  SiC  short  fiber/Al  matrix  composite  with  AT=-200oC,  k/|i=0.001. 


(a)  (b) 

Figure  3.3  Stresses  are  calculated  at  points  A,  B,  and  C  for  (a)  relaxed  and  (b)  unrelaxed  stages. 
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Dunand  and  Mortenson  [13]  modified  the  above  model  to  account  for  partial  relaxation  by 
dislocation  punching  along  the  fiber  axis.  The  eigenstrains  based  on  this  model  are  given  by  Eqs. 
(3.5)  and  (3.6)  where  (3,3)  element  in  e*1  is  replaced  by  (l-%)  ex  and  that  in  e  *2  by  XeT^~j> 

where  %  denotes  the  fraction  of  the  misfit  strain  to  be  relaxed  by  dislocation  punching.  The 
modified  model  of  Dunand  and  Mortenson  appears  to  explain  well  their  experimental  results. 

3.3  Dislocation  Punching  From  a  Spherical  Particle 

The  dislocation  punching  from  a  spherical  particle  was  studied  both  theoretically  and 
experimentally  by  Taya  et  al.  (1991).  It  is  expected  that  the  direction  of  dislocation  punching  from 
a  spherical  particle  is  spherically  symmetric.  Fig.  3.4(a),  (b)  and  (c)  show  the  CTE  mismatch  strain 
adhered  to  the  particle-matrix  interface  as  dislocation  loops  before  punching,  after  punching,  and 
the  rearrangement  of  the  dislocation  loops,  respectively.  The  dislocation  punching  shown  by  Fig. 
3.4(b)  should  take  a  form  of  glide  motion  as  illustrated  by  Fig.  3.5,  i.e.  during  punching  an  edge 
dislocation  is  degenerated  to  two  paitial  edge  dislocations  which  move  along  the  glide  planes,  and 
after  punching,  two  edge  dislocations  are  formed  to  one  edge  dislocation  arrayed  in  a  spherically 
symmetric  manner.  The  eigenstrains  representing  the  case  of  Fig.  3.4(b)  (or  Fig.  3.4(c))  are 
isotropic  in  a  domain  of  punched  spherical  region  including  the  particle  (Q2)  ar>d  given  by 


■iyA  " 


e*  |  0  eT(%„ 


0  «T  (%"f 


The  criterion  for  dislocation  punching  is  given  by  Eq.  (3.7).  c'  and  c"  are  defined  by  Fig. 
3.4(b).  The  geometry  after  punching  in  a  spherical  particle  MMC  is  shown  schematically  in  Fig. 
3.6  where  (a)  and  (b)  denote  the  case  of  extensive  and  less  extensive  punching,  respectively.  It  is 
noted  here  that  in  the  case  of  extensive  punching,  Fig.  3.6(a),  the  punched  dislocations  located  at 
the  boundary  of  0.2  are  canceled  out  if  they  are  overlapped  with  adjacent  dislocations,  resulting  in 
no  contribution  of  geometrically  necessary  dislocations  to  increase  in  the  matrix  flow  stress. 
However,  as  the  geometrically  necessary  dislocations  are  punched,  statistically  stored  dislocations 
are  left  out  in  the  wake  region  (O2  -  Qi)  since  the  actual  mode  of  the  dislocation  punching  is  not 
spontaneous,  but  successive  punching  of  dislocation  loops.  Thus,  for  extensive  punching,  one  can 
expect  an  increase  in  matrix  flow  stress  due  to  the  statistically  stored  dislocations. 
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The  formulations  of  the  stress  and  strain  field,  and  also  the  strain  energy  of  a  composite,  are 
again  based  on  the  Eshelby’s  method  and  its  details  are  given  elsewhere  [5,14].  The  average 
dislocation  density  for  less  extensive  punching  mode,  Fig.  3.6(b)  is  given  by 

P=iRT^f)  <3-10> 

where  f  is  the  volume  fraction  of  particles,  ej  is  CTE  mismatch  strain,  a  is  the  radius  of  a  spherical 
particle  and  b  is  Burgers’  vector.  A  comparison  between  Eqs.  (3.3)  and  (3.10)  reveals  that  the 
averaged  dislocation  density  in  the  matrix  predicted  by  the  Arsehault  and  Shi  model  is  twice  as 
much  as  that  by  the  Taya  et  al.  model,  if  t  is  set  equal  to  a.  In  the  case  of  extensive  punching,  Fig. 
3.6(a),  the  average  dislocation  density  is  computed  by  Shibata  et  al.  (1992)  and  given  by 

p  =  ^Iln(f-1/3)  (3.11) 

It  is  noted  that  Eq.  (3.1 1)  was  derived  by  accounting  for  only  statically  stored  dislocations. 

In  order  to  predict  the  yield  stress  of  a  particle  MMC,  one  should  consider  another 
strengthening  mechanism,  i.e.  back  stress  strengthening  which  is  to  account  for  resistance  of 
elastic  particles  against  plastically  deforming  matrix  metal.  Hence  the  increase  in  the  composite 
yield  stress  over  the  unreinforced  metal  is  due  to  two  mechanisms:  dislocation  punching  as  a  result 
of  relaxation  of  CTE  mismatch  strain  at  the  interface  (AGcte)  and  back  stress  (AOb).  In  order  to 
test  the  validity  of  the  above  strengthening  mechanisms,  Taya  et  al.  [5]  conducted  experiment  on 
SiC  particle/6061  Al  composite-T4  (SiCp/Al).  The  mechanical  properties  of  SiCp/Al  composite  are 
given  in  Table  3.2. 


Table  3.2.  Material  Properties 


Parameter 

Unit 

Al  Matrix 

SiC  Particle 

Young’s  modulus 

GPa 

68.3 

427 

Yield  strength 

MPa 

97.0 

— 

Poisson’s  ratio 

1 

0.33 

0.17 

CTE 

x  10-^C 

23.6 

4.3 

Burgers’  vector 

nm 

0.283 

— 

Average  particle  radius 

|im 

— 

5.0 
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Figure  3.4  Dislocation  punching  model  used  by  Taya  et  al.  [5]:  (a)  CTE  mismatch  strain  ei 
adhered  to  the  particle/matrix  surface  as  prismatic  dislocation  loops,  (b)  after  punching  to  c',  (c) 
rearrangement  of  loops  for  spherical  symmetry. 
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Figure  3.5  Dissociated  dislocations  make  glide  motion  and  are  combined  to  form  a  new 
dislocation. 


(a)  (b) 

Figure  3.6  Punched  regions  in  a  particle  MMC:  (a)  for  extensive  punching  resulting  in  overlapping 

of  the  frontal  boundaries  and  (b)  for  less  extensive  punching. 


SiCp/Al  specimens  were  solution-treated  at  530°C,  thus  eliminating  much  of  the  dislocations 
induced  by  processing,  then  quenched  to  several  different  temperatures  (Tq):  room  temperature,  0, 
-64,  and  -196°C,  followed  by  tensile  testing  at  these  temperatures  to  measure  the  yield  stresses. 
The  purpose  of  quenching  is  to  generate  given  CTE  mismatch  strain  at  the  interface.  Fig.  3.7(a) 
shows  the  composite  yield  stresses  experimentally  obtained  (open  circle)  and  predicted  by  the 
dislocation  punching  and  backstress  strengthening  models  (filled  circle)  as  a  function  of 
temperature  change  by  quenching,  AT  =  Tq  -  530.  In  order  to  examine  as  to  which  strengthening 
mechanism  is  dominant,  AGb  and  AgCte  are  plotted  in  Fig.  3.7(b),  indicating  that  AgCte  is  larger 
than  AGb-  Namely  strengthening  due  to  dislocation  punching  is  dominant. 


(a)  (b) 

Figure  3.7  Composite  yield  stress  vs.  temperature  change  by  quenching:  (a)  a  comparison 
between  the  experiment  and  the  prediction  based  on  the  dislocation  punching  model,  (b) 
breakdown  of  the  increase  in  the  composite  yield  stress  indicating  strenthening  due  to  dislocation 
punching  (Agcte)  dominant  over  that  due  to  back  stress  (AGb). 

3.4  Dislocation  Punching  From  a  Continuous  Fiber 

As  the  results  of  Fig.  3.2  suggest,  the  dislocation  punching  from  a  continuous  fiber  in  a 
continuous  fiber  MMC  is  more  likely  to  take  place  along  the  transverse  direction  (transverse  to  the 
fiber  axis).  This  was  studied  analytically  by  Shibata  et  al.  [15].  The  eigenstrains  after  punching  are 
given  by 
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(3.12) 


(3.13) 


where  Qj  and  Q2  denote  the  domains  of  the  fiber  and  the  punched  region  encompassing  the  fiber, 
respectively;  a  and  R  are  the  radii  of  the  fiber  and  the  punched  region,  respectively.  By  using  the 
Eshelby’s  model  to  compute  the  total  elastic  energy  of  the  composite  and  the  criterion  for 
punching,  Eq.  (3.7),  and  the  dissipation  energy  for  punching,  one  can  obtain  the  punching 
distance  (radius  of  the  punched  region)  R  as 


f  (1  +  2v)  |ieT  1  !/2 
1  2  (1  -  v)  k  J 


(3.14) 


where  |i  and  k  are  the  shear  modulus  and  friction  stress  of  the  matrix  metal,  respectively.  Eq. 
(3.14)  was  obtained  for  a  single  continuous  fiber  embedded  in  an  infinite  matrix.  The  result  of  the 
punching  distance  R  for  a  continuous  fiber  MMC  with  finite  volume  fractions  of  fiber  can  be 
obtained  similarly  [15].  It  is  noted  here  that  the  yield  criterion  based  on  Hill’s  model 

I  Gr  -  oe  I  =  2k  (3.15) 

is  not  satisfied  within  the  plastic  domain  (O2  -  &  l)  where  Gr  and  Gq  are  radial  and  circumferential 
stress  components.  However,  the  average  stress  across  the  boundary  of  02  is  defined  by 

ay  =  ^  {  g jj  (R  -  e)  +  Gy  (R  +  e) }  /2  (3.16) 

satisfy  the  yield  condition,  Eq.  (3.15)  which  gives  rise  to  Eq.  (3.14). 
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3.5  Dislocation  Punching  From  a  Planar  Interface 


The  generation  of  dislocations  at  the  interface  of  a  thin  coating/substrate  system  has  been 
studied  in  terms  of  critical  thickness  of  the  coating-misfit  strain  relation  [16,17],  These  studies  are 
aimed  at  epitaxial  layers  and  the  order  of  the  critical  thickness  is  normally  much  smaller,  10  ~  100 
nm.  For  a  coated  substrate  system  with  non-epitaxial  interface  and  intermediate  coating  thickness, 
the  results  of  the  above  studies  are  not  applicable. 

Here  we  consider  a  metal  thin  coating/ceramic  substrate  with  CTE  driven  misfit  strain  at  the 
interface.  The  CTE  misfit  strain  is  defined  by 


e*  =  (occ  -  as)AT 


(3.17) 


where  Oc  and  as  are  CTE  of  the  coating  and  substrate,  respectively,  8ij  is  knocker’s  delta  and  AT 
is  temperature  change.  First  we  examine  the  stress  field  before  punching  Fig.  3.8(a)  and  then  the 
case  of  punching  to  some  distance,  Fig.  3.8(b). 


-L  -L  -L.  J _ L 

metal 

ceramic 

(a)  (b) 

Figure  3.8  Dislocation  punching  model  for  a  thin  metal  coating/ceramic  substrate  system:  (a) 
before  punching,  (b)  after  punching  to  distance  z. 
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It  is  assumed  in  this  model  that  the  thickness  of  coating  is  reasonably  thin,  so  that  out-of¬ 
plane  stress  component  in  the  coating  (033)  can  be  ignored,  and  also  that  the  plane  strain  condition 
be  held  for  in-plane  directions,  xi-  and  x2-axes  within  the  coating. 

Before  Punching 

By  the  assumption  of  plane  strain  along  the  xj-  and  x2-axes, 


en  =e22  =  -e 


(3.18) 


where  ejj  is  elastic  strain  component.  Hooke’s  law  with  033  =  0  gives  rise  to 


.  _£ll  v 
el  1  —  £  ~  £  ^22 


022  V 

622  ~  £  ~  £  Oi  l 


(3.19) 


e33  =  -  £  (Oil  +  O22) 

where  v  is  the  Poisson’s  ration  of  the  coating. 

From  Eqs.  (3.18)  and  (3.19)  we  obtain  the  in-plane  stresses  and  out-of-plane  strain  as 


Oil  =  022  =  - 


E  e 


(1  -v) 


(3.20) 


v  * 

633  =  “ - e 

1  -  v 


Punching  to  Distance  z 

Assuming  that  the  dislocation  punching  took  place  with  its  frontal  boundary  located  Z  from 
the  interface.  Fig.  3.8(b),  we  shall  calculate  the  stress  and  strain  field  within  the  punched  region 
(Ql)  and  elastic  region  (Qq)  and  also  the  condition  for  punching. 
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in  Elastic  Region  (C2q): 


Stress  Gij  and' strain  ey  are  the  same  as  those  in  the  coating  before  punching,  i.e.  given  by 
Eq.  (3.20). 


in  Plastic  Region  (Qi): 


Denoting  plastic  strain  in  the  punched  region  by  e|j  and  using  the  plane  strain  condition  in  the 
xi  -  X2  plane  and  incompressibility  of  plastic  strain,  one  arrives  at 


(3.21) 


and 

ey  =  0 

It  is  noted  that  stresses  in  Qi  are  all  zero  due  to  complete  relaxation. 
Computation  of  Elastic  Energy  U  and  Plastic  Work  W 

The  elastic  energy  per  unit  width  along  the  X2*axis  is  given  by 

U  =  j  ay  ey  (£  -  z) 

A  substitution  of  Eqs.  (3.18)  and  (3.20)  into  (3.23)  yields 

Eel  (l-  z) 

U  =  - 1 - 

.  ( 1  —  v ) 

The  plastic  work  per  unit  width,  W  is  given  by 

W  =  Gjj  ep.  z 
J  y 


(3.22) 


(3.23) 


(3.24) 


(3.25) 


44 


by  using  Eq.  (3.21)  and  on  =  022  =  <7y  (yield  stress),  W  is  reduced  to 


W  =  4kejz 


where  Gy  =  2k  was  used  and  k  is  yield  stress  in  shear.  The  punching  criterion  is  given  by 


au  aw 
az  "  az 


(3.26) 


where  the  left  term  represents  the  driving  force  for  punching  while  the  right  is  the  retarding  force 
against  punching.  By  substituting  Eqs.  (3.24)  and  (3.25)  into  (3.26)  we  obtain  the  condition  for 
punching 


ey  > 


4k  (1-v) 
E 


(3.27) 


Inequality  (3.27)  can  be  derived  by  invoking  the  yield  criterion  across  the  frontal  boundary  of 
the  punched  region,  i.e.  the  interface  between  Qq  and  Qj,  Fig.  3.8(b).  The  average  stresses 
across  the  boundary  are  calculated  by  using  Eqs.  (3.16),  (3.20)  and  (3.22). 


<711  =<722=- 


E  ej 
2(1  -v) 


(3.28) 


which  is  substituted  into  Hill’s  yield  criterion 

I  033-011  I  >  2k  (3.29) 

A  substitution  of  Eq.  (3.28)  and  033  =  0  into  (3.29)  leads  to  (3.27). 
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4.  Analytical  Study  II:  Dimensional  Change  of  MMCs 


4.1  Analytical  Model  for  Short  Fiber  Composites 

4.1.1  Introduction 

Dimensional  stability  in  high  temperature  use  environments  is  crucial  for  many  metal  matrix 
composite  (MMC)  applications.  While  MMCs  reinforced  by  ceramic  fibers  offer  creep  resistance 
superior  to  the  unreinforced  metal,  it  is  well  known  that  the  combination  of  relatively  small  applied 
stresses  and  cyclic  thermal  loading  can  result  in  creep  rates  far  in  excess  of  those  observed  during 
isothermal  high  temperature  creep.  Even  in  the  absence  of  an  applied  stress,  a  detrimental 
dimensional  change  has  been  observed  in  a  number  of  thermal  cycled  metal  matrix  composite 
systems.  These  effects  are  attributed  to  internal  stresses  that  are  developed  during  thermal  cycling 
due  to  the  coefficient  of  thermal  expansion  mismatch  between  the  matrix  and  fibers.  Creep  and 
thermal  cycling  creep  are  among  the  least  understood  aspects  of  the  deformation  of  MMCs,  but  are 
of  utmost  importance  for  the  development  of  MMCs  for  high  temperature  structural  applications. 

Creep  in  short  fiber  metal  matrix  composites  has  been  studied  both  experimentally  and 
theoretically  by  numerous  researchers.  A  comprehensive  review  is  not  attempted  here,  but  the 
reader  is  referred  to  recent  reviews  by  Taya  [1]  and  Taya  and  Arsenault  [2],  Recent  works 
regarding  modeling  of  creep  of  short  fiber  composites  include  Zhu  and  Weng  [3-5],  Wang  and 
Weng  [6],  and  Pan  and  Weng  [7].  In  these  works,  Eshelby's  equivalent  inclusion  method  [8]  was 
combined  with  the  Mori-Tanaka  mean  field  approach  [9]  to  estimate  the  stress  redistribution  due  to 
an  incremental  creep  of  the  metal  matrix  and  the  subsequent  creep  deformation  of  the  composite.  A 
similar  approach  was  used  by  Taya  and  Mori  [10]  to  model  creep  of  a  short  fiber  MMC  during  the 
high  temperature  hold  stage  of  a  cyclic  thermal  loading.  With  regards  to  thermal  cycling  and 
thermal  cycling  creep  of  short  fiber  MMCs,  Derby  [11]  and  Taya  et  al.  [12]  have  recently  provided 
a  comprehensive  review  of  the  literature  and  a  list  of  references  with  regards  to  both  experimental 
and  theoretical  aspects  can  be  found  therein.  Other  recent  contributions  have  been  directed  at 
understanding  the  effects  of  reinforcement  geometry  on  the  thermal  cycling  creep  of  short  fiber 
MMCs  through  the  development  of  idealized  models  and  the  application  of  finite  element  methods 
[13-15]. 

During  isothermal  creep  of  a  MMC,  degradation  of  the  fiber  matrix  interface,  often  in  the 
forms  of  void  nucleation  and  growth  followed  by  debonding,  has  been  observed  to  lead  to  tertiary 
creep  and  creep  fracture  [16-19].  The  accumulation  of  inelastic  strain  during  thermal  cycling  has 
also  been  observed  to  be  accompanied  by  the  development  of  interfacial  damage  in  the  forms  of 
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porosity,  debonding,  and  an  attached  reaction  zone  containing  radial  cracks  [20].  These  damage 
processes,  in  both  isothermal  and  thermal  cycling  creep  further  accelerate  the  inelastic  strain 
accumulation  and  thus  it  is  important  to  understand  the  evolution  of  damage  and  its  effect  on  the 
overall  response  of  the  MMC. 

This  work  is  a  continuation  of  the  authors'  study  of  creep,  thermal  cycling,  and  creep  thermal 
cycling  (thermal  cycling  with  an  applied  stress)  of  metal  matrix  composites.  In  particular,  a 
micromechanics-based  analytical  approach  is  proposed  to  model  creep  and  thermal  cycling  creep  of 
short  fiber  MMCs  with  interfacial  damage.  The  approach  is  an  extension  of  the  analyses  of  Taya  et 
al.  [12]  for  creep,  Taya  and  Mori  [10]  for  thermal  cycling,  and  Dunn  and  Taya  [21]  for  creep 
thermal  cycling  of  MMCs  and  is  based  on  the  combination  of  Eshelby's  equivalent  inclusion 
method  [8]  and  the  Mori-Tanaka  mean  field  approach  [9]  (see  for  example  [22]).  Interfacial 
damage  is  modeled  through  a  hybrid  approach  where  the  composite  is  assumed  to  contain  both 
perfectly  bonded  fibers  and  fibers  with  damaged  interfaces.  The  evolution  of  damage  is  described 
by  an  arbitrary  statistical  distribution.  Efforts  to  correlate  the  statistical  distribution  of  interfacial 
damage  with  experiment  are  currently  underway  and  the  results  will  be  subsequently  reported. 

4.1.2  Analytical  Model 

It  is  assumed  that  initially  all  fibers  are  perfectly  bonded  to  the  matrix.  At  some  point  during 
the  creep  or  thermal  cycling  creep  loading,  though,  interfacial  damage  begins  to  occur.  The  effect 
of  the  damage  is  a  loss  in  the  ability  to  effectively  transfer  load  from  the  matrix  to  the  fibers.  When 
damage  has  occurred,  it  is  assumed  that  the  aligned  short  fiber  and  matrix  are  still  in  mechanical 
contact  except  at  and  near  the  poles  of  the  fibers.  This  characterization  is  based  on  numerous 
reports  from  the  literature  regarding  debonding  during  creep  and  a  microscopy  study  of  thermal 
cycled  W/FeCrAlY  short  fiber  composites  of  which  a  typical  micrograph  is  shown  in  Fig.  4.1 
[23].  The  fiber-matrix  interface  is  assumed  to  be  stress  free  in  the  x3  direction  but  not  in  the  Xj  and 
Xo  direction  except  in  the  vicinity  of  the  fiber  ends.  This  assumption  seems  reasonable  when  the 
interface  between  the  short  fiber  is  in  tension,  but  may  not  be  as  valid  when  the  stress  field  at  the 
interface  is  compressive.  These  implications  will  be  further  discussed  in  the  following  section. 
The  stress  field  in  the  debonded  fibers  is  analyzed  by  following  the  approach  proposed  by  Tanaka 
et.  al.  [24]  to  model  the  stress  distribution  in  a  spherical  particle  after  cavitation.  The  stress  field  in 
the  debonded  fibers  is  simulated  by  that  of  a  fictitious  fiber  with  anisotropic  moduli  of  Cijkl  =  0 

with  the  exception  of  Cdnn  =  C2222  and  cti22  =  C22if  lt  is  noted  that  this  reduction  in  moduli 
is  only  valid  for  the  case  considered  here,  aligned  short  fibers  subjected  to  uniaxial  mechanical 
loading  along  the  axis  of  the  fiber. 


48 


In  this  section,  the  response  of  a  composite  containing  two  types  of  short  fibers  (perfectly 
bonded  and  debonded)  will  be  obtained.  The  debonded  fibers  are  assumed  to  be  described  as 
discussed  in  the  previous  section  and  the  criterion  for  progressive  debonding  will  then  be 
discussed  in  the  following  section.  The  domain  of  the  analytical  model,  shown  in  Fig.  4.2, 
consists  of  an  infinite  metal  matrix  containing  a  random  distribution  of  both  perfectly  bonded  and 
debonded  aligned  short  fibers.  The  domains  of  the  entire  composite  and  of  the  bonded  and 
debonded  fibers  are  denoted  by  D,  Qb,  and  Qd  respectively,  and  that  of  the  matrix  is  thus  denoted 
by  D  -  Qb  -  Qd.  The  elastic  moduli  of  the  matrix  and  two  types  of  fibers  are  Cm,  Cb,  and  Cd, 
respectively.  It  is  noted  that  Cm  and  Cb  are  defined  in  the  standard  manner  while  Cd  has  been 
defined  in  Section  2.  In  this  work,  bold  and  underscored  symbols  represent  tensorial  quantities. 
Both  types  of  fibers  are  modeled  as  ellipsoidal  inhomogeneities  (prolate  spheroids)  and  for 
simplicity  both  the  matrix  and  fibers  are  assumed  to  be  initially  isotropic  in  stiffness  and  thermal 
expansion.  To  model  the  creep  response  of  the  composite,  a  constant  applied  stress,  a0,  is 

assumed  and  in  addition  the  idealized  time-temperature  thermal  cycle  shown  in  Fig.  4.3  is  assumed 
during  thermal  cycling  creep.  It  is  assumed  that  the  fibers  deform  elastically  and  the  matrix  is 
capable  of  elastic/plastic/creep  deformations.  The  material  properties,  with  the  exception  of  the 
yield  stress  and  creep  properties  of  the  matrix,  are  assumed  independent  of  temperature.  This 
assumption  can  be  relaxed,  however,  and  the  analysis  can  be  easily  implemented  in  an  incremental 
manner. 

Creep  of  Short  Fiber  Metal  Matrix  Composites 
Response  to  the  Applied  Stress 

The  overall  response  of  the  composite  to  the  uniform  applied  stress,  a0,  (with  corresponding 

strain  e°)  can  be  estimated  by  the  use  of  Eshelby's  equivalent  inclusion  method  [8]  coupled  with 
the  Mori -Tanaka  mean  field  approach  [9]  as  applied  to  a  three-phase  composite  [25].  Through  the 
equivalent  inclusion  method,  the  stress  in  a  single  representative  perfectly  bonded  fiber  is  given  by: 

<j°  +  cfi  =  Cb-[e°  +  e  +  eb]  =  Cm-[e°  +  e  +  eb  -  e*^]  (4.1) 


That  in  a  single  representative  debonded  fiber  is  given  by: 


49 


Area  of  Damage  (square  microns) 


l 


Figure  4.2  Analytical  model. 


Figure  4.3  Idealized  temperature-time  relation 


(4.2) 


g°  +  (3d  =  Cd-[e°  +  e  +  ed]  =  Cm-[e°  +  e  +  ed  -  e*°d] 

where  e  is  the  volume  averaged  disturbance  strain  in  the  matrix  and  e1  and  e*ai  are  the  disturbance 
strain  and  fictitious  eigenstrain  in  the  ith  domain  (i=b,d)  required  for  the  simulation  of  the 
inhomogeneities  by  the  equivalent  inclusion,  respectively.  These  can  be  expressed  as: 

<£>m=Cm-e  (4.3) 

e>  =  S-e*0’  '  (4.4) 

where  S  is  Eshelby's  tensor  [8]  and  <  a  >  m  is  the  volume  average  of  the  disturbance  stress  (due  to 

both  £2b  and  Qd)  in  the  matrix.  Eshelby's  tensor  is  a  function  only  of  the  shape  of  the 
inhomogeneity  and  the  elastic  moduli  (Poisson's  ratio  for  isotropic  materials)  of  the  matrix  and 
thus  is  the  same  for  both  the  bonded  and  debonded  fibers.  It  is  not  difficult  to  show  that  the 
volume  average  of  c*  over  the  entire  composite  vanishes  which  results  in: 

e  =  (S  -  I>[fb  e*^  +  fd  e*°d]  (4.5) 

where  f  is  the  volume  fraction  of  fibers.  After  some  manipulation,  eqs.  (4.1)  -  (4.5)  yield  <  o  >m 

which  can  be  used  to  obtain  the  average  stresses  and  strains  in  each  phase  and  in  the  composite. 

c° 

Of  particular  interest  are  the  average  stress  in  the  matrix ,  <om,  and  the  average  strain  of  the 
a° 

composite  ,  <e>  c  ,  which  can  be  expressed  as: 

<o>°m°  =  R(Cm,  Cb,  Cd,  S,  fb,  fd)-tf>  (4.6) 

<e>ac°  =  .  H(Cm,  Cb,  Cd,  S,  fb,  fd)- <f  (4.7) 


52 


Response  to  an  Incremental  Creep 


An  incremental  creep  of  the  matrix,  dec,  results  in  a  redistribution  of  stress  from  the  creeping 
matrix  to  the  elastic  fibers,  but  as  noted  by  Weng  [26],  is  itself  a  stress-free  process.  Thus, 
Eshelby’s  equivalent  inclusion  method  can  be  used  to  estimate  the  stress  redistribution  due  to  an 
incremental  creep: 


dab  =  Cb  ’[de  +  deb  +  dec  ]  =  Cm-[de  +  deb  -  de*Cb]  (4.8) 

dcd  =  Cd  -[de  +  ded  +  dec]  =  Cm-[de  +  ded  -  de*Cd]  (4.9) 


dec  and  de*c  1  are  the  incremental  creep  of  the  matrix  and  the  corresponding  incremental  eigenstrain 
and  all  other  stress  and  strain  increments  are  analogous  to  their  previous  definitions  in  eqs.(4.1) 
and  (4.2).  By  a  development  similar  to  that  used  to  obtain  eqs.  (4.6)  and  (4.7),  the  average  stress 
increment  in  the  matrix  and  strain  increment  of  the  composite  due  to  an  incremental  creep  of  the 
matrix,  dec,  can  then  be  expressed  as: 

<do>!n  =  U(Cm,  Cb,  Cd,  S,  fb,  fd)-dec  (4.10) 


ec 

<de>c 


V(Cm,  Cb,  Cd,  S,  fb,  fd)-dec 


(4.11) 


where  due  to  the  assumption  of  incompressibility  and  symmetry  of  deformation,  dec  is  given  by 
dec  =  [  -  d£c/2,  -  dec/2,  d£c,  0,  0,  0  ]T. 

For  simplicity,  the  creep  rate  of  the  matrix  is  assumed  to  be  described  by  the  Mises  J2  creep 
flow  rule  Odqvist  [27]: 


®Fj  =  k  J2  <  CT-j  >m  or  £c  =  ±A  I  am  ln  (4.12) 

where  d£C  =  £Cdt.  In  eq.  (4.12),  J2,  <  g-j  >m,  and  om  are  the  second  invariant,  deviatoric 

components,  and  the  flow  stress  of  the  average  stress  in  the  matrix  and  are  all  computed  in  the 
standard  manner  and  as  usual  the  (+)  or  (-)  sign  corresponds  to  the  sign  of  the  flow  stress.  The 
second  of  eqs.  (4.12)  results  when  axisymmetric  deformation  is  present  as  will  be  assumed  here. 
During  creep  of  the  matrix,  the  average  stress  in  the  matrix  can  be  determined  by  integration  of  eq. 
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(4.12)  with  the  initial  condition  of  eq.  (4.6).  From  the  average  stress  in  the  matrix,  the  flow  stress 
is  determined  as  a  function  of  £c,  which  in  turn  is  a  function  of  time.  The  flow  stress  can  then  be 
substituted  into  the- power  law  of  eq.  (4.12)  to  yield  a  first  order  nonlinear  ordinary  differential 
equation  in  ec  of  the  form: 

£c  =  A[(aa°  +  bec)]n  (4.13) 

For  the  applied  uniaxial  stress  along  the  fiber  (x3)  axis,  the  solution  to  the  ODE  of  eq.  (13)  is: 

£c(t)  { (a  o°)(1'n^  +  (1  -  n)  Abt}  (  (4.14) 

where  a  and  b  are  functions  of  Cm,  Cb,  Cd,  S,  fb,  and  fd.  It  is  noted  that  n  is  assumed  to  be 
independent  of  stress.  This  assumption  can  easily  be  relaxed  to  include  the  functional  dependence 
of  n  on  stress,  however  the  resulting  differential  equation  must  then  be  integrated  numerically. 
Finally,  the  creep  strain  of  the  composite  is  related  to  that  of  the  matrix,  £c(t),  through  eq.  (4.1 1). 

Progressive  Debonding 

The  model  developed  thus  far  is  valid  for  an  arbitrary  combination  of  short  fibers  with 
bonded  and  damaged  interfaces.  The  goal  of  this  approach  is  to  predict  the  progression  of  damage 
and  its  subsequent  effect  on  the  behavior  of  the  composite.  In  other  words,  it  is  desired  to 
understand  how  fd,  which  is  initially  zero,  evolves  with  the  accumulation  of  inelastic  strain.  Here 
an  attempt  is  made  to  relate  fd  to  the  strength  of  the  fiber-matrix  interface  through  a  statistical 
distribution  function,  the  form  of  which  will  be  obtained  from  careful  experiment  Heie  it  is 
assumed  that  the  evolution  of  fd  can  be  described  by  the  use  of  a  three-parameter  Weibull 
distribution  as: 


fd  -  F(^ib)  ftot 


(4.15) 


where  ftot  is  the  total  volume  fraction  of  fibers,  i.e.,  ftot  =  fb  +  fd  and  F(Cjb)  is  the  cumulative 
density  function: 


F(£tot)  =  1  -  exp 


fo  -  Cjb 
e  -  °ibj 


G><7ib 


(4.16) 


where  b,  oib,  and  0  are  the  shape  parameter,  guaranteed  value  of  a  (the  stress  O33  at  the  interface), 
and  the  scale  parameter  respectively.  Thus,  fd  is  a  function  of  (733  which  is  in  turn  a  function  of 
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time  as  G33  increases  as  stress  is  redistributed  from  the  creeping  matrix  to  the  elastic  fibers. 
Equation  (4.13)  is  still  valid,  but  now  it  must  be  integrated  numerically  as  a  and  b  are  now 
functions  of  time.  Although  a  Weibull  distribution  is  used,  here,  it  is  emphasized  that  the  proposed 
model  is  easily  able  to  handle  an  arbitrary  statistical  distribution. 

Thermal  Cycling  Creep  of  Short  Fiber  Metal  Matrix  Composites 

In  this  section  the  development  of  an  approach  to  model  thermal  cycling  creep  of  short  fiber 
metal  matrix  composites  is  outlined.  The  model  extends  that  of  Dunn  and  Taya  [21]  for  thermal 
cycling  creep  of  undamaged  short  fiber  composites  to  model  the  effects  and  progression  of 
interfacial  damage  in  the  same  manner  as  was  done  with  regards  to  isothermal  creep  in  the  previous 
section.  Thus,  only  a  brief  outline  of  the  approach  will  be  given  and  the  reader  is  referred  to  Dunn 
and  Taya  [21]  for  additional  details. 

The  composite  subjected  to  thermal  cycling  creep  is  assumed  to  be  subjected  a  constant 
applied  stress  and  then  the  idealized  thermal  cycle  shown  in  Fig.  4.3.  The  initial  response  to  the 

applied  stress  at  Aj  of  Fig.  4.3  is  identical  to  that  described  in  Section  4.2.1.  During  cooling  (At 

to  B^,  the  composite  is  subjected  to  the  uniform  temperature  drop  AT  =  TL  -  TH  which  induces  a 
non-uniform  thermal  stress  field  in  the  composite  due  to  the  mismatch  in  coefficients  of  thermal 
expansion  (CTEs)  between  the  matrix  and  fibers.  By  use  of  the  equivalent  inclusion  method,  the 
inhomogeneities  with  thermal  strain  eATcan  be  simulated  by  equivalent  inclusions  with  fictitious 
eigen  strains  e*AT  to  yield: 

Gb  =  Cb  ’[®"  +  eb  ‘  eAT  ]  =  Cm'[e  +  eb  -  e*ATb]  (4.17) 

2d  =  cd  [e  +  ed  -  eAT  ]  =  Cm*[e  +  ed  -  e*ATd]  (4.18) 

where  eAT  =  (ab  -  am)  AT  is  a  result  of  the  mismatch  between  the  CTEs  of  the  matrix,  am,  and  the 
bonded  and  debonded  fibers,  ab  =  ad . 

In  a  manner  similar  to  that  used  to  obtain  eqs.  (4.6)  and  (4.7),  the  average  thermoelastic 
AT  AT 

stress,  <G>  m  ,  and  strain,  <e>  c  ,  due  to  eATcan  be  obtained  as: 

AT 

<g>  m  =  P(Cm,  Cb,  Cd,  S,  fb,  fd,  ab,  aj  AT  (4. 19) 

AT 

<e>  c  =  Q(Cm,  Cb,  Cd,  S,  fb,  fd,  ab,  aj  AT  (4.20) 
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It  is  assumed  that  the  temperature  drop  from  TH  to  TL  is  sufficient  to  initiate  uniform  yielding 
of  the  matrix  metal.  The  critical  temperature  change  at  which  yielding  of  the  matrix  begins,  ATP  = 
Th  -  Tyl,  can  be  found  by  subjecting  the  average  thermoelastic  stress  in  the  matrix  to  the  yield 
criterion: 


3/2  <  G  ’>m-<  g 


(4.21) 


where  <  (J  '>m  is  the  deviatoric  portion  of  the  thermoelastic  stress  and  0YL  is  ^e  yield  stress  of  the 
matrix  at  TL.  When  the  applied  stress  is  uniaxial  along  the  x3-axis,  the  critical  temperature  change 
required  for  yielding  of  the  metal  matrix  can  then  be  expressed  as: 


AT 


cr 


gYL  ~  (^33  ~  ^3l)  a° 

P3-Pl 


(4.22) 


The  temperature  drop  from  TYL  to  TL  will  then  result  in  plastic  deformation  of  the  matrix. 

To  compute  the  plastic  deformation,  it  is  assumed  that  the  temperature  drop  from  T^  to  TL 
induces  a  uniform  plastic  strain,  ePL,  in  the  matrix  and  that  the  matrix  is  a  non-hardening  material. 
The  assumption  of  a  uniform  plastic  strain  in  the  matrix  neglects  the  effects  of  microyielding  that  is 
known  to  occur  near  the  fiber-matrix  interface  upon  relatively  low  temperature  changes.  It  is  the 
average  stress  field  in  the  matrix,  however,  that  is  believed  to  influence  the  macroscopic 
deformation  of  the  composite  and  this  is  the  basis  for  the  proposed  approach.  The  uniform  plastic 
strain,  ePL,  is  assumed  to  satisfy  the  incompressibility  requirement  e^L  =  0  and  as  a  result  of 
symmetry  (of  both  loading  and  microstucture)  can  be  expressed  as  ePL  =  ePL  [-1/2,- 1/2,  1,  0,  0,  0 
]T  =  ePL  K.  The  stress  field  due  to  ePLcan  be  computed  by  use  of  Eshelby's  equivalent  inclusion 
method  where  eAT  and  e*ATiin  eqs.  (4.17)  and  (4.18)  are  replaced  by  -ePLand  e*PL  1  respectively. 
The  plastic  strain,  ePL,  can  be  determined  by  applying  the  yield  criterion  to  the  unknown  total  stress 
existing  at  TL  (B  t).  With  the  plastic  strain,  ePL,  known,  the  average  stress  in  the  matrix  and  strain 
of  the  composite  at  can  be  expressed  as: 

d  rrO  /XT  0^ 

<<7>rn  =  <0>m  +  <0>  m  +  <cr>  m  =  R’g°  +  P'AT  +  U-£PL  (4.23) 

R  n°  AT  ePL 

<e>c1  =  <e>m  +  <e>m  +  <e>  m  =  H- a°+ Q-AT  +  V-ePL  (4.24) 


where  U  and  V  are  defined' by  eqs.  (4.10)  and  (4.11).  During  the  low  temperature  plateau,  B j-* 
Cj,  no  further  deformations  are  assumed  to  occur  and  thus  eqs.  (4.8)  and  (4.9)  describe  the  state 
of  the  composite  at  Q. 
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The  thermoelastic  stress  developed  during  the  heating  process,  Cj— >  Dj,  is  equal  in 
magnitude  but  of  opposite  sign  to  that  generated  upon  cooling.  Reverse  yielding,  resulting  in 
plastic  strains  ePH,  may  or  may  not  occur  depending  on  the  values  of  aYL,  the  yield  stress  at  TH 
(<7YH),  ar,d  AT.  The  state  of  the  composite  at  Dj  can  be  determined  in  the  same  manner  as  that 
upon  cooling  to  yield: 


<a>m1  =  Rg°  +  U-[ePL  +  ePH] 


<e>^  =  H-  g°+  V-[epl  +  ePH] 


(4.25) 

(4.26) 


During  the  high  temperature  plateau,  Dj— >  A2,  it  is  assumed  that  the  flow  stress -existing  at 
Th  (Dj)  is  further  relaxed  by  bulk  creep  of  the  matrix.  For  simplicity,  it  is  assumed  that  the  creep 
properties  of  the  composite  matrix  are  the  same  as  those  of  the  unreinforced  matrix  alloy. 
Eshelby’s  equivalent  inclusion  method  can  be  used  to  estimate  the  stress  redistribution  due  to  an 
incremental  creep  in  the  same  manner  as  described  in  Section  4.2. 1 .  Again  the  creep  rate  of  the 
matrix  is  assumed  to  be  described  by  eq.  (4. 12)  which  is  integrated  in  the  same  manner  outlined  in 
Section  4.2.1  but  now  with  the  initial  conditions  of  eq.  (4.26)  to  provide  the  creep  strain  at  the  end 
of  the  high  temperature  plateau  (t  =  x).  Once  ec(x)  is  obtained  the  overall  state  of  the  composite  is: 

R-  o°  +  U-[ePL  +  eph  +  ec(x)]  (4.27) 

H-  o°+  V-  [ePL  +  ePH  +  ec(x)]  (4.28) 

The  stress-strain-time  histories  during  the  second  and  subsequent  thermal  cycles  are  computed  in 
the  same  manner  as  that  of  the  first  thermal  cycle  except  for  the  provision  of  a  possible  non- zero 
residual  stress  at  A2  which  is  easily  accommodated.  It  turns  out  that  the  flow  stress  at  the  end  of 
each  subsequent  thermal  cycle  is  then  the  same  as  that  at  the  end  of  the  first. 

Finally,  a  statistical  distribution  can  be  implemented  to  describe  the  progression  of  damage  in 
a  manner  similar-  to  that  in  Section  4.2.1.  For  the  calculations  presented  in  the  following  section,  a 
three-parameter  Weibull  distribution  is  again  used  where  here  the  cumulative  density  function  is 
assumed  to  be  a  function  of  the  accumulated  inelastic  strain.  It  is  reiterated  that  the  form  of  the 
statistical  distribution  has  only  been  assumed  and  its  precise  form  will  be  determined  from  careful 
experimental  efforts  which  are  presently  underway.  It  is  noted  that  the  fiber-matrix  interface  is  in 
compression  during  cooling  and  tension  during  heating  [28].  Thus  as  damage  evolves,  it  is 
probably  more  appropriate  to  use  ^tot -  f b  +  f d  during  the  heating  portion  of  the  thermal  cycle  and 
ftot  =  fb  during  the  cooling  portion. 


A, 


<o>, 


<e> 


A2 
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4.1.3  Sample  Calculations 


To  illustrate,  numerical  computations  have  been  performed  using  material  properties  typical 
of  SiC  whisker  reinforced  6061-T6  A1  as  given  by  Morimoto  et  al.  [17].  Unless  otherwise  noted, 
the  creep  thermal  cycling  parameters  used  are:  2x  =  200  s,  AT  =  350K,  c°  =  10  MPa,  a  (aspect 
ratio)  =  5,  and  ftot  =  .2.  The  material  properties  used  are  EA1  =  67.6  GPa,  vA1  =  .33,  aA1  = 
24.7xl0-6  IK,  ESiC  =  427  GPa,  vSiC  =  .17,  <xA1  =  4.3xl0'6  [K,  Aa1  =  1.26xl0-16  Pa"/S,  nA1  = 
4.88. 

Fig.  4.4  shows  the  predicted  creep  strain  versus  time  for  a  total  volume  fraction  of 
reinforcement,  ftot  -  .15,  at  various  levels  of  debonding,  ftot-  .15.  It  is  seen  that  even  though  the 
matrix  is  assumed  to  be  in  steady-state  creep,  the  composite  exhibits  a  period  of  transient  creep.  It 
is  readily  apparent  that  the  presence  of  debonded  fibers  dramatically  affects  the  creep  behavior  of 
the  MMC.  Also  shown  in  Fig.  4.4  are  experimental  results  for  a  SiC  whisker/606 1-T6  Al  short 
fiber  composite  which  fall  between  the  predictions  for  significant  debonding.  The  analytical 
prediction  for  the  stress  at  the  fiber-matrix  interface  at  time,  t=0,  is  =  257  MPa  which 
significantly  exceeds  the  ultimate  strength  of  unreinforced  6061-T6  Al  at  573K.  It  is  thus  likely 
that  significant  debonding  existed  from  the  onset  which  agrees  with  the  predictions  of  the  proposed 
model.  In  Fig.  4.5,  the  capabilities  of  the  present  model  with  regards  to  the  effects  of  fiber-matrix 
debonding  are  illustrated.  Fiber-matrix  debonding  is  assumed  to  be  described  by  the  Weibull 
statistics  of  eqs.  (4.15)  and  (4.16)  where  it  is  assumed  that  o^0  =  180  MPa  (roughly  twice  the 
ultimate  strength  of  6061-T6A1  at  573K)  corresponding  to  a  fairly  strong  bond.  The  effects  of  the 
Weibull  shape  parameter  on  the  creep  strain  of  the  composite  are  shown  in  Fig.  4.5  where  the 
corresponding  Weibull  probability  and  cumulative  de'nsity  functions  are  shown  in  the  inset.  It  is 
readily  seen  that  the  proposed  approach  allows  simulation  of  the  entire  creep  strain  vs.  time  curve 
of  an  aligned  short  fiber  composite  with  the  exception  of  the  creep  rupture  strain. 

In  Fig.  4.6,  the  accumulated  strain  as  a  function  of  the  number  of  thermal  cycles  is  shown  for 
various  values  of  fy.  for  a  SiC/Al  MMC  subjected  to  thermal  cycling  creep  loading.  The  highly 
nonlinear  response  as  a  function  of  f^  is  clearly  evident.  The  progression  of  damage  appears  as  a 
jumping  from  the  bottom  line  in  Fig.  4.6  (fj  =  0)  up  to  higher  values  of  with  increasing  N 
resulting  in  an  overall  nonlinear  response  of  the  composite.  This  is  in  qualitative  agreement  with 
experimental  findings  of  Yoda  et  al.[29]  and  Armstrong  et  al.  [23]  for  thermal  cycling.  The 
predicted  accumulated  strain  of  the  composite  subjected  to  creep  only  at  Tjj,  based  on  the  present 
theory  (eqs.  (4.8)  -  (4.14))  is  far  less  than  those  under  thermal  cycling  creep  loadings,  thus  cannot 
be  plotted  in  Fig.  4.6  (it  would  be  almost  zero  in  this  scale  of  Fig.  4.6).  It  is  clear  that  thermal 
cycling  creep  results  in  accumulated  strains  far  in  excess  of  those  from  isothermal  creep  which  is 
consistent  with  numerous  experimental  reports.  Finally,  the  effects  of  progressive  damage  are 
illustrated  in  Fig.  4.7  by  employing  the  statistical  description  assumed  in  the  previous  section  for 
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Figure  4.4  Predicted  creep  strain-time  relation  of  a  SiC  whisker/6061  A1  matrix  composite  for 
various  volume  fraction  of  debonded  fiber  (fd). 


Time  (xlO6  s) 

Figure  4.5  Predicted  creep  strain-time  relation  of  a  SiC  whisker/6061  A1  matrix  composite  for 
various  values  of  shape  parameter  of  Weibull  function  Eq.  (4. 16). 
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Figure  4.6  Predicted  accumulated  strain-thermal  cycles  relation  of  a  SiC  whisker/6061  A1  matrix 
composite  subjected  to  combined  creep  and  thermal  cycling. 


Number  of  Cycles,  N 

Figure  4.7  Predicted  accumulated  strain-thermal  cycles  relation  of  a  SiC  whisker/6061  A1  matrix 
composite  subjected  to  combined  creep  and  thermal  cycling  for  various  values  of  shape  parameter 
b. 
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various  values  of  the  shape  parameter,  b,  at  a  guaranteed  ecr  =  .01.  The  effect  of  b  is  to  broaden 
the  distribution  (as  b  increases)  so  that  the  development  of  damage  after  ecr  is  less  abrupt  as  is 
evident  in  Fig.  4.7!  Although  the  example  shown  in  Fig.  4.7  is  primarily  an  illustration  of  the 
capabilities  of  the  proposed  model,  it  is  evident  that  the  experimentally  observed  phenomena  of 
increasing  nonlinear  strain  accumulation  with  progressive  damage  can  be  described  with  this 
model.  Efforts  to  correlate  the  interfacial  damage  with  observed  strain  accumulation  are  currently 
underway. 
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4.2.  Analytical  Model  for  a  Particulate  Composite 
4.2.1  Introduction 

A  metal  matrix  composite  (MMC)  consists  of  a  matrix  metal  and  reinforcements  and  its 
mechanical  properties  are  normally  superior  to  those  of  the  unreinforced  metal  except  for  the 
fracture  toughness  and  the  thermal  cycling  creep  resistance.  The  main  reason  for  this  defect  of 
MMC  is  the  mismatch  of  the  coefficient  of  thermal  expansion  (CTE)  between  the  matrix  metal  and 
reinforcements  which  is  an  inevitable  characteristic  of  MMCs  [1].  It  is  well  known  that  the 
combination  of  relatively  small  applied  stress  and  cyclic  thermal  loading  can  result  in  increased 
creep  rates  and  superplastic  deformations  in  metal  matrix  composites  [2].  Even  without  an  applied 
stress,  a  detrimental  dimensional  change  has  been  observed  in  a  number  of  thermal  cycled  metal 
matrix  composite  systems  [3,4],  These  effects  are  important  aspects  for  study  in  the  assessment  of 
MMCs  as  candidates  for  high  temperature  structural  applications  requiring  dimensional  stability. 
Although  a  series  of  experiments  on  the  thermal  cycling  creep  behavior  of  discontinuous  MMCs 
have  been  made  [2,5-7],  the  analytical  model  to  predict  dimensional  change  still  remains  to  be 
made. 

For  the  analysis  of  the  enhanced  creep  rate  under  creep/thermal  cycle  loading,  the  determination  of 
internal  stress  is  prerequisite.  The  existence  of  internal  stress  has  been  proved  by  the  use  of  X-ray 
and  neutron  diffraction  methods  and  the  measurements  of  differences  between  compressive  and 
tensile  mechanical  properties  [8].  High  internal  stresses  generated  at  the  interface  between  the  fiber 
and  matrix  are  readily  relaxed  to  a  lower  stress  state  by  a  dislocation  punching.  This  phenomenon 
has  been  observed  in  a  SiC  whisker/Al  composite  in  situ  transmission  electron  microscopy  study 
[9,10], 

Hong  et  al.  studied  the  thermal  cycling  creep  behavior  of  SiC  whisker/2024  A1  composite 
both  analytically  and  experimentally,  and  computed  strain  per  thermal  cycle  by  using  the  matrix 
creep  equation  of  a  hyperbolic  sine  type  [5].  The  average  internal  stress  in  the  matrix  induced  by 
temperature  change  was  estimated  by  using  equation  to  fit  the  thermal  cycling  data  at  low  stresses 
for  the  composite  and  was  considered  to  be  the  local  yield  stress  of  the  matrix.  The  matrix  creep 
equation  was  constructed  from  the  isothermal  creep  data  of  the  composite  in  the  transverse 
direction  instead  of  unreinforced  matrix.  Although  an  applied  stress  is  distributed  to  both  the 
matrix  and  fiber,  they  considered  the  matrix  stress  to  be  the  applied  stress. 

Daehn  et  al.  also  studied  the  thermal  cycling  creep  behavior  of  SiC  whisker/6061  Al 
composite  both  analytically  and  experimentally,  whose  model  predicted  strain  rate  by  using  the 
matrix  creep  equation  of  power  law  type  [6].  The  matrix  stress  consists  of  the  shear  stress  from  an 
applied  stress  and  an  axial  stress  from  a  temperature  change  and  the  resultant  stress  in  the  matrix 
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and  its  direction  were  computed  by  using  2-D  Mohr  circle  approach.  The  elongation  of  the 
composite  is  assumed  to  result  from  the  formation  of  voids  at  the  ends  of  fiber.  Although  the 
model  composite  uses  the  fiber  aspect  ratio  of  10,  the  thermal  stress  in  the  matrix  was  computed  by 
using  the  concept  of  continuous  fiber  reinforced  MMC  and  the  isostrain  assumption.  In 
constructing  the  power  creep  law  for  the  matrix,  the  exponent  of  power  law,  n,  is  selected  as  1 1, 
which  may  be  too  high  for  the  A1  matrix. 

Pickard  and  Derby  extended  the  Greenwood  et  al.  model  to  predict  the  thermal  cycling  creep 
rate  of  SiC  particle/  Al  1100,  where  the  strain  per  thermal  cycle  was  computed  by  using  the 
uniform  time-independent  plastic  deformation  of  the  matrix  [7].  Greenwood  et  al.  model  was 
originally  developed  for  the  study  on  the  phase  transformation  plasticity  of  metals  at  temperatures 
low  compared  with  the  melting  temperature  of  the  materials  and  small  applied  stress  so  that  the 
creep  of  the  matrix  was  neglected  [11],  which  may  not  reflect  the  actual  situation  in  the  MMC 
subjected  to  thermal  cycling  creep  loading  at  higher  temperature  than  the  homologous  temperature. 
When  the  amplitude  of  temperature  cycling  is  not  large  enough,  the  matrix  yielding  is  localized 
around  the  reinforcement  so  that  uniform  yielding  of  the  matrix  is  not  physically  realistic.  The 
critical  temperature  change  required  to  initiate  such  a  full  plastic  flow  on  cooling  and  heating  was 
derived  using  yield  criterion,  but  the  interaction  effect  between  particles  was  not  included. 

Dunn  and  Taya  constructed  a  model  for  the  perfectly-bonded  and  debonded  MMCs  under  the 
condition  of  large  temperature  cycling,  where  strain  per  thermal  cycle  was  computed  by  using  the 
time-independent  plastic  deformation  of  the  matrix  during  cooling  and  heating,  and  the  creep  of  the 
matrix  during  high  temperature  plateau  and  debonding  at  the  fiber-matrix  interface  was  used  to 
explain  the  non-linear  strain  accumulation  of  MMCs  associated  with  the  later  stage  of  thermal 
cycling  [12],  For  the  perfectly-bonded  interface,  their  prediction  of  the  strain  per  thermal  cycle 
becomes  zero  except  the  first  cycle.  If  the  CTE  mismatch  strain  generated  by  temperature  change  is 
small  so  that  the  yielding  of  the  matrix  metal  is  localized  around  fibers,  then  the  assumption  of 
general  yielding  of  the  matrix  in  the  entire  matrix  domain  is  not  valid.  Therefore,  the  internal  stress 
needs  to  be  estimated  based  on  the  localized  plastic  zone. 

The  purpose  of  this  work  is  multi-fold:  to  observe  the  dimensional  change  of  a  MMC  under 
combined  creep  and  thermal  cycling  with  small  temperature  change  where  the  localized  plastic 
deformation  around  particulates  is  expected  to  take  place,  to  develop  an  analytical  model  to  predict 
the  dimensional  change  and  to  compare  the  analytical  predictions  with  the  experimental  results.  In 
this  model,  we  focus  on  the  combined  effects  of  the  internal  stress,  externally  applied  stress  and 
the  relaxation  of  the  stress  induced  by  creep  strain.  The  temperature-induced  internal  stress  is 
calculated  based  on  stress-free  temperature  by  using  dislocation  punching  model  [13]  and  the 
matrix  stress  carried  by  externally  applied  stress  is  also  calculated,  where  the  Eshelby's  equivalent 
inclusion  method  [14]  and  Mori-Tanaka's  mean  field  theory  [15]  are  used.  Finally,  the  predictions 
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will  be  compared  with  the  experimental  results  by  the  use  of  material  properties  of  each  constituent 
and  temperature  profile  as  model  inputs. 


4.2.2.  Analytical  model 

It  is  assumed  in  the  present  model  that  prolate  particulates  (hereafter  particulates)  are  purely 
elastic  and  the  matrix  is  elastic-plastic-creeping.  The  key  analysis  is  the  computation  of  the  average 
stresses  carried  by  the  matrix.  For  the  prediction  of  the  dimensional  change  of  a  MMC  under 
creep/thermal  cycle  loading,  those  stresses  are  induced  due  to  the  three  driving  sources:  (1) 
temperature  change,  (2)  applied  stress  and  (3)  the  creep  and  diffusional  relaxation  of  the  matrix. 
The  first  is  relaxed  by  dislocation  punching,  based  on  which  the  average  matrix  stresses  are 
calculated  at  each  temperature  for  a  temperature  profile  of  thermal  cycle  by  using  temperature- 
dependent  yield  stress.  The  amount  of  the  second  carried  by  the  matrix  is  a  function  of  material 
properties  and  particulate  shape  and  calculated  by  using  Eshelby  type  model.  The  last  is  assumed 
to  be  relaxed  by  a  relaxation  mechanism  such  as  interfacial  diffusion  and  calculated  by  Eshelby 
type  model  as  well.  Finally,  the  strain  per  cycle  is  numerically  calculated  by  using  the  equation  of 
matrix  creep,  the  power  law  of  Dorn  type,  where  the  above  stresses  are  used  as  a  driving  source. 

Matrix  stress  due  to  temperature  change 


Upon  cooling  by  AT  from  a  stress-free  temperature  (T0)  to  a  lower  temperature  (T),  the  CTE 
mismatch  strain  will  be  developed  at  the  particulate-matrix  interface  and  given  by  a*=(af-am)AT 
where  ctf  and  am  are  the  CTEs  of  the  particulate  and  matrix,  respectively  and  AT=T-T0.  These 
strains  are  simulated  by  arrays  of  prismatic  dislocation  loops  attached  to  the  particulate-matrix 
interface  where  a  particulate  is  modeled  by  a  prolate  spheroid  (Qi).  These  loops  are  assumed  to  be 
punched  out  symmetrically  to  reduce  the  otherwise  localized  high  stress  field  at  the  particulate- 
matrix  interface  and  the  punched-out  dislocation  loops  are  represented  by  another  prolate  spheroid 
(Q2)  w^h  e2*  (={1  1  1  00  0}a*/r|)  where  the  aspect  ratio  of  Q2  is  assumed  to  be  the  same  as  Qp 

The  illustration  of  the  relaxed  and  non-relaxed  states  are  given  in  Fig.  4.8. 

The  computational  procedure  used  here  is  the  same  as  that  of  Taya  and  Mori  [16]  except  for  a 
minor  modification,  thus  only  its  brief  summary  is  given  in  Appendix  A.  The  stress  fields  in  the 
donut-shaped  domain  of  (f^'^l)  and  in  the  domain  of  particulate  (f2i)  can  be  expressed  as: 


a2  =C 


m 


n3f)e2>-fel“ 


in  Q2-^l 


(4.29) 
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a2+c‘  =cm  -(S-I)-[(l-Ti3f)e2*+(l  - 


in  £2i 


(4.30) 


where  e1**  is  the  yet  unknown  fictitious  eigensttain  inti'oduced  into  Q\  and  is  solved  as: 

e‘"  =  A“1  (l-'n3f)-(Cm-Cf)-(S-I)-e2* 

where  A  is  given  by  Eq.  (A- 8)  in  Appendix  A. 


(4.31) 


Figure  4.8  Prismatic  dislocation  loops  generated  by  the  mismatch  of  CTEs  between  the  particulate 
and  matrix  and  temperature  change:  (a)  before  punching  and  (b)  after  punching. 


The  normalized  punching  distance,  r|  (=c’/c),  is  determined  by  using  both  the  driving  and 
retarding  forces  of  dislocation  motion.  The  driving  force  is  a  result  of  the  decrease  in  the  total 
potential  energy  as  the  loops  are  punched  out  and  the  retarding  force  is  a  resultant  force  opposing 
both  climbing  and  gliding  motions  of  dislocations.  The  criterion  to  determine  the  punching 
distance  is  proposed  by  Tanaka  et  al.[17]: 


au  aw 

ar|  3r| 


(4.32) 
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where  U,  and  W  denote  the  total  potential  energy  of  the  composite  and  the  total  work  dissipated  by 
the  motion  of  the  dislocation  loops  over  the  punching  distance,  respectively.  The  U  and  W  are 
calculated  from  the  stress  and  strain  fields  in  the  composite  and  the  geometrical  consideration  of 
dislocation  movement  and  expressed  as: 


„3  2  2*  1  2* 

T|  a  e  +g  e 


(4.33) 


W  = 


(P+4) +A/2(l+p-2)( 4(3+1)  f(Tl-l)ka* 


(4.34) 


where  P  and  k  denote  the  particulate  aspect  ratio  and  the  yield  stress  of  the  matrix  in  shear, 
respectively.  The  detailed  derivation  of  U  and  W  is  given  in  Appendix  B. 

By  using  Eq.  (4.30)  and  the  facts  that  the  stress  in  the  particulate  is  proved  to  be  uniform  by 
Eshelby  and  the  integral  of  internal  stress  over  the  composite  domain  vanishes, 
(l-f)om+fof=0,  the  average  stress  in  the  matrix  due  to  the  temperature  change  ,g  ,  is 

computed  as: 


f 

1-f 


gW 


(4.35) 


Matrix  stress  due  to  applied  stress 

The  stress  externally  applied  to  the  composite  material  G°  shown  in  Fig.  4.9  is  shared  by  the 

matrix  and  particulate.  The  stresses  carried  by  each  component  are  determined  by  its  mechanical 
properties  and  particulate  shape.  By  using  the  Eshelby's  equivalent  inclusion  method,  the  average 
stress  carried  by  the  matrix  under  applied  stress  G° ,  g£  ,  can  be  calculated  as[12] 

<  =[l  +  f  C„  •  (S  - 1)  ■  A-1  •  (C,  -  Cm )  ■  c;1]  •  <J°  (4.36) 
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Figure  4.9  Idealized  model  for  SiCp/Al  composite  with  perfectly  aligned  particulates  to  compute 
the  matrix  stress  due  to  an  applied  stress  o° . 

Matrix  stress  due  to  matrix  creep  and  interfacial  diffusion 


As  creep  strain  is  accumulated  in  the  matrix  domain,  the  matrix  is  hardened  by  the  back  stress 
due  to  geometrically  necessary  dislocations  around  particulates  resisting  against  the  otherwise 
homogeneous  plastic  deformation  of  the  matrix.  However,  it  has  been  also  recognized  that  the 
back  stress  state  is  not  stable  and  can  be  relaxed  by  some  mechanisms[18].  It  is  assumed  that 
interfacial  diffusion  takes  place  and  relaxes  some  of  the  buildup  of  back  stress[19].  Since  the 
stress  field  in  the  particulate  due  to  back  stress  and  its  diffusional  relaxation  becomes  hydrostatic,  it 
does  not  affect  further  plastic  deformation  of  the  matrix[20]. 

The  small  amount  of  uniform  creep  strain  in  the  matrix,  dem(=Bdem)>  and  its  diffusional 
strain  in  the  particulate,  ded(=Bded),  induce  the  stresses  in  the  particulate,  doc 
(=-(l-f)U  dem/f)  and  dgd  (=(l-f)U-ded/f),  respectively,  where  B={  -1/2,  -1/2,  1,  0,  0, 
0 }  and  U  =  f  Cm  •  (S  - 1)  •  A-1  •  Cf  •  B .  The  detailed  derivations  of  these  stresses  are  given  in 
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Appendix  C.  From  the  requirement  of  the  diffusional  relaxation,  the  1-  and  3-  components  of 
stresses- satisfy  the  following[20]: 

ciG^  +da^  =  daf  +dof  (4.37) 

By  inserting  Eqs.  (C-7)  and  (C-9)  into  Eq.  (4.37),  the  diffusional  strain  in  the  particulate  becomes 
equal  in  magnitude  to  the  creep  strain  in  the  matrix,  i.e., 

dem=ded  (4.38) 

By  using  the  requirement  of  internal  stress  in  a  composite:  (1  —  f  )o  +f  <7  f=0,  the  resultant  matrix 
stress  dom  generated  by  the  incremental  creep  and  diffusional  strains  is  calculated  and  canceled 
out: 

- 

dom  =  U  •  (dem  -ded)=0  (4.39) 

Strain  per  thermal  cycle 

Since  the  creep  resistance  of  silicon  carbide  is  much  higher  than  that  of  aluminum  due  to  its 
high  modulus  and  low  diffusion  rate,  it  can  be  assumed  that  SiC  particulates  do  not  deform 
plastically  while  A1  matrix  does  during  thermal  cycling/creep.  The  creep  rate  of  a  composite 
depends  only  on  the  matrix  creep  rate  which  is  assumed  to  obey  the  power  law  of  Dom  type[21]. 

From  Eqs.  (4.35),  (4.36)  and  (4.39),  the  matrix  creep  rate  is  obtained  as 

em  =  +  Cm]"  exP(“^>  <4-40) 

where  A0,  T,  R  and  Q  are  Dom  constant,  temperature  in  Kelvin,  gas  constant  and  activation  energy 
respectively  and  matrix  stresses  and  strain  are  expressed  in  terms  of  equivalent  quantities  based  on 
the  Von-Mises  yield  criterion.  Once  the  matrix  creep  rate  is  given,  the  composite  strain  rate,  ec,  is 

calculated  by  adding  the  creep  and  diffusional  strain  rates: 

ec  =  ( 1-f)em  +fed  •  (4-41) 

By  using  Eqs.  (4.38)  and  (4.41),  we  arrive  at 
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(4.42) 


ec  —  em 

Finally,  the  strain  per  thermal  cycle  of  the  composite  can  be  computed  by  numerically  integrating 
Eq.  (4.40)  over  one  thermal  cycle: 

ec=jemdt  (4.43) 

o 


where  x  is  duration  of  one  thermal  cycle. 

4.2.3.  Results  and  discussion 

The  isothermal  creep  data  of  the  matrix  are  required  as  input  data  for  the  analytical  model,  but 
only  limited  creep  data  are  available.  Using  the  creep  data  at  the  temperatures  of  300  °C,  the  power 
law  creep  equation  in  the  form  of  Eq.  (4.40)  is  constructed,  where  Q  is  the  activation  energy  of 
pure  aluminum  for  lattice  diffusion  and  is  obtained  from  the  literature [21],  The  temperature- 
dependent  yield  stress  of  the  matrix  was  measured  and  the  results  are  shown  in  Fig.  4.10.  The 
properties  and  parameters  of  each  constituent  material  used  in  the  present  analysis  are  summarized 
in  Table  4.1[22,23],  where  the  stress-free  temperature  is  assumed  to  be  the  solution-treatment 
temperature  T0  (=530°C). 

Table  4. 1  Material  properties  and  parameters  of  the  particulate  and  matrix  used  as  model  inputs. 


Units 

A1 

SiC 

E 

GPa 

70 

427 

V 

- 

0.33 

0.17 

a 

xlO'6 

23.6 

4.3 

f 

- 

0.8 

0.2 

Aq 

- 

1522 

n 

- 

4.73 

Q 

KJ/mole  • 

142 

T0 

P 

.  °C 

530 

2.25 
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Figure  4.10  The  predictions  of  thermal  residual  stress  in  the  matrix  at  each  temperature  based  on 
the  temperature-dependent  yield  stress  of  the  matrix. 


c 

© 


Is 

■o 

1  C 
£  .2 

§  O 

u 

©  4> 

.c  n 


(« 

B 

0> 

Q 


-90  -45  0  45  90 

Fiber  orientation  angle  (deg.) 


Figure  4.1 1  Results  of  microstructural  studies  on  SiCp/Al  composite  in  terms  of  the  density  of  the 
normal  distribution  of  particulate  orientation:  (a)  as-received  state,  (b)  as-tested  under  the  applied 
stress  of  o°= 24  MPa  and  Tmax=300  °C  and  (c)  Tmax=350  °C,  respectively.  Angle  0°  means  that 
the  particulate  major  axis  coincides  with  the  directions  of  the  extrusion  and  the  applied  stress. 
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Since  SiC  particulates  in  the  composite  are  not  spherical  and  are  misoriented,  the  distributions 
of  their  aspect  ratio  and  orientation  angle  are  needed  to  predict  dimensional  change  per  cycle.  The 
results  of  the  aspect  ratio  of  SiC  particulates  and  the  particulate  orientation  angle  distributions  of  the 
as-received  composite  are  shown  in  Figs.  2. 1 1  (a)  and  4. 1 1  (a),  respectively.  These  figures  reveal 
that  particulates  are  mostly  oriented  in  the  direction  of  extrusion  and  their  average  aspect  ratio  is 
2.15,  where  angle  0°  denotes  the  directions  of  the  extrusion  and  applied  stress  loading  axes.  For 
the  simplification  of  the  analysis,  we  assume  that  particulates  are  perfectly  aligned  in  the  direction 
of  extrusion.  Fig.  2.11  also  confirms  that  the  average  aspect  ratio  hardly  changes  during  test. 
Takao  et  al.  have  examined  the  effect  of  variable  aspect  ratios  on  the  stiffness  and  thermal 
expansion  coefficients  of  a  short  fiber  composite  and  concluded  that  the  results  predicted  by  using 
the  average  value  of  aspect  ratio  are  not  much  different  those  of  using  variable  aspect  ratio  unless 
the  distribution  of  aspect  ratio  is  extraordinarily  scattered[24].  Thus,  the  use  of  the  average  aspect 
ratio  of  2.25  in  the  model  calculations  is  justified.  Furthermore,  the  microstructural  change  due  to 
the  fracture  of  particulates  during  testing  was  examined  by  using  micrographs  shown  in  Fig.  2.10. 
The  micrograph  of  Fig.  2.10  (c)  showed  the  higher  possibility  of  the  fractured  particulates  than 
Figs.  2.10  (a)  and  2.10  (b),  while  the  micrograph  of  Fig.  2.10  (b)  showed  the  less  possibility  of 
the  fractured  particulates  than  Fig.  2.10  (a).  It  can  be  concluded  from  these  observations  that  the 
fracture  of  particulates  does  not  take  place  during  the  creep  thermal  cycling,  but  during  the 
processing  of  the  composite. 

The  predicted  internal  stresses  in  the  matrix  due  to  the  temperature  change  are  plotted  as  a 
function  of  temperature  in  Fig.  4.10,  which  are  expressed  as  the  equivalent  stresses  (ce),  <*3-01, 
based  on  Von-Mises  yield  criterion.  Comparing  them  with  the  yield  stress  of  the  matrix  shown  in 
Fig.  4.10,  it  can  be  concluded  that  the  temperature-induced  stress  strongly  depends  on  the  yield 
stress  of  the  matrix. 

The  strains  per  thermal  cycle  at  various  applied  stresses  observed  from  the  SiCp/Al 
composites  that  were  thermal-cycled  with  two  different  temperature  profiles  of  Tmax=350  °C  and 
300  °C  are  plotted  in  Fig.  4.12.  The  experimental  results  show  that  the  strain  per  thermal  cycle  is 
an  increasing  function  of  an  applied  stress  for  both  temperature  profiles.  Using  the  input 
parameters  tabulated  in  Table  4.1,  the  predictions  of  the  strain  per  thermal  cycle  of  the  composite 
under  combined  creep/thermal  cycling  by  the  present  model  are  also  plotted  in  Fig.  4.12. 

For  the  creep  thermal  cycle  of  Tmax=350  °C,  the  predicted  results  agree  well  with  the 
experimental  results.  Although  there  exists  a  reasonably  well  agreement  between  the  predictions 
and  experiment,  the  predictions  of  the  strain  per  cycle  for  the  creep  thermal  cycle  of  Tmax=300  °C 
overestimate  the  experimental  results.  This  discrepancy  can  be  explained  as  follows.  First,  it  is 
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Figure  4.12  Comparison  between  the  analytical  and  experimental  results  of  strain  per  thermal  cycle 
of  SiCp/Al  composite  subjected  to  temperature  cycling  and  various  applied  stresses  o°:  (1)  the  solid 
and  open  circles  represent  the  experimental  data  for  Tmax=350  °C  and  300  °C,  respectively  and  (2) 
the  solid  and  dotted  lines  denote  the  analytical  results. 

assumed  that  the  interfacial  diffusion  takes  place  almost  instantaneously  at  each  time  during  thermal 
cycling.  According  to  the  models  developed  by  Mori  and  his  co-workers,  the  characteristic  time  of 
the  interfacial  diffusion  is  strongly  dependent  on  the  temperature  and  the  calculation  of  the 
characteristic  times  for  T=350,  300,  250  and  200  °C  gives  0.5,  1.5,  5  and  20  seconds, 
respectively [20, 25].  Thus,  the  recovery  of  back  stress  by  diffusional  relaxation  is  limited  at  low 
temperatures  so  that  the  unrelaxed  back  stress  in  the  matrix  lowers  the  extent  of  the  plastic 
deformation  of  the  matrix.  Secondly,  it  has  been  also  known  that  although  SiC  particulates  are 
randomly  distributed  in  the  matrix  before  test,  they  are  reoriented  in  the  direction  of  an  applied 
stress  during  testing[2,26].  This  phenomenon  was  also  observed  by  the  microstructural  studies  on 
the  as-received  state  and  as-tested  states,  whose  results  are  depicted  in  terms  of  the  density  of  the 
normal  distribution  of  particulate  orientation  and  summarized  in  Fig.  4.1 1.  The  angle  of  particulate 
orientation  is  defined  as  the  angle  between  the  particulate  major  axis  and  the  applied  stress  or 
extrusion  direction.  The  composite  with  perfectly  aligned  particulates  gives  rise  to  higher 
dimensional  change  (strain  per  cycle)  as  compared  with  the  composite  with  misaligned  particulates. 
Because  of  less  particulate  reorientation  at  Tmax=300oC,  the  degree  of  misalignment  of  SiC 
particulates  in  a  composite  tested  at  Tmax=300°C  is  larger  than  that  at  Tmax=350°C.  Therefore,  the 
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Volume  fraction  of  fiber,  f 


Figure  4.13  The  predicted  results  of  the  effect  of  particulate  volume  fraction  f  on  strain  per  thermal 
cycle  of  SiCp/Al  composite  subjected  to  the  temperature  cycling  of  Tmax=350  °C  and  the  applied 
stress  of  o°=24  MPa. 


Fiber  aspect  ratio,  (3 


Figure  4. 14  The  predicted  results  of  the  effect  of  particulate  aspect  ratio  (3  on  strain  per  thermal 
cycle  of  SiC/Al  composite  subjected  to  the  temperature  cycling  of  Tmax=350  °C  and  the  applied 
stress  of  c°=24  MPa  by  the  present  model  and  macroscopic  model^l. 
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predictions  of  strain  per  thermal  cycle  for  the  lower  thermal  cycle  of  Tmax=300  °C  become  higher 
than  the  experimental  values. 

A  parametric  -study  is  conducted  to  examine  the  effects  of  particulate  volume  fraction  f  and 
particulate  aspect  ratio  for  the  case  of  o°= 24  MPa  and  Tmax=350°C.  The  results  of  strain  per 
thermal  cycle  as  a  function  of  f  are  shown  in  Fig.  4.13,  where  the  maximum  strain  per  cycle  is 
obtained  around  f=0.2.  This  agrees  with  the  observations  by  Pickard  et  al.  where  the  peak  value 
of  strain  per  cycle  takes  place  between  f=0.2  and  0.3[7].  The  effect  of  fiber  aspect  ratio  on  the 
strain  per  thermal  cycle  of  the  composite  was  examined  by  using  two  models,  the  present  model 
and  macroscopic  model.  The  macroscopic  model  is  the  extension  of  the  model  for  the  pure  thermal 
cycling  of  MMCs  made  by  Taya  and  Mori[27]  and  assumes  the  idealized  temperature  profile  of 
step-function  type,  a  plastic  deformation  of  the  matrix  during  cooling  and  heating  and  interfacial 
diffusion  during  high  temperature  plateau  process.  The  strain  per  thermal  cycle  by  the  present 
model  results  in  much  higher  strain  per  thermal  cycle  than  the  macroscopic  model,  as  fiber  aspect 
ratio  increases.  This  can  be  explained  by  the  fact  that  the  mode  of  dislocation  punching  considered 
in  the  present  model  was  the  simultaneous  punching  along  both  longitudinal  and  transverse 
directions,  which  is  considered  to  be  valid  only  for  small  particulate  aspect  ratio.  Since  the 
dislocation  punching  has  preferred  direction  along  the  fiber  major  axis  as  fiber  aspect  ratio 
increases^18],  the  equivalent  stress  (cre)  of  the  matrix  induced  by  temperature  change  is 
overestimated  by  the  present  model,  and  the  matrix  creep  rate  is  strongly  dependent  on  ae.  While 
the  macroscopic  model  gives  strain  per  thermal  cycle  only  above  the  critical  fiber  aspect  ratio  (see 
Fig.  4.14).  The  assumption  of  uniform  yielding  of  the  matrix  during  cooling  is  more  compatible 
with  a  MMC  with  larger  aspect  ratio  of  fibers.  Therefore,  the  validity  of  the  present  model  is 
restricted  to  small  aspect  ratios  up  to  3,  beyond  which  the  other  models  such  as  our  previous 
model[27]  should  be  applied.  For  the  range  between  3  and  the  critical  aspect  ratio,  the  present  and 
macroscopic  models  can  not  predict  the  dimensional  change  accurately.  It  should  be  noted  here 
that  as  f  increases,  this  discrepancy  of  the  predicted  dimensional  change  between  the  present  and 
the  macroscopic  models  decreases,  resulting  in  more  smooth  transition. 
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Appendix  A:  stress  fields  in  Qi  and  ^2-^1 

The  composite  domain  is  divided  into  three  regions:  (1)  the  matrix  domain  without 
dislocations  (D-Q2X  (2)  the  punched-out  domain  (Q2X  and  (3)  fiber  domain  (Qi),  Fig.  4.8.  With 
equation  (3),  the  problem  is  changed  to  that  of  an  inhomogeneous  inclusion  Qj  with  Cf  and  e2* 
surrounded  by  a  matrix  with  Cm  and  e2*  which  in  turn  is  embedded  in  an  infinite  matrix  domain. 
First  apply  Eshelby's  equivalent  inclusion  method  to  a  hypothetical  composite  which  contains  the 
punched-out  domains  (0-2)  with  uniformly  distributed  eigenstrain  e2*  (i.e.,  fibers  are  replaced  by  a 
material  having  matrix  elastic  constants).  Then  the  resulting  internal  stress,  o2 ,  in  Q.2  is  given  by 

o2=Cm-(e  +  e2-e2*)  (A-l) 

where 

e2  =  S2  •  e2*  (A-2) 
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where  e2  and  S2  are  the  internal  strain  caused  by  inclusion  Oi  and  Eshelby's  tensor  related  to  O2 
and  superscript  2  denotes  the  quantity  related  to  Cb- 

Next  introduce  fiber  (Qj)  into  the  inside  of  ^2-  The  stress  field  in  the  fiber  domain  caused 
by  this  introduction,  a2  +  a1 ,  is  given  by  Eshelby's  equivalent  inclusion  method: 


a2  +G1  =Cf  •(e  +  e2-e2*+e1 ) 

=  C  •(e  +  e2-e2H'+e1-e1**) 

mV  / 


(A-3) 


where  Cf  and  Cm  are  the  stiffness  tensor  of  the  short  fiber  and  the  matrix,  respectively.  Since  e1 
is  the  disturbance  strain  due  to  the  introduction  of  into  O2,  it  is  related  to  the  fictitious 
eigenstrain  (e1**)  in 


e1  =  S1  •  e1** 


(A-4) 


where  S1  is  Eshelby's  tensor  of  Qi. 

Since  o1  and  G2  are  internally  induced  in  the  composite,  the  volume  average  of  the  internal 
stress  over  the  entire  composite  domain  must  vanish: 

j(o1+a2)dV=0  (A-5) 

From  this  requirement  with  Eqs.  (A- 1~4),  e  can  be  expressed  as: 

e  =  -f2(s2-l).e2*-f(s1-l)e1".  (A-6) 


Here,  f  and  f2  denote  the  volume  fractions  of  fibers  and  the  punched-out  domains,  respectively  and 
are  related  as  f2=ri3f.  S2  denote  Eshelby's  tensor  of  0.2  and  is  equal  to  S=Si.  From  Eqs.  (A- 

1~6),  the  unknown  e1**  can  be  solved  as 


e‘"  =  A-'.(l-f2)(Cm-Cf).(S2-l)-e2* 

(A-7) 

where 

A  =  {(C,-C„)-[(l-f)S1+fl]+C„} 

(A-8) 
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By  inserting  Eqs.  (A-2),  (A-4)  and  (A-7)  into  Eqs,  (A-l)  and  (A-3),  the  stresses  in  each 
domain  result  in  the  followings: 


g2  =C  •  ( S-I) ■  I  ( 1— r|3 f )e2* -f  e1**  in  Q2-^i 


g2+g‘  =  Cm  ■  ( S-I) •  [ ( l-r|3 f)e2*  +( 1-  fje1**]  in  Qi 


(A- 10) 


nx  n:  work  aissn 


In  order  to  estimate  the  dissipated  energy  during  the  motion  of  dislocation  loops,  we 
decompose  the  dislocation  loops  into  two  groups.  The  first,  designated  by  L,  is  to  relax  the  misfit 
strain  in  longitudinal  direction;  and  the  other,  designated  by  T,  is  to  relax  the  misfit  strain  in 
transverse  direction.  The  dissipated  energy  due  to  the  dislocation  punching  of  the  L-type 
dislocation  loops  will  be  derived  in  detail. 


Figure  B:  A  schematic  figure  depicts  that  a  dislocation  loop  Di  is  punched  out  to  a  dislocation  loop 
Do  mathematically. 

The  dislocation  punching  considered  in  Figure  B  is  based  on  the  geometrical  observation  that  the 
dislocation  located  at  Di  is  punched  out  to  D2.  In  reality,  the  dislocation  punching  takes  place  by 
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the  combination  of  gliding  and  climbing  motions,  among  which  dislocation  climb  usually  occurs  at 
high  temperatures,  typically  above  half  the  melting  temperature  of  the  matrix[28].  Since  the 
frictional  stresses  -for  gliding  and  climbing  motions  of  dislocations  are  defined  by  k  and  2k 
respectively,  the  resultant  friction  stress  along  the  0  direction  for  L-type  dislocation  loops,  G^,  can 

be  expressed  as 


Gq  =  k(  2cos0+sin0) 


(B-l) 


where  k  is  the  yield 


stress  of  the  matrix  material  in  shear  deformation  and  tan  0  = 


IH 

a/W 


The 


number  of  dislocation  loops  generated  at  the  fiber-matrix  interface,  Nl,  can  be  obtained  using  the 
relationship  proposed  by  Ashby[29],  NAb  =  AV,  and  is  given  by  NL=2ca*/b,  where  N,  A,  AV 
and  b  denote  the  number  and  area  of  dislocation  loops,  the  volume  change  due  to  the  CTE 
mismatch  and  the  magnitude  of  the  Burgers  vector  of  the  matrix,  respectively.  In  addition,  the 
distance,  Dl,  over  which  the  dislocation  located  at  Di  moves  to  D2  and  the  length  of  the  prismatic 
dislocation  loop,  Ll,  are  expressed  as: 


DL=(n-l)a^2+P2M2) 


(B-2) 


LL=27taVl-^2  (B-3) 

where  £,  =  X3/C  and  (3  denotes  the  aspect  ratio  of  fiber.  Once  the  number  of  dislocation  loops,  the 
resultant  friction  stress,  the  distance  moved  and  the  length  of  dislocation  are  known,  the  dissipated 
work  per  the  unit  volume  of  the  composite  for  the  L-type  dislocation  loops,  Wl,  can  be  expressed 
as 


WL  =f((3+4)(ri-l)ka*  (B-4) 

Using  the  same  method  used  in  the  calculation  of  the  dissipated  work  due  to  the  punching  process 
along  the  longitudinal  direction,  the  dissipated  work  per  unit  volume  of  the  composite,  Wj,  due  to 
the  transverse  punching  process  can  be  calculated  as: 

WT  =  V2fVl+(3-2(4p+l)(Tl-l)ka*  (B-5) 
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Therefore,  the  total  dissipated  work  due  to  both  longitudinal  and  transverse  punching  processes  is 
the  sum  of  Eqs.  (B-4)  and  (B-5)  and  is  equal  to: 


W  = 


(P+4)+V2(l+p-2)(4p+l)  f(T]-l)ko* 


(B-6) 


Appendix  C:  fiber  stress  due  to  the  creep  and  diffusional  strains 

The  small  amount  of  uniform  creep  strain  in  the  matrix,  dem ,  and  its  diffusional  strain  in  the 
fiber,  ded ,  can  be  represented  by  using  the  assumption  of  incompressibility  and  axis-symmetry  of 
deformation  mode: 

dem  =  Bdem  (C-l) 


ded  =Bded 


(C-2) 


where  B={  -1/2,  -1/2,  1,  0,  0,  0}.  By  applying  the  Eshelby’s  equivalent  inclusion  method  to  a 
composite  whose  matrix  has  dem,  the  stress  in  the  fiber,  doc ,  due  to  the  creep  strain,  dem,can  be 

expressed  as: 


dac  —  Cf  •  (e +e+dem ) 
=Cm-(e+e-e*) 


(C-3) 


where  the  disturbed  strain  e  is  related  to  the  unknown  eigenstrain  e*  defined  only  in  the  fiber  by 

e  =  S  ■  e*  (C-4) 

Requirement  of  the  integration  of  the  internal  stress  over  the  entire  composite  domain  being  zero 
leads  to 


e+f (e-e*  )=0 

Eliminating  e  and  e  from  Eq.  (C-3)  by  using  Eqs.  (C-4)  and  (C-5),  e*  is  given  by 

e*  =  ~  A-1  •  Cf  •  dem 


(C-5) 


(C-6) 
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where  A  is  given  by  Eq.  (A-8). 

By  inserting  Eqs.  (C-4-6)  into  Eq.  (C-3),  the  stress  in  the  fiber  due  to  the  creep  strain  can  be 
expressed  as 


dem 


(C-7) 


where 


U  =  f  Cm  •  (S  - 1)  •  A-1  •  Cf  •  B 


(C-8) 


By  replacing  dem  with  -ded  in  Eq.  (C-7),  the  stress  in  the  fiber  due  to  the  diffusional  strain,  ded, 
can  be  obtained  as: 


5.  Analytical  Study  III:Stress  Relaxation  by  Plastic  Flow,  Interfacial  Sliding, 

and  Diffusion  in  an  Inclusion  Bearing  Material 


5.1  Introduction 

.  Many  problems  in  materials  science  and  mechanics  are  concerned  with  the  relaxation  of 
stresses,  externally  imposed  or  internally  produced.  Relaxation  is  a  phenomenon  of  decreasing 
total  energies  by  the  occurrence  of  inelastic  deformation,  and  the  state  of  a  minimum  energy  is 
called  the  complete  relaxation.  This  subject  was  studied  for  several  cases  in  terms  of  elastic  states 
and  kinetics  in  the  past,  interfacial  or  boundary  sliding  [1,2],  interfacial  diffusion  on  th.e  interface 
of  or  bulk  diffusion  around  an  inclusion  (second  phase  particle)  [3-8],  and  local  plastic  flow  such 
as  dislocation  punching  [9].  Interfacial  diffusion  was  also  used  to  study  creep  of  a  composite  [10, 
11]. 

The  present  paper  discusses  general  conditions  of  elastic  states  after  the  complete  relaxation 
caused  by  plastic  flow,  diffusion,  and  interfacial  sliding  or  the  combination  of  the  latter  two.  The 
analysis  is  based  on  variational  methods,  with  which  the  elastic  state  after  the  complete  relaxation 
can  be  clearly  perceived.  Each  mode  of  the  above  mentioned  deformation  of  plastic  character 
requires  a  certain  physical  condition.  It  is  a  subsidiary  condition  imposed  on  finding  the  minimum 
state  of  energy.  Examining  subsidiary  conditions,  one  can  see  whether  a  mode  is  more  effective  in 
the  relaxation  than  another.  Some  examples  will  be  given  to  a  specific  mode  of  the  relaxation. 

5.2.  Elastic  State  After  Complete  Relaxation 


Let  Gy  be  the  stress  present  before  the  occurrence  of  stress  relaxation.  For  simplicity,  we 
consider  that  Gy  is  caused  by  an  external  traction  on  the  boundary  of  a  body  D.  However,  Gy  can 

be  an  internal  stress. 

5.2.1  Plastic  flow 

Suppose  that  plastic  strain  £y 

displacement  and  stress  caused  by 
* 

£y  is  written  as 


is  introduced  in  a  domain  Q  to  relax  Gy\  Let  uj  and  Gy  be  the 
£y,  respectively.  The  total  potential  energy  change,  F,  due  to 
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“  ~  |^ij(ui,j  eij)dv  j^ij  ^ijdv, 


(5.1) 


D 


Q 


after  applying  Colonetti's  theorem.  The  first  term  in  Eq.  (5.1)  is  the  elastic  energy  induced  by  Ey 
and  the  second  term  is  the  change  in  the  potential  energy  of  the  external  traction  induced  by  £y.  If 
Gy  is  an  internal  stress,  the  second  term  is  the  interaction  elastic  energy  between  Gy  and  £y. 
Since -ey  is  of  plastic  character,  we  have  the  condition 

£”=0  inQ.  (5.2) 


The  minimum  of  Eq.  (5.1)  under  the  condition  Eq.  (5.2)  is  found  by  minimizing 

=A  Jcy(Ujj  -Ey)dv-  Joy  £ydv  +  j  A,(x)ejdv  (5.3) 

2d  a  a 

where  X(x)  is  the  Lagrange  multiplier.  The  variation  of  Eq.  (5.3)  with  respect  to  £y  is  calculated 
as 


5G 


J(Gy  +  Gy  -X5y)8£ydv 
£2 


(5.4) 


Here,  5y  is  the  Kronecker  delta.  Since  §£y  is  now  arbitrary  after  X  has  been  introduced,  G  or  F  is 
minimized  or  equivalently  the  complete  relaxation  occurs  when  the  total  stress,  Gy  +  Gy  ,  becomes 


Gy  +  Gy  =  X8y  in  Q. 

Further,  the  equilibrium  condition  leads  to 

X  =  const. 


(5.5) 


(5.6) 


That  is,  the  complete  relaxation  is  attained  when  the  total  stress  in  a  domain  undergoing  plastic 
deformation  becomes  uniformly  hydrostatic,  as  schematically  shown  in  Fig.  5.1  (a).  This  holds 
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If  Q  is  an  ensemble  of  unconnected  subdomains,  Qi,  0.2, ...,  the  uniform  hydrostatic 
stress  changes,  in  general,  from  subdomain  to  subdomain.  In  particular,  if  a  part  of  the  surface  of  a 
subdomain  is  on  a  normal  traction  free  surface  of  the  body  D,  the  stress  in  the  subdomain  vanishes 
after  the  complete  relaxation. 

We  can  exchange  the  role  of  inclusions  with  that  of  their  surrounding  matrix.  If  only  the 
matrix  deforms  plastically,  the  equilibrium  state  is  achieved  when  the  total  stresses  in  the  maxtrix 
become  uniform  and  hydrostatic.  As  will  be  discussed  later,  a  practical  subject  covered  by  this 
case  is  creep  in  a  composite  or  dispersion  hardened  alloy.  Reinforcements  in  these  materials  are 
plastically  strong. 

Here,  we  will  concentrate  our  discussion  on  a  metal  matrix  composite.  A  composite  under  a 
constant  external  stress  stops  creep  eventually,  as  the  matrix  approaches  a  uniform  hydrostatic  state 
with  an  increase  in  plastic  (creep)  strain  in  it.  (Of  course,  this  holds  as  long  as  other  relaxation 
mechanisms,  such  as  those  to  be  discussed  later,  do  not  operate.)  Exactly  speaking,  a  creep  rate 
approaches  zero  after  a  sufficient  time  of  creep.  Using  the  average  field  theory  [12],  this 
conclusion  has  been  shown  to  be  correct  [13]. 

Also,  we  will  show  a  result  of  two-dimensional  FEM  analysis  to  further  substantiate  this 
conclusion.  An  element,  as  shown  in  Fig.  5.2,  contains  a  quarter  of  a  cylindrical  inclusion 
surrounded  by  a  creeping  matrix.  Young's  modulus  of  the  cylinder  is  2  x  10^  MPa  and  its 
Poisson  ratio  0.33.  The  matrix  creep  rate  was  assumed  to  be  proportional  to  the  cube  of  a  tensile 
stress  under  a  uniaxial  loading  condition.  A  constant  load,  corresponding  to  the  average  external 
tensile  stress  of  50  MPa,  was  applied  along  the  vertical  direction  (X3-axis).  The  Xi-axis  was 
parallel  to  the  horizontal  direction.  The  origin  of  the  coordinate  system  was  in  the  center  of  the 
inclusion.  The  following  boundary  conditions  were  used:  ui=0  on  Xi=0  and  U3=0  on  X3=0.  For 
the  other  boundaries,  the  periodic  boundary  condition  was  imposed,  by  tying  the  nodes,  i.e., 
keeping  the  total  resultant  force  on  the  Xy -surface  (boundary  perpendicular  to  the  Xi-axis)  zero 
and  the  displacement  some  constant  value.  The  actual  FEM  calculation  was  earned  out,  using  the 
software,  MARC,  Version  6.  We  found  that  as  time  passed,  the  creep  rate  monotonically 
decreased.  We  believe  that  it  must  have  approached  zero,  if  a  sufficiently  long  computing  time  had 
been  allowed.  At  this  stage,  the  von  Mises  stress  became  almost  uniform  and  negligibly  small  in 
the  creeping  matrix,  as  shown  in  Fig.  5.2.  Because  of  the  creep  law,  proportional  to  the  cube  of 
the  stress,  further  creep  deformation  became  extremely  slow  in  the  stage  corresponding  to  Fig. 
5.2.  The  details  of  this  study  will  be  reported  in  a  separate  paper. 

Contrary  to  the  above  conclusion,  the  FEM  studies  of  Dragon  and  Nix  [11]  and  Kim  and 
McMeeking  [14]  on  creep  of  a  composite  ascertained  the  existence  of  a  stationary  state.  This  is 
very  puzzling,  in  the  light  of  the  above  theoretical  analysis,  which  has  also  been  numerically 
confirmed  by  our  own  FEM  analysis.  This  is  not  due  to  the  difference  in  the  creep  law,  as  can  be 
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Figure  5.2.  Two-dimensional  FEM  analysis.  An  elastic  cylindrical  inclusion  is  surrounded  by  a 
creeping  matrix.  The  figure  shows  the  map  of  the  von  Mises  stress  distribution  when  a  creep  rate 
becomes  negligibly  small  after  a  sufficient  creep  time. 
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seen  in  Ref  [14],  but  probably  due  to  different  boundary  conditions.  Kim  and  McMeeking  used 
the  boundary  conditions:  U3=0,  C]3=0  on  X3=0,  ui=0,  G]3=0  on  X3=0  and  G|3=0  on  the  other 
boundaries.  Dragon  and  Nix  employed  the  conditions  ui=0,  G|2=0  on  Xi=0,  U2=0,  Gi2=0  on 
X2=0,  U2,i=0,  G|2— 0  on  the  other  boundary  normal  to  the  fiber  axis  and  up2=0,  G|2=0  on  the 
boundary  parallel  to  the  fiber  axis.  The  adequacy  of  these  boundary  conditions  should  be 
reexamined. 

Approximate  analytical  modeling  of  creep  of  a  short  fiber  composite  also  led  to  stationary 
creep  [15].  The  following  discussion  shows  that  this  result  is  incorrect.  A  composite  or 
dispersion  strengthened  alloy  shows  larger  work  hardening  than  a  single  phase  material.  Analysis 

based  on  continuum  mechanics  or  dislocation  theory  explains  the  larger  work  hardening  [16,17]. 
At  the  same  time,  a  hardening  rate  is  shown  to  be  independent  of  yield  strength  Gy  as  long  as  the 

amount  of  plastic  strain  is  small.  This  is  valid  in  the  extreme  case  of  Gy  — >0.  Thus,,  we  have  a 

one-to-one  correspondence  between  the  external  stress  and  plastic  strain.  If  plastic  deformation  is 
time  dependent,  the  equilibrium  plastic  strain  for  a  given  external  stress  will  be  reached  eventually 
after  a  long  time.  This  means  that  creep  of  a  composite  tends  to  stop,  without  attaining  a  stationary 
state. 

Example  1:  dislocation  punching  from  a  spherical  inclusion 


An  infinitely  extended,  elastically  isotropic  and  uniform  body  contains  a  spherical  inclusion 
Q0  with  the  radius  a,  Fig.  5.3.  g,j  is  thermal  stress  caused  by  the  mismatch  of  thermal  expansion 

between  the  inclusion  and  surrounding  matrix.  The  pure  shear  state  of  G^  outside  the  inclusion 

can  be  relaxed  by  dislocation  punching  or  equivalently  by  plastic  deformation.  In  fact,  it  has  been 
shown  that  if  the  plastic  strain  characterized  by  the  strain  [9] 


* 


3  T 
=  a  £ 


§ii 

-JL-3 

3  ^ 


xiXj 


(5.7) 


occurs  in  Q,  the  total  stress  in  Q.  becomes  hydrostatic;  that  is 


+  aA=_M±^eTW38ij  inn. 

1  lJ  3(1 -v)  U )  1J 


(5.8) 


Here,  Q  stands  for  the  spherically  symmetric  domain  surrounding  Qq  with  the  radius  R.  The  center 
of  Q  coincides  with  that  of  Q0.  p  is  the  shear  modulus  and  v  is  the  Poisson  ratio.  eT  is  the 
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difference  of  thermal  expansion  between  the  inclusion  and  matrix,  causing  the  eigenstrain 
ejj  =  eT8y  in  Q.0.  Although  the  hydrostatic  stress  in  Q  is  constant  as  required  by  Eq.  (5.6),  Gy 

and  similarly  Gy  depends  on  the  position  x  in  Q.  Gy  is  written  as 

a_2|!(1  +  v)  i  t 

O  —  c  a 

1J  3(1 -v) 

The  reason  why  £y  is  written  in  the  form  of  Eq.  (5.7)  is  given  in  Appendix  1. 


°ii 

_y_  t. 

3  -*  5 

XJ  YJ 


XiXj 


in  Q. 


(5.9) 


Figure  5.3.  Inifinite  matrix  contains  a  spherical  inclusion  (Q0)  surrounded  by  concentric  spherical 
domain  (Q)  where  discloations  are  punched  from  the  boundary  of  Q0- 

5.2.2  Sliding  on  the  interface  of  an  inclusion 

This  subject  was  extensively  studied  by  Mura  and  his  co-workers  [18,  19],  assuming  that  the 
total  tangential  traction  on  the  interface  |  Q I  of  an  inclusion  vanishes  after  sliding  on  |  Q I .  The 
following  analysis  justifies  this  assumption  and  indeed  shows  that  it  is  the  energy  minimum 
condition.  The  pertinent  quantities  appearing  in  the  following  discussions  are  defined  in  Fig.  5.4. 
Let  the  sliding  on  |  Q.  \  be  bp  b;  is  the  jump  in  displacement  uj  at  I Q I ,  that  is 

bi  =  ui  (in)  -  uj  (out)  on  |q|.  (5.10) 
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bi  defined  by  Eq.  (5.10)  is  generally  called  Somigliana  dislocations. 

Ini 


Figure  5.4.  An  inclusion  Q,  its  surface  |  Q  | ,  a  unit  normal  vector  ^  and  the  jump  in  displacement 
b  are  defined  in  the  figure. 

The  total  energy  F[b]  is  written  as 

F[b]  =  —  |  cr^njbids  +  jcjjnjbjds  (5.11) 

2  Id  lol 

where  Ojj  is  the  stress  due  to  bj  and  nj  is  the  unit  outward  normal  vector  on  |  Q  | .  The  two  terms  in 
the  right  hand  side  of  Eq.  (5.1 1)  have  identical  meanings  to  those  in  Eq.  (5.1).  bi  must  satisfy 

-njbi  =  0  on  |  £2  | .  (5.12) 

This  equation  means  that  the  Somigliana  dislocations  under  consideration  are  of  glide  type. 

Under  Condition  Eq.  (5. 12),  Eq.  (5. 1 1)  becomes  mimimum  when 

G[b]  =  - |<7ijnjbids+ joy  njbjds- j^njbids  (5.13) 

2  Id  Id  Id 
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reaches  a  minimum  state.  Taking  the  variation  of  G  with  respect  to  bj,  the  condition  for  the 
minimum  of  G  is  found  as 

Gynj  +  Gy nj  =  X(x)nl  on  |q|-.  (5.14) 

That  is,  the  complete  relaxation  by  interfacial  sliding  occurs  when  the  total  traction  on  the  interface 
becomes  normal  to  the  interface  or  the  tangential  traction  vanishes.  This  result  is  shown  in  Fig.  5.1 
(b). 

Example  2:  a  spherical  inclusion  in  an  isotropic  body 

This  subject  was  studied  by  Ghahremani  [20]  and  Mura  and  Furuhashi  [18].  Mura  and 
Furuhashi  took  an  approach  similar  to  the  present  one  but  gave  an  incorrect  factor  in  the  expression 
for  bj.  This  was  corrected  later  [21]  and  is  recorded  in  Appendix  2  in  the  present  paper. 

Example  3:  a  circular  cylinder  in  an  isotropic  body 

A  circular  cylinder  Q.  of  the  radius  a  exists  parallel  to  the  X3  axis.  An  external  stress  Gy  (1,  j 
=  1,  2)  is  applied  at  infinity.  The  tangential  traction  due  to  Gy  is 

TjA  ='Gy  nj  — 'G^nkn^ni  on  |  Q  \ .  (5.15) 

where  'Gy  is  the  deviatoric  part  of  Gy  .  Since  the  hydrostatic  stress,  (Gfj  +  G22)/2,  does  not 

cause  the  tangential  traction,  it  can  be  excluded  in  analysis.  Noting  the  form  of  Eq.  (5.15),  we 
postulate  that  bi  is  given  as 


bj  =Bya2Xj-B^xkx^xi.  (5.16) 

Using  the  paper  by  Asaro  [22]  and  Asaro  and  Barnett  [23],  the  stress  Gy  due  to  Eq.  (5.16)  is 
calculated  as 

aij=_2(^tij5,jBk'XkX'“2a^oB|j(a2'3x2)  in£2'  <517) 

Therefore,  if  By  is  taken  as 


92 


(5.18) 


Cl- v) 
'  M-a2 


the  total  tangential  traction  |  £2  |  vanishes  as  verified  by 


aijnj  -^Wnkn^ni  +  '^ij  nj  -  'a^nkn^ni  =0. 


(5.19) 


5.2.3  Local  diffusion  around  an  inclusion 

What  we  are  concerned  with  here  is  the  relaxation  achieved  by  the  removal  of  some  mass 
from  certain  areas  and  the  deposit  of  them  on  other  areas  (mass  transport).  At  the  moment,  we  do 
not  specify  the  path  of  diffusion.  This  subject  was,  to  our  knowledge,  first  studied  by  Koeller  and 
Raj  [3]  for  a  cylindrical  inclusion,  confining  the  diffusion  path  on  the  interface.  Their  solution  is 
incomplete  in  the  sense  of  mechanics.  Thus,  we  will  show  our  analysis  here. 

Eq.  (5. 1 1)  is  again  used.  Instead  of  Eq.  (5.12),  the  conditions  described  by 

bj  =a(x)nj  on  |Q|  (5.20) 

and 

-  Jnjbjds  =  -  jads  =  0  (5.21) 

Ini  Ini 

are  imposed.  The  removal  or  deposit  of  some  mass  on  |  £2  |  leads  to  Eq.  (5.20)  and  the  local 
conservation  of  mass  to  Eq.  (5.21).  Although  the  diffusion  path  is  not  specified,  the  removal  or 
deposit  of  mass  occurs  only  on  |  £2  | .  The  Somigliana  dislocations  defined  by  Eq.  (5.20)  are  of 
prismatic  type. 

The  total  energy  change  F  is  minimized  when 

G[a]  =  —  j<7jjnjanjds+  JaynjOuqds-X  J  ads  (5.22) 

2  Ini  fol  Ini 

reaches  a  minimum.  Taking  the  variation  of  Eq.  (5.22)  with  respect  to  a,  we  find  that  G  is 
minimized  when 
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Since  the  Lagrange  multiplier  X  is  constant  in  this  case,  it  is  clear  that  the  relaxation  is  completed 
when  the  total  normal  traction  on  |  &  |  becomes  constant.  In  other  words,  if  the  normal  traction  is 
position  dependent,  a  mass  transport  occurs  further  to  decrease  the  energy.  As  formulated  by 
Koeller  and  Raj  [3],  we  may  state  that  the  gradient  of  chemical  potential  of  an  atom  along  the 
interface  drives  diffusion.  As  the  above  analysis  indicates,  the  chemical  potential  should  be  taken 
proportion  to  the  normal  traction,  regardless  of  hoop  stresses.  Although  it  appears  redundant,  we 
may  add  that  the  tangential  traction  still  remains  and  its  magnitude  is  position  dependent.  The 
traction  on  |  Q  |  after  the  different  relaxation  is  schematically  shown  in  Fig.  5.1  (c). 

Example  4:  a  spherical  inclusion  in  an  isotropic  body 

The  normal  traction  NA  due  to  Cy  is  given  by 

NA=Gj^nkn^  on  |  Q.  \ .  (5.24) 

This  encourages  us  to  postulate  that  bj  is  written  in  a  form  of 

t>i=Bfa?xkW  (5.25) 

which  satisfies  the  conditions  Eq.  (5.20).  We  also  anticipate  that  is  proportional  to  the 
deviatoric  part  of  Oy\  since  the  hydrostatic  stress  is  uniform  on  |  Q I  and  does  not  induce  a  mass 
transport.  The  above  anticipation  means  that 


Bk)c  =  0  (5.26) 

Note  that  Eq.  (5.26)  with  Eq.  (5.25)  assures  the  condition  Eq.  (5.21).  By  conducting  calculation 
similar  to  that  involved  in  Example  3,  the  stress  CTjj  caused  by  Eq.  (5.25)  is  calculated  as 


_  148pv  . 

°1J_35^  U  W  k  ‘ 

+  35^  v)  {~2(7  - 4v)BijX2  -  12v(B^x^Xj  +  B^xjXj) 
+2(7  -  10v)8ijBwxkx*  + 1(5  -  v)a2Bjj 


in  Q. 


If  Bij  is  taken  as 


105(1 -v)  /A 
,J  4|i(35  +  llv)a2  1J’ 


(5.27) 


(5.28) 


the  total  normal  traction  becomes  independent  of  the  position,  that  is 

aijnjnj  +  ay  njnj  =-0^/3  on|Q|.  (5.29) 

Example  5:  a  circular  cylinder  in  an  isotropic  body 

'  Similar  to  Example  4,  we  postulate  that  bi  has  a  form 

bi=Bk^xkx^xi  (i,  k,  1  =  1,2).  (5.30) 

By  conducting  a  similar  calculation  to  that  having  led  to  Eq.  (5.27),  it  can  be  shown  that  the  stress 
ojj  due  to  Eq.  (5.30)  is  given  by 

°ii  =  2fl^){35ijBkWf  +  (2a2'3x2)Blj}  ‘na  <531) 

The  position  independent  condition  of  the  total  normal  traction  on  |  Q  |  gives 

Bij  =-^-^  'ay-  (5.32) 

Interestingly,  this  is  exactly  the  same  as  Eq.  (5.18).  Combining  Eqs.  (5.31)  and  (5.32),  we  have 
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(5.33) 


\ 

'Gy  in  Q., 

J 

after  the  complete  relaxation.  The  total  normal  traction  is  now 

Gynjnj  +  Gy  njnj  =  0^/2  on|Q|.  (5.34) 

Eq.  (5.32)  together  with  Eq.  (5.30)  satisfies  Eqs.  (5.20)  and  (5.21). 

5.2.4  Combination  of  interfacial  sliding  and  local  diffusion 

The  total  energy  change  is  given  by  Eq.  (5.11).  However,  this  case  requires  the  subsidiary 
condition  for  bi  written  as 

—  J  nj  bids  =  0 .  (5.35) 

Ini 

Thus,  our  problem  is  to  minimize 

G[b]  =  —  |ai:njbids+  jay  njbjds-A,  J njbjds.  (5.36) 

2ld  Id  Id 

For  G  to  be  minimum,  the  condition 

Gynj  +Oy  nj  =  )uij  on  |q|  (5.37) 

must  be  satisfied.  Here  X  is  constant,  that  is,  the  stress  is  completely  relaxed  when  the  total  traction 
on  |  Q  |  becomes  constant  and  normal  to  |  Q I ,  as  shown  in  Fig.  5.1  (a).  The  total  stress  in  an 
inclusion  is  not  necessarily  uniform  nor  hydrostatic.  To  emphasize  this  point,  we  drew  Fig.  5.1  (d) 
differendy  from  Fig.  5.1  (a),  where  the  uniform  hydrostatic  state  is  schematically  shown.  The  total 
stress  in  Fig.  5.1  (d)  becomes  uniformly  hydrostatic  at  any  point  within  Q  for  an  elastically 
uniform  inclusion.  This  particular  case  was  noted  by  Eshelby  [24]. -He  mentioned  that  for  an 
ellipsoidal  inclusion  Q  with  an  incoherent  interface,  having  the  transformation  strain  (eigenstrain) 

£^n ,  the  total  stress  Gy  in  the  inclusion  is  found  by  solving 


2a2°ij 


9  ^ij  ^k£xkx<? 
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(5.38) 


^ij  “  ^ijk^(^k£mnemn  —  £k^) 


=  P5ij 


in  Q 


with 


eE=eiSk.  (5-39) 

Here,  Cy^  is  an  elastic  constant  and  S^mn  the  Eshelby  tensor  [25].  The  reasoning  for  Eqs. 
(5.38)  and  (5.39)  is  as  follows: 

The  original  e^n  >  in  general,  results  in  a  stress  Gy  which  is  not  hydrostatic  but  uniform  in 
Q.  Thus,  normal  and  tangential  tractions  due  to  Gy  are  not  uniform  on  |  Q  | .  This  drives  sliding 
and  mass  transport  on  |  Q.  \ .  (Note  that  sliding  and  diffusion  can  occur  quite  easily  on  an 
incoherent  interface.)  The  Somigliana  dislocations,  bi,  satisfy  the  condition  Eq.  (5.35).  After  Gy 
is  completely  relaxed,  the  total  traction  on  |  Q  |  becomes  uniform  and  normal  to  |  Q  | .  Since  only 
an  elastically  uniform  inclusion  is  considered  here,  the  stress  in  the  inclusion  becomes  hydrostatic, 
Pdij,  see  Eq.  (5.38).  bi  can  be  extended  to  the  interior  of  £2,  as  shown  by  Asaro  [22).  As  will  be 
shown  later,  the  dilatation  of  the  eigenstrain  £y  which  gives  the  identical  stress  to  that  due  to  bi  and 
is  defined  by  the  negative  of  the  gradient  of  bi  inside  Q.  is  zero.  Thus,  the  dilatation  of  the  total 
eigenstrain,  £y  +  £y,  is  equal  to  that  of  £y ,  see  Eq.  (5.39). 

Example  6:  a  spherical  inclusion  in  an  isotropic  body 

Since  the  hydrostatic  component  of  Gy  causes  a  uniform  normal  traction  and  no  tangential 
traction  on  |  Q  | ,  we  will  only  consider  the  case  of  Gj2  =  G0  .Apparently  we  can  take 

£  || 
bj=-£X2,  b2=-exi,  other  components  are  zero,  on  |Q|.  (5.40) 

As  will  be  shown  later,  the  stress  oy  due  to  Eq.  (5.40)  can  be  calculated  by  giving  Q  the 
* 

eigenstrain  £y 


Gy  is  now 


£j2  =  £*,  other  components  are  zero. 


(5.41) 
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°12 


2|i(7-5v) 
15(1- v) 


* 

e  , 


other  components  are  zero,  in  Q. 


Thus,  if  £*  is  taken  to  satisfy 


(5.42) 


* 

£ 


15(1 -v) 
2ji(7  -5v) 


"O’ 


in  Q 


(5.43) 


the  total  stress,  Oy  +  Gy,  becomes  zero  in  Q,  a  special  case  of  the  present  subsection.  Note  that 
Eq.  (5.40)  satisfies  Eq.  (5.35)  and  that  bi  is  not  always  normal  nor  parallel  to  I  (2 1 .  This  requires 
that  bi  must  be  achieved  by  the  combination  of  mass  transport  and  sliding. 

What  Mori,  Okabe,  and  Mura  [4]  and  Onaka,  Okada,  and  Kato  [8]  considered  as  the  stress 
relaxation  around  an  ellipsoidal  inclusion  is  not  solely  due  to  interfacial  diffusion  but  also 
accompanies  sliding,  correction  of  Mori's  earlier  paper  [4]. 

5.2.5  The  combination  of  sliding  on  |  Q  |  and  mass  transport  to  I  Q  I  from  an 
outer  free  surface 

This  case  also  holds  for  the  case  that  vacancies  flow  to  I Q I  from  an  outer  free  surface  of  the 
body  D.  No  subsidiary  condition  is  required  when  we  minimize  Eq.  (5.1 1).  It  is  apparent  that  the 
total  energy  is  minimized  when 


Oynj  +  Gy nj  =0  on  |q|.  (5.44) 

That  is,  when  the  inclusion  under  consideration  becomes  free  from  the  traction  on  its  interface.  For 
all  practical  purposes,  one  can  say  that  the  inclusion  becomes  unstressed.  This  is  the  case  studied 
by  Srolovitz  et  al.  for  a  spherical  inclusion  [5].  The  stress  free  state,  the  condition  of  the  most 
effective  relaxation,  holds  for  an  inclusion  of  any  shape.  Note  that  interfacial  sliding  must  occur  in 
addition  to  diffusion. 

If  we  impose  the  subsidiary  condition  of 

bj  =  ani  on  |  Q  | ,  (5.45) 

by  excluding  sliding,  the  minimum  of  Eq.  (5.1 1)  is  achieved  when 
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(5.46) 


.Gjjiijnj  +  Gy  n jiij  =0  on  |  ^2 1 . 

This  is  the  special  case  of  subsection  2.3.  However,  the  total  normal  traction  is  zero,  but  the 
tangential  traction,  in  general,  remains  to  be  further  relaxed. 


5.3.  Discussion 


The  results  of  the  variational  analysis,  Eqs.  (5.5),  (5.14),  (5.23),  (5.37),  (5.44),  and  (5.46), 
are  necessary  conditions  to  describe  an  energy  minimum  state.  Adopting  the  line  used  by  Huang 
and  Mura  [26],  we  can  see  that  these  conditions  are  also  sufficient  ones. 

Let  us  take  interface  related  phenomena.  Let  bj  and  Gy  be  the  Burgers  vectors  of  the 
Somigliana  dislocations  and  those  stresses  due  to  them  which  satisfy  the  necessary  conditions. 
Suppose  that  Sbi  is  introduced,  resulting  in  changes  8Gy  in  the  stresses.  Then,  the  change  in  the 
total  potential  energy  is  expressed  as 


F[b+8b]-F[b] 


=  \  J(Gij+8Gij)nj(bi+5bi)ds+  Jg£ n^bj  +  8bj)ds-i  ^n^ds-  jGynjbjds 

|n|  Ini  ln|  Ini 

=  —  | Gyn jSbjds  -i-  —  | SG^n jbjds  +  —  8Gynj8bjds  +  |g^  njSbjds  . 

2  Ini  2  Ini  2  Ini  Ini 


(5.47) 


Because  of  the  Maxwell-Betti  reciprocal  theorem 

jGynj8bjds=  |  SGyiijbids,  (5.48) 

ln|  Ini 

Eq.  (5.47)  is  rewritten  as 


F[b+Sb]-F[b] 

=  |(Gy  +  Gy  )nj8bids  + 

Ini 


7  ^VjSb.ds 


Ini 


(5.49) 


Because  of  the  conditions,  Eqs.  (5.14),  (5.23),  (5.37),  (5.44),  and  (5.46),  together  with  the 
corresponding  subsidiary  conditions  Eqs.  (5.12),  (5.20),  (5.21),  (5.35),  and  (5.45),  the  first  term 
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in  Eq.  (5.49)  vanishes.  The  second  term  is  the  elastic  energy  of  the  Somigliana  dislocations  of  the 
Burgers  vector  5bj.  It  must  be  positive  for  non-vanishing  5bj.  Thus, 

F[b+8b]  >  F[b],  (5.50) 


which  proves  that  the  necessary  conditions  are  also  the  sufficient  ones.  Almost  identical  discussion 
can  be  given  to  the  case  of  plastic  flow. 

The  energy  expressions,  Eqs.  (5.1)  and  (5.11),  do  not  contain  an  energy  dissipation  term. 
The  energy  dissipation  in  plastic  flow  depends  on  the  path  or  mode  of  deformation.  It  is  caused  by 
non-zero  resistance  to  plastic  flow.  Thus,  the  direct  application  of  Eq.  (5.5)  is  only  valid  for  creep 
deformation,  since  no  threshold  stress  exists  in  creep.  Similarly,  the  interface  related  relaxations, 
sliding  and  diffusion,  occur  with  no  threshold  stress.  Thus,  the  resulting  expressions,  Eqs.  (5.14), 
(5.23),  (5.37),  (5.44),  and  (5.46),  are  always  valid.  Diffusion  is  caused  by  a  gradient  of  a 
chemical  potential  of  an  atom;  the  chemical  potential  is  proportional  to  stresses.  Thus,  diffusion 
occurs  if  a  stress  is  position  dependent  even  when  its  magnitude  is  small.  Similarly,  interfacial 
sliding  occurs  proportional  to  the  tangential  traction  on  an  interface,  as  evidenced  experimentally  by 
grain  boundary  sliding  [27]  and  analyzed  theoretically  [28]. 

Now,  the  strength  of  relaxations  due  to  specific  modes  examined  in  the  present  paper  is 
discussed.  The  strength  here  means  the  effectiveness  of  the  relaxation  modes:  the  larger  the  energy 
decrease,  the  more  effective  the  process.  The  sliding  and  diffusion  processes  cannot  be  compared, 
since  the  conditions  for  the  two  processes,  Eqs.  (5.12),  (5.20)  and  (5.21),  are  independent  of  each 
other.  Mathematically,  this  is  expressed  by  the  fact  that  Eq.  (5.12)  is  not  a  part  of  Eqs.  (5.20)  and 
(5.21),  or  (5.20)  and  (5.21)  is  not  of  Eq.  (5.12).  On  the  other  hand,  the  combination  of  sliding 
and  diffusion  (S  u  D)  has  the  condition  Eq.  (5.35)  which  includes  the  condition  Eq.  (5.12)  for 
sliding  (S)  and  those  Eqs.  (5.20)  and  (5.21)  for  diffusion  (D).  Thus,  we  have 


(5.51) 

Fmin  (S  U  D)^  Fmin  (D), 

(5.52) 

that  is,  the  combined  sliding  and  diffusion  process  is  more  effective  in  relaxation  than  the 
separately  occurring  sliding  and  diffusion  processes.  In  this  respect,  one  should  recall  that  the 
climb  and  annihilation  process  of  an  Orowan  loop  (shear  loop)  also  accompanies  diffusion  (pipe- 
diffusion)  and  glide  motion  [29].  This  is  the  process  to  relax  internal  stresses  developed  by  plastic 
deformation  in  an  inclusion  bearing  material.  When  all  the  loops  are  annihilated,  the  total  energy 
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becomes  zero;  the  lowest  energy  state.  Similarly,  if  bulk  diffusion  from  or  to  a  normal  traction  free 
surface  is  allowed,  we  have  a  more  effective  way  of  relaxation.  This  is  because  no  condition  for 
mass  conservation  such  as  Eqs.  (5.21)  and  (5.35)  is  imposed  in  a  process  allowing  bulk  diffusion. 
It  is  clear  that  bulk  diffusion  with  sliding  reduces  energies  more  effectively  than  that  without 
sliding. 

The  combined  process  of  sliding  and  diffusion  can  be  compared  with  plastic  flow.  As  shown 
by  Asaro,[22]  the  stress  caused  by  bj  on  the  interface  |  Q  |  of  an  inclusion  Q.  can  be  reproduced  by 
giving  the  eigenstrain 

£*j  =  — (b|  j  +  bj  i )  /  2  inQ  (5.53) 

*  *  * 

to  Q.  If  we  understand  £jj  such  that  £y  is  produced  by  plastic  flow  (£jj  =  0),  we  have 

-bi,i  =  0  inD.  (5.54) 

Integrating  Eq.  (5.54)  in  Q  and  applying  the  Gauss-Green  theorem,  we  see 

-Jbinids  =  0,  (5.55) 

Ini 

which  is  identical  to  Eq.  (5.35).  That  is,  the  condition  for  plastic  flow  (P)  is  included  in  the 
condition  for  the  combination  of  sliding  and  diffusion.  On  the  other  hand,  the  condition  Eq.  (5.35) 
does  not  necessarily  lead  to  that  in  Eq.  (5.2)  for  every  point  in  Q.  Thus, 

Fmin  (SuD)<Fmin(P).  (5.56) 

In  short,  the  most  effective  process  to  achieve  the  stress  relaxation  is  the  one  allowing  bulk 
diffusion  and- interfacial  sliding  among  all  the  processes  examined  in  the  present  study. 

However,  what  occurs  in  an  actual  situation  also  depends  on  kinetics.  Bulk  diffusion  takes  a 
long  time  and  can  be  excluded  at  low  temperatures  when  interfacial  sliding  and  diffusion  can  occur 
with  a  sufficient  rate.  The  kinetics  cannot  be  discussed  in  a  general  manner  even  when  the 
relaxation  process  is  limited  to  interfacial  phenomena,  the  kinetics  depends  on  the  shape  of  an 
inclusion  [8].  Confining  to  interfacial  phenomena  and  separately  examining  sliding  and  diffusion, 
one  can  see  from  other  studies  [2-5,  8]  that  the  relaxation  occurs  in  the  first  order  kinetics.  For  a 
round-shaped  or  cylindrical  inclusion,  the  relaxation  time  is  proportional  to  the  cube  of  the 
inclusion  radius  for  interfacial  diffusion  (Xq)  and  to  the  radius  for  interfacial  sliding  (x§).  Thus, 
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above  a  critical  radius,  the  rate-controlling  process  in  the  interfacial  relaxation  is  diffusion.  We  can 
roughly  estimate  the  critical  radius.  Ignoring  the  detailed  geometrical  factor,  the  relaxation  times  are 
given  as  [2-4] 


_  kT  a3 
tD_  D  h(iQ 


(5.57) 


Ts  = 


an 

M-h 


(5.58) 


Here,  D  is  the  interfacial  diffusivity,  a  the  radius,  h  the  thickness  of  the  interface  layer,  Q.  the 
atomic  volume,  h  the  interface  viscosity  and  kT  has  the  usual  meaning.  For  a  disordered  boundary 
such  as  an  incoherent  inclusion/matrix  interface  Ashby  [28]  estimated 


kT 

8bD 


(5.59) 


where  b  is  the  atomic  size.  Eqs.  (5.58)  and  (5.59)  result  in 

,  .kT  a 
Xs  "  D  8|J.bh 


(5.60) 


Since  h  =  b  and  Q  =  b3,  Eqs.  (5.57)  and  (5.60)  lead  to  the  critical  radius  of  the  order  of  the 
atomic  size.  Thus,  for  an  ordinary  inclusion  to  be  discerned  as  a  second  phase,  sliding  occurs  first, 
followed  by  interfacial  diffusion.  It  was  in  fact  evidenced  by  several  experiments  [30-32]  with  the 
observations  that  the  annihilation  rate  of  internal  stresses  in  Fe-particles  embedded  in  Cu  is 
inversely  proportional  to  the  cube  of  the  particle  radius.  Of  course,  the  annihilation  must  be 
accompanied  by  sliding.  This  interfacial  diffusion  will  also  accompany  sliding,  occurring  fast,  to 
achieve  the  condition,  Eq.  (5.37).  The  coexistence  of  sliding  and  diffusion  is  necessary,  similar  to 
diffusional  creep  of  a  polycrystal  which  accompanies  grain  boundary  sliding  [1]. 

Finally,  we  will  briefly  discuss  the  case  that  Gy  is  produced  by  an  external  load.  When  the 

load  is  removed,  the  internal  stress  due  to  £y  on  bi  remains.  This  means  that  the  stress  is  again 

relaxed  by  a  reverse  process.  The  reverse  process  follows  the  path  of  a  forward  process  or  a 
different  path.  In  the  former  case,  a  body  behaves  as  an  anelastic  material.  The  latter  case  is 
realized,  for  example,  when  a  forward  process  is  plastic  deformation  in  a  matrix  surrounding 
inclusions  and  a  reverse  process  is  achieved  by  the  combination  of  local  diffusion  and  sliding  on 
the  interface  of  the  inclusions. 
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5.4.  Summary 


The  stress  relaxation  in  an  inclusion  bearing  material  has  been  examined  from  a  mechanics 
point  of  view  with  the  aid  of  variational  approach.  Four  processes  have  been  dealt  with:  plastic 
flow,  interfacial  sliding,  local  diffusion,  and  bulk  diffusion.  Some  of  them  have  been  combined. 
Each  process  or  combination  involves  a  physical  process  which  restricts  the  plastic  strain  or 
sliding.  The  restriction  involves  a  subsidiary  condition  when  one  seeks  the  minimum  state  of  the 
total  energy  achieved  by  a  particular  relaxation  process.  The  minimum  state  is  described  by  an 
elastic  state  designated  by  the  relaxation  process  under  consideration.  By  examining  the  subsidiary 
conditions  the  effectiveness  of  two  relaxation  processes  can  be  compared. 
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Appendix  1:  Dislocation  punching  from  a  spherical  inclusion  referred  to  Fig.  5.3. 

Because  of  the  spherical  symmetry,  we  can  assume  that  the  displacement,  uj ,  the  derivative 
of  which  gives  the  plastic  strain  is  expressed  by 

m=—  f(x)  (Al) 

x 

Here  f(x)  is  a  function,  to  be  determined,  of  the  distance  from  the  center  to  the  point  x, 
x  =  -s/xpt” .  e*i  =  u*  j  =  0  gives 

-  — -2-  =  0  (A2) 

dx  x 


This  is  solved  as 

f(x)  =  4- 

xz 


(A3) 
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with  a  constant  c.  Inserting  Eq.  (A3)  into  Eq.  (Al)  and  differentiating  tq  ,  we  have  £y 

*  '  8;;  X;X;  A 

eij  =  c  — 5 —  3 — 5^  in  (A4) 

Vx  x 

The  eigenstrains  inside  Q0  and  in  Q  jump  on  |  Q0  |  =  |  Q  | .  Thus,  surface  dislocations  are 
defined  [33].  The  vanishing  condition  for  the  surface  dislocations  on  |  £2 1  gives 

c  =  eT/a3  .  (A5) 

In  this  calculation,  we  excluded  impotent  dislocations  [33],  since  they  produce  no  elastic  field. 

Appendix  2:  Spherical  sliding  inclusion. 

The  stress  field  due  to  sliding  bj  on  the  interface  of  a  spherical  inclusion  in  an  isotropic  body. 
The  hydrostatic  component  of  Gy1  does  not  induce  the  tangential  traction  TA  on  |  Q  | .  Thus,  only 

the  deviatoric  part  'ey  is  included  in  the  calculation.  TA  is  written  as 

TiA  =  'Oy  n  j  -  'a&nkn^ni ,  (A6) 


The  stress  outside  Q  is  very  complicated  and  is  not  recorded  here. 
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The  tangential  traction  Tj  due  to  Gy  on  |  Q  |  is 

Ti  =  aijnj-ak^nkn£ni- 


(A9) 


The  condition 


'  Tj  +  TjA  =  0 


(A10) 


gives 


35(1 -v)  A 
1J  2(i(21 +  v)a2  1J' 


(All) 


106 


6.  Analytical  Study  IV:  Dislocation  Punching  in  Functionally-Graded  Materials 


6.1  Introduction 

As  a  metal/ceramic  composite  is  subjected  to  temperature  change  AT,  a  stress  field  of  large 
magnitude  is  induced  due  to  the  difference  in  the  coefficients  of  thermal  expansion  (CTE)  and  AT. 
Some  relaxation  processes  must  operate  to  reduce  the  large  stresses.  Dislocation  punching  is  one 
of  them  and  has  successfully  been  applied  to  dispersion  hardened  alloys  [1-3]  and  metal  matrix 
composites  [4-9].  Consequences  of  the  punching  have  been  examined  in  terms  of  a  strength 
decrease  or  increase. 

When  a  metal/ceramic  interface  is  planar,  a  similar  relaxation  mechanism  by  dislocation 
punching  from  the  interface  is  expected  to  occur.  Shieu  and  Sass  [10]  studied  both  theoretically 
and  experimentally  the  dislocation  mechanisms  to  relax  thermal  stress  in  a  Pt  film  on  NiO 
substrate.  They  suggested  that  as  a  sample  is  cooled  from  high  temperature  to  room  temperature, 
the  dislocation  punching  takes  place  first  in  a  climb  mode,  then  followed  by  a  glide  mode  process 
during  further  cooling.  Parallel  to  the  dislocation  punching  from  a  thin  metal  coating/ceramic 
substrate  system,  the  generation  of  dislocations  at  a  metal/ceramic  interface  has  also  been  studied  in 
terms  of  critical  thickness  [11,12],  These  studies  were  aimed  at  epitaxial  layers.  Taya  and  Mori 
[13]  recently  examined  a  dislocation  punching  mechanism  from  a  single  interface  of  a  thin  metal 
coating  over  a  thick  ceramic  substrate,  using  a  simple  analytical  model.  They  examined  a  criterion 
for  the  punching  to  occur  and  determined  a  critical  magnitude  of  the  CTE  mismatch  strain,  beyond 
which  the  punching  occurs.  A  thin  metal  coating  over  a  ceramic  substrate  is  widely  used  in 
electronic  packaging,  where  the  microstructure  of  a  thin  metal  coating  bonded  to  a  ceramic 
substrate  can  be  tailored  to  form  a  functionally  graded  material  (FGM),  the  composition  being 
smoothly  graded  from  metal  to  ceramic  phases. 

In  this  report,  we  will  develop  a  new  dislocation  punching  model  for  FGM.  Using 
Eshelby's  inclusion  theory  [14]  and  the  average  field  approach  [15],  the  punching  condition  and 
plastic  strain  in  a  metal  phase  will  be  determined.  A  FGM  layer  consists  of  sublayers  with 
different  compositions  of  metal  and  ceramic.  The  relaxation  of  thermal  stress  in  these  sublayers, 
hereafter  called  simply  layers,  will  be  examined.  The  analysis  is  rather  simple,  since  the 
compositions  of  metal  and  ceramic  phases  changes  one-dimensionally  perpendicular  to  the  interface 
between  the  substrate  and  overlayer.  We  will  determine  the  critical  CTE  mismatch  strain  to  induce 
the  relaxation  as  a  function  of  the  volume  fraction  of  the  ceramic  phase. 
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6.2  Model  and  Analysis 


6.2.1  Elastic  analysis 

A  functionally-graded  material  (FGM)  consists  of  several  thin  coating  layers  and  a  thick 
ceramic  substrate  as  schematically  shown  in  Fig.  6.1.  The  i-th  layer  can  be  approximated  as  a  flat 
inclusion,  parallel  to  the  substrate/overlayer  interface.  The  direction  normal  to  the  interface  is  along 
the  X3-axis.  When  such  a  flat  inclusion  with  uniform  eigenstrains  exists  in  an  infinite  matrix,  all 
the  components  of  the  stress  and  total  strain  in  the  matrix  become  zero  [16].  Here,  the  eigenstrain 
and  total  strain  are  defined  with  respect  to  the  matrix.  Thus,  the  inclusion  causes  no  stress  field 
outside  it.  This  means  that  as  far  as  the  thickness  of  a  substrate  is  sufficiently  large,  the  elastic 
state  of  the  i-th  layer  can  be  examined,  independent  of  the  other  layers.  The  present  model 
assumes  that  the  thickness  of  the  i-th  layer  is  much  smaller  than  that  of  the  substrate.  We  further 
assume  that  the  i-th  layer  has  uniform  eigenstrain  as  a  whole  when  the  overall  elastic  state  is 
evaluated  for  the  layer.  The  overall  elastic  state  is  an  ensemble  of  those  of  the  constituting  phases. 
The  latter  can  be  calculated,  using  the  average  field  theory  [15]. 

Consider  that  FGM  is  subjected  to  a  temperature  change  and  then  a  thermal  mismatch  strain 
relative  to  the  ceramic  substrate  is  generated  in  the  metal  phase  due  to  CTE  mismatch.  There  are 
two  conceivable  structures  in  the  i-th  layer:  ceramic  inclusions  are  surrounded  by  a  metal  matrix 
(Model  A),  Fig.  6.2  (a),  and  metal  inclusions  are  present  in  a  ceramic  matrix  (Model  B),  Fig.  6.2 
(b).  In  Fig.  6.2  inclusions  are  shown  as  isolated  spherical  islands.  The  CTE  mismatch  strain  in 
the  metal  relative  to  the  ceramic  substrate  is  given  by 

eT  =(a  -a  )at8..  =  a*8..  (6.1) 

ij  \  m  c  /  ij  ij 

where  am  and  ac  are  the  CTEs  of  metal  and  ceramic,  respectively,  and  AT  is  the  temperature 
change. 

Suppose  that  the  relaxation  is  achieved  by  uniform  plastic  strain  £p  in  the  metal  phase,  as 
given  by 


sp  =ep[  1  1  -2  0  0  0] 


where  £ 

p 


is  the  magnitude  of  the  plastic  strain  along  the  xi  and  X2  directions. 


(6.2) 
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Figure  6. 1  Schematic  diagram  of  FGM  which  has  multiple  layers  from  a  ceramic  substrate  to  a 
metal  layer. 


Metal:  a*  =(<xm  -ac)AT  Metal:  a*  =(am  -ajAT 


Ceramic  Ceramic 

(a)  (b) 

Figure  6.2  Microstructure  of  the  i-th  layer:  (a)  Model  A  and  (b)  Model  B. 

Now,  we  will  examine  the  elastic  state  of  the  i-th  layer  when  a  temperature  change  is 
induced.  If  thermal  stress  is  large,  relaxation  due  to  plastic  deformation  in  the  metal  phase  occurs. 
This  relaxation  is  equal  to  dislocation  punching  from  the  interfaces  between  FGM  layers.  This 
problem  can  be  analyzed  by  conducting  the  following  hypothetical  procedure,  using  Model  A. 
This  hypothetical  procedure  is  employed  to  apply  the  average  field  theory  [15].  Model  B  is 
similarly  dealt  with,  by  changing  the  roles  of  elastic  constants  for  the  matrix  and  inclusions.  First, 
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the  i-th  layer  is  hypothetically  removed  from  FGM.  Next,  a  temperature  change  is  introduced, 
inducing  eigenstrains  in  the  metal,  Eq.  (6.1).  Further,  plastic  strains,  Eq.  (6.2),  whose  magnitude 
is  to  be  determined  later,  are  introduced  into  the  metal  matrix.  Now,  we  will  apply  the  Eshelby 
theory  [14]  to  the  isolated  i-th  layer.  Since  the  eigenstrains,  Eqs.  (6.1)  and  (6.2),  are  defined  in 
the  matrix  (metal),  we  will  superimpose  -eT  and  -ep  over  the  whole  domain  of  the  i-th  layer  so 

that  the  resulting  problem  can  be  reduced  to  the  Eshelby's  elasticity  problem  where  the  eigenstrains 
are  only  given  in  the  ceramic  inclusions.  This  induces  no  additional  stress.  By  using  the 

equivalency  condition  and  assuming  that  the  ceramic  particles  are  ellipsoidal,  the  average  stress 
in  the  ceramic  particles  is  calculated  by  solving 

o  =  C  •( e+e+£T+ep 

~c  c  \  ~ 

=c„,'(5+e‘e~')' 

Here  C  is  the  stiffness  tensor,  e  is  the  total  strain  in  a  single  particle,  e  is  the  average  elastic  strain 
in  the  matrix  domain,  e1*  is  the  equivalent  eigenstrain  of  the  equivalent  inclusion  yet  to  be 

determined,  subscripts  m  and  c  refer  to  the  matrix  metal  and  ceramic  particle,  respectively.  The 
integration  of  the  stress  over  the  whole  i-th  layer  vanishes,  thus  giving 

e=-f(e-ei*)  (6.4) 

The  total  strain  e  is  related  to  e1*  through  Eshelby's  tensor  S  as 

e  =  S-e‘*  (6.5) 

From  Eqs.  (6.3)~(6.5),  the  average  stress  in  the  metal  matrix  due  to  both  the  CTE  mismatch  strain 
and  the  plastic  strain,  om,T  and  Gm,p ,  respectively,  can  be  computed  as 

°n  T  =  °22T  =  a33'T  =  ~2f^a*  (6-6) 

G-’p=G-’p  =  2fnmY11ep,  a£p  =  2fnmY33ep  (6.7) 


(6.3) 
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where  f  and  |im  denote  the  volume  fraction  of  ceramic  particles  and  the  shear  modulus  of  the  metal, 
respectively.  Yu,  Y33  and  ^  are  given  by 


2ik 
Y  -Jim. 

hi  . 


>m+9)(C22+2C12)  +  K+1)(C21+2Cll) 
15(*  “  Vm  )(C1  J22  "  (212<221 ) 


(6.8) 


^33 


4^[-(5Vm+1)(C22+2C12)  +  4(C21+2Cll)‘ 


Vm)(^'l  J22  ^12^21) 


(6.9) 


2(1  +  V  ) 


(6.10) 


where  v  represents  the  Possion  ratio  of  the  matrix  metal. 


2(>-2vm)(a+2p) 

0-vm) 

gyre 

[l  +  v  ) 

1-V  ) 
m/ 

•(3X+2p) 

1  +  3^ 

f(3Xc+2nc)  +  (1-f)| 

(An+20]- 

X  =  X  -X  ,U=U  -U  .  (6.12) 


X  is  Lame's  constant  and  Cn,  C12,  C21  and  C22  are  tabulated  in  Appendix.  In  these  calculations, 

Eqs.  (6.8)~(6.1 1),  we  have  assumed  that  the  ceramic  particles  are  spherical. 

The  elastic  energy  change  AUi  per  unit  volume  of  the  i-th  layer  is  defined  as 

Jr  fo..(e?+eJdv  ,  where  D  and  Q  denote  the  volume  of  the  i-th  layer  and  the  domain  of  the 
2D  q  y  \  ‘J  •  ij  / 

ceramic  particles.  AUi  can  be  computed  from  Eqs.  (6.1),  (6.2),  (6.6)  and  (6.7),  and  is  calculated 
as 


AUj  =  - f(l - f)(im[2 AYe“  +  (2Yu  +  Y33)a  ep -3 


(6.13) 


with  AY=YirY33- 
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Finally,  the  freed  i-th  layer,  hypothetically  detached  in  the  previous  stage,  is  reattached  to 
FGM.  Since  the  i-th  layer  has  a  permanent  strain  due  to  the  plastic  and  CTE  mismatch  strains 
relative  to  the  ceramic  substrate,  the  whole  i-th  layer  should  be  elastically  strained  until  the  total 
strain  in  the  i-th  layer  becomes  equal  to  the  strain  of  the  substrate.  The  average  total  strain  of  the  i- 
th  layer  is  the  Sum  of  the  CTE  mismatch  strain,  plastic  strain  and  eignestrain  accounting  for 
different  elastic  moduli  of  the  i-th  layer  and  can  be  expressed  as 

e*a  =8T+£P+f£1*  (6.14) 

ij  ij  ij  ij 

where  ej*  is  defined  in  Eq.  (6.3).  The  first  two  terms  are  to  compensate  the  uniformly 
superimposed  eigenstrains,  -et  and  -ep,  to  establish  Eq.  (6.3).  The  relevant  components  of 
£*a  are  given  by 


en  ~e22  =(1-f)  (l-fPK+d>£p 


(6.15) 


where 


2(C22+2Cj2)  qc 

c„c22-c12c2,)^ 


(6.16) 


Since  the  average  stress  component  a°3  in  the  i-th  layer  is  zero  and  the  elastic  strain  needed 

for  the  reattachment  is  the  negative  of  Eq.  (6.15),  the  stress  of  the  i-th  layer  due  to  this  permanent 
strain  is  computed  as  by  using  the  Hooke’s  law 


o 


O 

11 


(6.17) 


where  Ej  and  V;  are  Young's  modulus  and  the  Poisson  ratio  of  the  i-th  layer,  respectively  and  are 
computed  by  the  average  field  method  combined  with  the  Eshelby  theory  [14,  15].  This  straining 
process  increases  the  strain  energy  density  of  the  i-th  layer  by  AU2,  which  is  defined  as  --jG°  £*a 

and  computed  by  Eqs.  (6.15)  and  (6.17)  as 


E. 


au2  =irtr[(1_f)(1-fp)a*+<I)ep 


(6.18) 
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The  change  in  the  total  elastic  energy  density  of  the  i-th  layer,  AU,  subjected  to  the 
temperature  change  with  its  metal  phase  undergone  plastic  deformation,  is  the  sum  of  Eqs.  (6.13) 
and  (6.18)  as  given  by 


AU  =  -f(l-f)|i 
E 


+ 


1  —  v. 


2Ayep.+  (2yu  +Y33)a*8p-3£(a*)~ 

(l-f)(l-fP)a*  +  <De 


(6.19) 


6.2.2  Plastic  work 

The  metal  matrix  is  assumed  to  be  elastic-perfectly  plastic,  while  ceramic  particles  deform 
only  elastically.  The  work  dissipated  by  the  plastic  deformation  of  the  metal  matrix  per  unit 
volume  of  the  i-th  layer  can  be  expressed  as 


W  =  4k(l-f)ep 

where  k  is  the  yield  stress  of  the  metal  matrix  in  shear  deformation. 


(6.20) 


6.2.3  Critical  CTE  mismatch  strain  for  relaxation 


The  minimum  CTE  mismatch  strain  which  can  initiate  the  dislocation  punching  (or  plastic 
relaxation)  is  defined  as  the  critical  CTE  mismatch  strain.  If  the  reduction  in  the  elastic  energy  of 
the  i-th  layer  due  to  the  plastic  deformation  is  larger  than  the  work  required  for  the  plastic 
deformation  of  the  metal  matrix  of  the  i-th  layer,  the  dislocation  punching  takes  place 
( AU  +  W  <  0).  Since  thermal  activation  is  ignored  in  the  present  case,  the  following  condition, 


5(AU) 


9e 


aw 

de 


£P=p 


<  0, 


(6.21) 


'VO 


in  addition  to  AU  +  W  <  0,  must  be  met  for  the  relaxation  to  occur.  Eq.  (6.21)  is  satisfied  when 


*  ^  *c 

a  <  a 


(6.22) 


113 


It  is  noted  in  Eq.  (6.22)  that  since  a*,  and  a*c  are  both  negative,  Eq.  (6.22)  can  be  rewritten  as 


la*  >  a*c 


where  the  critical  CTE  mismatch  strain,  a  ,  is  given  by 


a*c  = 


4k 


2E. 


(6.22)’ 


(6.23) 


6.2.4  Plastic  strain  in  the  metal 

When  the  absolute  magnitude  of  the  actual  CTE  mismatch  strain  is  larger  than  that  of  the 
critical  one,  the  dislocation  punching  takes  place  in  the  metal  domain  of  the  i-th  layer.  The 
magnitude  of  plastic  strain  is  determined  by  minimizing  the  sum  of  the  elastic  strain  energy  change 
and  plastic  work  with  respect  to  the  plastic  strain,  as  expressed  by 

0  (6.24) 

3e  de 

p  p 

Inserting  Eqs.  (6.19)  and  (6.20)  into  Eq.  (6.24),  the  plastic  strain  can  be  computed  as 


,  v  2E.  , 

(l-f). 

+ 

1 

_ j 

i 

a'  r 

j 

-  2  £ 

-4  f  ( 1  —  f )  q  Ay  + - Ld>2 

v  "in  1  —  V. 


6,2.5  The  average  stress  in  the  metal  and  ceramic 


(6.25) 


The  average  stresses  both  in  the  metal  and  ceramic  phases  can  be  given  by  combining  the 
stresses  due  to  the  CTE  mismatch  strain  and  plastic  strain  in  the  hypothetically  free  state  and  those 
due  to  the  reattachment  process.  The  stresses  by  the  first  two  sources  are  given  in  Eqs.  (6.6)  and 
(6.7).  The  stresses  in  ceramic  particles  by  the  last  origin  can  be  computed  by  Eshelby's  equivalent 
inclusion  method  as 
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(6.26) 


o,  =Cc-(e°  +e+e) 

=  Cm-(e° +e+e-e^*j, 

where  e°  =  C“l  •  a0  and  the  components  of  a0  are  defined  in  Eq.  (6.17).  By  solving  Eq.  (6.26) 
with  the  same  steps  as  Eqs.  (6.4)  and  (6.5),  the  stresses  in  the  metal  due  to  o°  can  be  expressed 
as 


°no=°so =2kk°°iK’ 


(6.27) 


where 


[-K  +9)(D,C  -  D  C12H5V„  ;,c,,  -  D.C,,)]  +  1  +  Vm  (6  28) 

Vm)(CllC22  C12C2l) 


P33  f 


-2(!  +  5VJ(D1C22  -  D2C12)  +  8(DlC21-D2C»l)] 


15(*  Vm)((ni<222  ^12^21) 


(6.29) 


and  where  D  is  tabulated  in  Appendix.  The  total  stress  in  the  metal  can  be  computed,  using  Eqs. 
(6.6),  (6.7)  and  (6.27). 


6.3  Discussion 

In  order  to  demonstrate  the  effectiveness  of  the  above  analysis,  we  will  discuss  here  the 
dislocation  punching,  equivalently  plastic  relaxation,  using  a  model  system  of  palladium  (Pd) 
coatings  on  an  AEO3  substrate.  The  thermomechanical  properties  of  Pd  and  AI2O3  are  given  in 
Table  6.1.  We  consider  that  Pd/AEC^  is  cooled  from  270  °C  to  room  temperature  (20°C),  thus 
AT=  -250  °C. 
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Table  6. 1  Material  properties  for  computation 


Pd 

AI2O3 

Young's  Modulus  (GPa) 

112 

363 

Poisson  ratio 

0.39 

0.25 

CTE  (xlO-6) 

11.6 

5.4 

Yield  stress  in  shear  (MPa) 

26 

- 

Temperature  difference  (AT) 

-250 

6.3.1  Critical  CTE  mismatch  strain 


The  critical  CTE  mismatch  strain  in  Model  A  is  shown  against  the  volume  fraction  of  ceramic 
particles  (f)  in  Fig.  6.3.  As  the  volume  fraction  of  the  ceramic  particles  increases,  the  absolute 


magnitude  of  the  critical  CTE  mismatch  strain  increases.  For  smaller  fs,  a 


a 


*c 


,  thus  the 


dislocation  punching  takes  place,  resulting  in  the  plastic  deformation  in  the  metal  phase.  As  f 
increases  beyond  the  intersecting  point  with  |a*|  =  1.55  x  1CT3,  shown  as  dash-dot  line  in  Fig.  6.3, 

the  dislocation  punching  can  not  take  place,  where  the  intersecting  points  of  the  two  critical  strains 
with  the  horizontal  line  (f*)  are  0.87  for  Model  A  and  0.81  for  Model  B.  If  a*  changes,  f* 
changes  according  to  Eqs.  (6.22)  and  (6.23). 

Fig.  6.3  also  shows  the  critical  CTE  strain  cx*c  computed  in  Model  B.  It  is  smaller  than  that 
in  Model  A.  This  result  is  not  surprising,  since  the  structure  differs  in  the  two  models.  We  now 
discuss  why  oc*c  is  smaller  in  Model  B  than  in  Model  A.  We  believe  that  the  reason  for  it  is 
similar  to  the  difference  between  the  upper  and  lower  bounds  in  the  Hashin-Shtrikman  theory  of 
the  elastic  constants  of  a  two  phase  material  [17].  As  proved  by  Wakashima  [18],  the  average  field 
theory  with  spherical  inhomogeneities  gives  the  overall  modulus,  exactly  coinciding  with  the  lower 
bound  in  the  Hashin-Shtrikman  theory,  when  the  inhomogeneities  are  elastically  harder  than  the 
matrix.  When  the  inhomogeneities  are  softer,  the  average  field  theory  gives  the  value,  equal  to  the 
upper  bound  in  the  Hashin-Strikman  theory.  Since  AI2O3  is  elastically  harder  than  Pd,  Model  A 
gives  a  result,  corresponding  to  the  lower  bound.  Model  B  corresponds  to  the  upper  bound.  The 
examination  of  an  extreme  case  makes  this  reasoning  clearer.  Suppose  that  a  ceramic  phase  is 
extremely  harder  than  a  metal  phase.  In  Model  B,  then,  the  overall  CTE  strain  is  almost  the  same 
as  that  of  the  substrate  unless  f=0.  The  stress  in  the  spherical  metal  particles  is  hydrostatic.  Then, 
a  large  a*  in  magnitude  is  needed  to  induce  the  plastic  deformation.  On  the  other  hand,  in  Model 
A,  the  overall  CTE  strain  differs  from  that  of  the  substrate.  Because  of  this,  the  stress  to  induce 
the  plastic  deformation  can  exist  when  a*  is  relatively  small  in  magnitude. 
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Volume  fraction  of  ceramic  (f) 

Figure  6.3.  The  prediction  of  the  critical  CTE  mismatch  strain  in  the  Pd-Al203  FGM  by  Models  A 
and  B  as  a  function  of  the  volume  fraction  of  AI2O3.  The  actual  CTE  mismatch  strain  is 
-1.55  xlO-3. 

6.3.2  Plastic  strain  in  the  metal 

The  plastic  strain  in  the  metal  phase  after  the  relaxation  is  also  computed  from  Eq.  (6.25)  and 
the  results  are  shown  as  a  function  of  the  volume  fraction  of  the  ceramic  phase  in  Fig.  6.4.  The 
smaller  the  volume  fraction  of  the  ceramic  phase,  the  larger  plastic  deformation  the  metal  phase 
undergoes.  This  is  because  a  larger  thermal  stress  is  induced  in  the  metal  phase  as  the  amount  of 
the  ceramic  phase  decreases.  This  explanation  also  holds  for  the  calculated  difference  in  the  plastic 
strains  in  Models  A  and  B. 

6.3.3  Stress  state  in  metal 


The  present  analysis  employs  the  energy  approach,  in  which  the  expressions  such  as  Eq. 
(6.19)  contain  no  terms  explicitly  expressed  by  the  stresses  in  the  metal  phase.  This  is 
advantageous  because  an  energy  is  a  scalar  quantity  and  the  parameters  to  determine  are  also 
scalars  such  as  a*  and  e  ,  owing  to  a  certain  symmetry  requirement.  However,  we  can  conduct  a 

similar  analysis,  using  a  stress  approach  by  examining  the  magnitude  of  Von-Mises  stress, 
=  |<5™  -  G™  I.  If  this  is  larger  than  2k,  the  plastic  deformation  occurs  until  it  becomes 
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equal  to  2k,  as  £p  increases.  To  demonstrate  this  point,  we  have  drawn  Fig.  6.5,  indicating  that 
when  f  is  smaller  than  f*,  oj"  -  before  the  relaxation  is  large  than  2k.  After  the  relaxation, 
becomes  2k.  It  is  noted  in  Fig.  5  that  the  intersecting  points  with  lines  are  the  same  as 

those  in  Fig.  3,  indicating  that  the  model  based  on  the  stress  criterion  gives  rise  to  the  same  results 
by  the  energy  approach. 
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Volume  fraction  of  ceramic  (f) 

Figure  6.4.  Plastic  strain  in  Pd  in  Pd-Ab03  FGM  calculated  by  Models  A  and  B  as  a  function  of 
volume  fraction  of  AI2O3. 
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6.3.4  Effect  of  the  heterogeneous  structure  of  a  layer 

To  calculate  the  elastic  state,  the  i-th  layer  is  hypothetically  removed  from  FGM.  After  this 
stage  and  before  the  reattachment  process,  the  traction  on  its  surfaces  is  zero,  regardless  of  the 
presence  of  ej  and  eP.  When  these  are  introduced,  the  originally  planar  surfaces  become  slightly 

roughened,  because  of  the  traction  free  condition.  The  wavelength  of  the  roughness  is  in  the  order 
of  the  particle  spacing.  Thus,  when  the  layer  is  reattached,  stress  due  to  the  surface  roughness  is 
induced,  in  addition  to  Eq.  (6.17).  However,  the  additional  stress  diminishes  to  a  negligible 
degree  at  a  few  wavelengths  below  the  surfaces,,  the  Saint  Venant  principle.  If  the  thickness  of  the 
layer  is  sufficiently  larger  than  the  particle  spacing,  the  effect  of  the  surface  roughness  on  the 
energy  evaluation  can  be  ignored. 

6.3.5  General 

Recently,  extensive  studies  were  reported  for  the  nucleation,  arrangement  and  motion  of 
crystal  dislocations  in  multi-layered  structures  [19-24],  In  these  studies,  the  structures  consist  of 
single  crystal  layers,  which  are  coherendy  connected  to  each  other.  On  the  contrary,  in  the  present 
study  a  layer  is  polycrystalline,  containing  embedded  particles.  Thus,  a  continuum  approach  as 
adopted  here  is  adequate.  In  this  case,  the  process  of  dislocation  generation  is  not  involved, 
because  the  composition  and  structure  do  not  change  significantly  across  the  interface  of  two 
neighboring  layers.  In  particular,  most  of  a  metal  phase  is  connected  at  the  interface  with 
neighboring  layers  in  Model  A.  In  Model  B,  the  relaxation  occurs  inside  isolated  metal  particles, 
whose  interface  with  a  ceramic  matrix  cannot  be  coherent  everywhere.  Thus,  the  ceramic/metal 
interface  can  be  assumed  to  be  dislocated  before  the  relaxation.  Similarly,  no  change  in  interfacial 
energies,  to  be  calculated  essentially  by  the  discrete  dislocation  core  energies,  is  involved. 

Finally,  we  would  like  to  mention  that  if  diffusion  and  sliding  on  the  metal/ceramic  interfaces 
can  occur,  a  more  relaxed  state  is  achieved  [25],  Roughly  speaking,  this  is  because  the  combined 
process  of  sliding  and  diffusion  on  a  ceramic/metal  interface  is  equivalent  to  or  more  effective  than 
the  plastic  deformation  occurring  inside  the  metal.  It  is  noted  that  sliding  occurs  as  long  as  a 
tangential  traction  does  not  vanish,  while  usual  plastic  deformation  encounters  a  certain  resistance, 
yield  strength.  Similarly,  interfacial  diffusion  occurs  whenever  the  gradient  of  a  normal  traction 
along  the  interface  is  not  zero.  However,  these  sliding  and  diffusional  relaxation  processes  occur 
with  a  significant  rate  only  at  high  temperatures. 
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6.4  Summary 


The  criterion  to  initiate  dislocation  punching  in  the  metal  of  FGM  consisting  of  metal  and 
ceramic  phases  was  examined.  Two  models  were  employed:  ceramic  particles  embedded  in  metal 
matrix  (Model  A)  and  metal  particles  embedded  in  ceramic  matrix  (Model  B).  The  theory  was 
applied  to  a  model  Pd-Al203  FGM  system,  where  100%  Pd  is  located  atop  and  graded  smoothly 
toward  100%  AI2O3  ceramic  substrate.  When  the  volume  fraction  of  the  ceramic  is  larger  than 
0.87  and  0.81  for  Models  A  and  B,  respectively,  the  dislocation  punching  does  not  take  place.  For 
the  lower  volume  fraction  of  AI2O3,  the  plastic  strains  by  Models  A  and  B  were  also  computed, 
which  increase  as  the  volume  fraction  of  the  ceramic  decreases.  The  Von-Mises  stresses  in  the 
metal, 


la™  -  a™  ,  before  and  after  the  relaxation  were  calculated,  showing  that  it  is  relaxed  to  the 


yield  stress  of  the  metal  when  the  punching  criterion  is  met.  The  critical  CTE  mismatch  strains, 
plastic  strains  and  stresses  before  and  after  the  relaxation  by  both  Models  A  and  B  coincide  each 
other  relatively  well.  Thus,  the  present  model  can  predict  the  whole  range  of  the  volume  fraction 
of  ceramic  phase. 
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Appendix 


Since  the  calculation  for  Q;  and  D;  is  rather  involved,  only  its  results  are  given  below. 
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7.  Conclusion 


The  major  results  of  this  project  are  summarized  as 

(1)  The  dimensional  change  (e)  of  a  metal  matrix  composite  subjected  to  creep/thermal  cycling 
increases  with  maximum  temperature  (Tmax),  applied  stress  (a)  while  the  volume  fraction  of  filler 
(f)  dependence  of  e  is  that  e  increases  with  f  and  peaked  at  a  certain  value  (f*),  then  e  decreases 
with  f  beyond  f*.  The  shape  of  thermal  cycling,  particularly  duration  of  heating  has  a  strong  effect 
on  the  dimensional  change  of  a  metal  matrix  composite.  For  example,  if  the  duration  is  short  even 
for  higher  maximum  temperature,  its  effect  on  the  composite  dimensional  change  is  much  less  than 
for  the  case  of  longer  duration  of  high  temperature  (Tmax)  exposure  and  lower  Tmax-  The  case  of 
short  duration  of  heating  with  relatively  high  Tmax  can  be  encountered  in  the  application  of  a  MMC 
to  brake  disc,  to  which  the  present  model  is  applicable. 

(2)  The  main  driving  force  for  the  degradation  of  metal  matrix  composites  after  subjected  to 
creep/thermal  cycling  is  the  CTE  misfit  strain  e*=AaAT  where  a,  T  are  the  CTE  and  temperature 
and  A  is  change.  When  e*  is  not  large  so  that  the  local  yielding  around  ceramic  particulates  takes 
place,  the  model  of  dislocation  punching  is  needed  to  simulate  the  relaxation  precisely  while  £*  is 
large,  the  assumption  of  general  yielding  in  the  matrix  phase  is  valid. 

(3)  Relaxation  phenomena  in  a  metal  matrix  composite  can  take  many  forms:  interfacial  sliding, 
interfacial  diffusion,  and  bulk  diffusion.  Complete  relaxation  can  be  found  by  minimizing  the  total 
potential  energy  by  using  the  variational  principle  developed  in  this  project.  For  example,  the 
complete  relaxation  of  a  creeping  composite  is  that  the  Von-Mises  stress  in  the  metal  matrix 
becomes  zero,  i.e.  hydrostatic  state  of  stress  which  can  be  obtained  as  a  part  of  the  present 
analysis. 

(4)  In  a  functionally  graded  material  (FGM)  metal  phase  is  smoothly  changed  from  metal  phase  to 
ceramic  phase  through  a  number  of  layers  consisting  of  metal-ceramic  composite.  The  relation  of 
CTE  misfit  strain  in  a  FGM  by  dislocation  punching  is  studied  with  findings:  there  exists  a 
threshold  FGM  layer,  beyond  which  toward  ceramic  phase,  no  dislocation  punching  takes  place, 
and  below  which  toward  metal  phase,  dislocation  punching  takes  place,  resulting  in  the  plastic 
strains  (ep)  in  the  metal  phase.  The  magnitude  of  ep  is  obtained  by  the  present  model. 
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